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Preface 


This is the abstract sequel booklet monograph to the recently published mono- 
graphs, by the same authors, titled: “Intelligent Numerical Methods: Applications to 
Fractional Calculus,” Studies in Computational Intelligence 624, and “Intelligent 
Numerical Methods II: Applications to Multivariate Fractional Calculus,” Studies in 
Computational Intelligence 649, both in Springer Heidelberg New York, 2016. It is 
regarding applications of Newton-like and other similar methods for solving 
abstract functional equations, which involve abstract Caputo and Canavati type 
fractional derivatives. The functions we are dealing with are Banach space valued 
of a real domain. These are studied for the first time in the literature, and chapters 
are self-contained and can be read independently. In each chapter, the first sections 
are prerequisites for the final section of abstract fractional calculus applications. 
This short monograph is suitable to be used in related graduate classes and research 
projects. We exhibit the maximum of these numerical methods at the abstract 
fractional level. 

The motivation to write this monograph came by the following: Various issues 
related to the modeling and analysis of fractional order systems have gained an 
increased popularity, as witnessed by many books and volumes in Springer’s 
program: 


http://www.springer.com/gp/search ?query=fractional&submit=Prze%C5 %9Blij 


and the purpose of our book is to provide a deeper formal analysis on some issues 
that are relevant to many areas for instance: decision making, complex processes, 
systems modeling and control, and related areas. The above are deeply embedded in 
the fields of engineering, computer science, physics, economics, social and life 
sciences. 

The list of covered topics here follows: 


explicitimplicit methods with applications to Banach space valued functions in 
abstract fractional calculus, 

convergence of iterative methods in abstract fractional calculus, 

equations for Banach space valued functions in fractional vector calculi, 
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iterative methods in abstract fractional calculus, 

semi-local convergence in right abstract fractional calculus, 

algorithmic convergence in abstract g-fractional calculus, 

iterative procedures for solving equations in abstract fractional calculus, 
approximate solutions of equations in abstract g-fractional calculus, 
generating sequences for solving equations in abstract g-fractional calculus, 
and numerical optimization with fractional invexity. 

An extensive list of references is given per chapter. 

This book’s results are expected to find applications in many areas of applied 
mathematics, stochastics, computer science, and engineering. As such, this short 
monograph is suitable for researchers, graduate students, and seminars of the above 
subjects, also to be in all science and engineering libraries. 

The preparation of this book took place during the academic year 2016-2017 in 
Memphis, Tennessee, and Lawton, Oklahoma, USA. 

We would like to thank Prof. Alina Lupas of University of Oradea, Romania, for 
checking and reading the manuscript. 


Memphis, USA George A. Anastassiou 
Lawton, USA Ioannis K. Argyros 
June 2017 
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Chapter 1 

Explicit-Implicit Methods with Applications 
to Banach Space Valued Functions 

in Abstract Fractional Calculus 


Explicit iterative methods have been used extensively to generate a sequence approx- 
imating a solution of an equation on a Banach space setting. However, little attention 
has been given to the study of implicit iterative methods. We present a semi-local con- 
vergence analysis for a some general implicit and explicit iterative methods. Some 
applications are suggested including Banach space valued functions of fractional 
calculus, where all integrals are of Bochner-type. It follows [5]. 


1.1 Introduction 


Sections 1.1—1.3 are prerequisites for Sect. 1.4. 

Let B,, B, stand for Banach spaces and let & stand for an open subset of By. 
Let also U (z, p) := {u € By: ||u —z|| < p} and let U (z, p) stand for the closure of 
U (z, p). 

Numerous problems in Computational Sciences, Engineering, Mathematical 
Chemistry, Mathematical Physics, Mathematical Economics and other disciplines 
can be brought in a form like 

F (x) =0 (1.1.1) 


using Mathematical Modeling [1-16], where F : 82 — Bz is a continuous operator. 
The solution x* of Eq. (1.1.1) is sought in closed form, but this is attainable only 
in special cases. That explains why most solution methods for such equations are 
usually iterative. There is a plethora of iterative methods for solving Eq. (1.1.1). We 
can divide these methods in two categories. 

Explicit Methods [6, 7, 11, 14, 15]: Newton’s method 


Mat Sie FP Ga) Fe). (1.1.2) 
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Secant method: ; 
Xn+1 = Xn — [tet Xn; F| F (Xn) ; (1.1.3) 


where [-, -; F'] denotes a divided difference of order one on Q x Q [7, 14, 15]. 
Newton-like method: 
Xn4+1 = Xn — E,'F (Xn) ey (1.1.4) 


where E,, = E (F) (x,) and E : Q — L (Bj, B>) the space of bounded linear oper- 
ators from B, into By. Other explicit methods can be found in [7, 11, 14, 15] and 
the references there in. 

Implicit Methods [6, 9, 11, 15]: 


F (Xn) + Ay (Xn41 = Xn) =0 (1.1.5) 
Xntl = Xn — A, F (Xn) (1.1.6) 

where A, = A (Xn41,Xn) = A(F) (41, Xn) and A: Q x Q—> L(B,, Bo). 
There is a plethora on local as well as semi-local convergence results for explicit 


methods [1-8], [10-15]. However, the research on the convergence of implicit meth- 
ods has received little attention. Authors, usually consider the fixed point problem 


P(x) =x, (1.1.7) 
where 
P(x) =x+ F(Z) +A(x, 2) (4-2) (1.1.8) 
or 
P(x) =z—A(x,z)7' F(z) (1.1.9) 


for methods (1.1.5) and (1.1.6), respectivelly, where z € Q is given. If P is a con- 
traction operator mapping a closed set into itself, then according to the contraction 
mapping principle [11, 14, 15], P, has a fixed point x* which can be found using the 
method of succesive substitutions or Picard’s method [15] defined for each fixed n 
by 
Yk+1n = Py (Ye,n) > Yon = Xn Xn41 = lim Ykyn+ (1.1.10) 
k—> +00 


Let us also consider the analogous explicit methods 
F (Xn) + A (Xn, Xn) Xn41 — Xn) = 0 (1.1.11) 
Ane = te =A Caste) F GG) C13) 


F (Xn) + A (rn, Xn—1) (Xn+1 = X4) =0 (1.1.13) 
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and 
Xn = Xn — A Kn» Xn-1) F Gn). (1.1.14) 
In this chapter in Sect. 1.2, we study the semi-local convergence of method (1.1.5) 
and method (1.1.6). Section 1.3 contains the semi-local convergence of method 
(1.1.11), method (1.1.12), method (1.1.13) and method (1.1.14). Some applications 


to Abstract Fractional Calculus are suggested in Sect. 1.4 on a certain Banach space 
valued functions, where all the integrals are of Bochner-type [8]. 


1.2 Semi-local Convergence for Implicit Methods 


The semi-local convergence analysis of method (1.1.6) that follows is based on the 
conditions (#): 
(h;) F:QcC B, > By is continuous and A (F) (x, y) € £(B), Bo) for each 
(x,y)EQxX Q, 
(hz) There exist 1 > 0 and Qo C B, such that A (F) (x, yy! € £L(Bo, B,) for 
each (x, y) € Qo x Qo and 
JA) @, yy" se. 


Set Q] = QN Qo. 
(h3) There exist real numbers a1, a2, a3 satisfying 


0 < az < aq; and0 < a3 <1 
such that for each x, y € Qy 
IF @) — F(y)-A(F)@, y)@—y)IlS 
1(S lle — yl +2 ly — xoll +08) Ix — yl. 
(h4) For each x € Qo there exists y € Qo such that 
y=x-A(y, x) FQ). 
(hs) For xo € Qo and x; € Qo satisfying (h4) there exists 7 > 0 such that 
||A CF) Ga, x0)! F (x0) || < 2. 
(he) h = ayn < }(1—a3)’. 


and _ 
(hz) U (xo, t*) C Qo, where 
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be | Ravina ay £0 
rere a, =0. 
Then, we can show the following semi-local convergence result for method (1.1.6) 
under the preceding notation and conditions (#7). 


Theorem 1.1 Suppose that the conditions (H) are satisfied. Then, sequence {xy} 
generated by method (1.1.6) starting at xo € Q is well defined in U (xo, t*), remains 
in U (x0, t*) for each n = 0, 1,2, ... and converges to a solution x* € U (xo, t*) 
of equation F (x) = 0. Moreover, provided that (h3) holds with A(F)(z,y) replacing 
A(F)(x,y) for each z € Q) , if a; 4 0, the equation F (x) = 0 has a unique solution 
x* in U, where 
a a | U Go. t*) N26, ifh = 5 (1 — 03)” 

~ | U Qo, #7) Qo, ifh < 4 (1-3) 


and, if a, = 0, the solution x* is unique in U (xo, 72), where t** = 


1—a3+4/ (1—a3)?—2h 


ay 


Proof Case a; #0. Let us define scalar function g on R by g(t)= of _ 
(1 — a3) t + 7 and majorizing sequence {t,} by 


to = 0, th = t&-1 + g (t%e-1) foreachk = 1,2,.... (1.2.1) 
It follows from (H6) that function g has two positive roots t* and t**, t* < t**, and 
ty < ty41. That is, sequence {t,} converges to 7*. 
(a) Using mathematical induction on k, we shall show that 


Xx41 — Xell S ei — tk. (1.2.2) 


Estimate (1.2.2) holds for k = 0 by (As) and (1.2.1), since ||x; — xo|| < 7 = t — fo. 
Suppose that for 1 <m<k 


Xm — Xm—1ll < tn — tn-1- (1.2.3) 


Them, we get |x, — xoll < tf —t = t% < t* and A (xg, xg_1) is invertible by (h2). 
We can write by method (1.1.6) 


Xep1 — Xe = Ag! (F (xn) — F (e1) — Ant (te — Xe-1))- (1.2.4) 


In view of the induction hypothesis (1.2.3), (h2), (43), (4), (1.2.1) and (1.2.4), we 
get in turn that 


Ica — Xell = Ag! F Gx) |] = Ag! (F Ge) — F Gri) — Agi Ore — xe-1))| 
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< Az! | LF Ge) — F Cet) — Ani Ge — xe S 


a 

aa (> ce — ae-all + 022 Ul xe—1 — x0ll + as) llxe — xe-1l] < (1.2.5) 
Q1 2 
= (th — the)” + Op (ty — th_1) th_-1 + 03 (te — th-1) = 


a 
(te = te1)° + 2 (te = th1) th_1 + 03 (te — th-1) — (te — th-1) + 8 (k-1) = 
& (te) — (1 — 2) he = te-1) te-1 S 
8 (tk) = thei — th, (1.2.6) 


which completes the induction for estimate (1.2.2). 
That is, we have for any k 


WXea1 — Xell S ter — (1.2.7) 


and 
xx — Xol] < th <1". (1.2.8) 


It follows by (1.2.7) and (1.2.8) that {x;,} is a complete sequence in a Banach space 


B, and as such it converges to some x* € U (xo, t*) (since U (xo, t*) is a closed 
set). By letting k — +00, using (h;) and (h2), we get 1 Jim || F (xx) || = 0, so 
> +00 


F (x*) =0. 
Let x** € U be such that F (x**) = 0. We shall show by induction that 
||x** — x, || <t* —% foreachk =0,1,2,.... (1.2.9) 


Estimate (1.2.9) holds for k = 0 by the definition of x** and U, Suppose that 
\|x** — xg|| < t* — t,. Then, as in (1.2.5), we obtain in turn that 


|x** — xegi |] = [[x** — xe + Ag! F (xe) — Ag! F (x*) |] = 


Act (Ac 2" — 2%) + F Gx) — F (x) 


IA 


[Aci T IF @*) - F Ge) — An (2 — x4) | 


IA 


(> ||x** — xe] + a llxe — xoll + as) ||x** — xz|| < 


(; (t* _ tk) + ot, + a3) (t* _ tk) = 
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> (r*)? + S (th)? — ayt,t* + a (t* — tk) te + 03 (t* —&) = 


a 
—n+ Ud — a3) a Zt —aytt® + Qt,t* _ ant; +a3t* — a3zt, 


=f = te, (1.2.10) 
which completes the induction for estimate (1.2.9). Hence, _jim x, = x**, But we 
—> +00 
showed that lim x; = x*, so x** = x*. 
k->+00 


Case a; = 0. Then, we have by (3) that a2 = 0 and estimate (1.2.5) gives 
Ilxe+1 — Xell < 03 lle — xXe-11] < ++ < Of |]x1 — xoll < a47 (1.2.11) 


and 
Wxxa1 — Xoll < Wxeg1 — xell + lle — xe-1 |] +--+ + llr — oll 


_ k+1 
2-3 pe (1.2.12) 
1 = 03 1 = a3 


Then, as in the previous case it follows from (1.2.11) and (1.2.12) that 


l-a , 
3 aby, (1.2.13) 


Ilxx+i — Xxll < 

1- a3 

so sequence {x,} is complete and x* solves equation F (x) = 0. Finally, the unique- 
ness part follows from (1.2.10) for ay = a2 = 0, since 


2k 


|x - xe+1| <ay ||x** — XE | < okt |e" — x0 | < a (1.2.14) 
l—-—ay 


which shows again that lim x, = x**. | 
k>+00 
Remark 1.2 (1) Condition (h2) can be incorporated in (3) as follows 
(h4) There exist real numbers @, @2, @3 satisfying 0 < @2 < a, and0 < a3 < | 
such that for each x, y € Q 


JAG "LF @)-FO)-AGy)@—-y)I] s 
(1/2) ll = yll + lly — x0ll +.) [lx — yl. 
Then, (4) can replace (M2) and (h3) in Theorem 1.1 for a) = @, @ = ao, 
Q3 = a3 and Qo = Q. Moreover, notice that @) < a}, @ < a; and @3 < a3, which 


play a role in the sufficient convergence criterion (/6), error bounds and the precision 
of t* and t**. Condition (43) is of Mysowksii-type [11]. 
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(2) Suppose that there exist J) > 0, ag > 0 and L € £(B,, Bo) with Ee 
L (Bo, By) such that || L~'|| < 19! 


|| A (F) (x, y) — L|| < a4 for each x, y € Q 


and 
as = Iga <i. 


Then, it follows from the Banach lemma on invertible operators [7, 9, 11, 14, 15] 
and 


|Z AG) @ y) - LI < Ig tag = a5 <1 


it 


that A (F) (x, y)! € £L(Bo, By). Set 17“! = Tras” Then, under these replacements, 
condition (42) is implied, so it can be dropped from the conditions (#7). 

(3) Clearly method (1.1.5) converges under the conditions (1), since (1.1.6) 
implies (1.1.5). 

(4) Let R > 0 and define Ro = sup {t € [0, R) : U (x, Ro) C D}. Set Qo = 
U (xo, Ro). Condition (3) can be extended, if the additional term a» |x — xo] 
is inserted inside the paranthesis at the right hand side for some a2 > 0. Then, 
the conclusions of Theorem 1.1 hold in this more general setting, provided that 
a3 = a2 Ro + a3 replaces a3 in conditions (he) and (h7). 

(5) Concerning the solvability of Eq. (1.1.6) (or (1.1.5)), we wanted to leave condi- 
tion (h4) as uncluttered as possible in conditions (H). We did this because in practice 
these equations may be solvable in a way other than using the contraction mapping 
principle already mentioned earlier. 


Next, we show the solvability of method (1.1.5) using a stronger version of the 
contraction mapping principle and based on the conditions (C) : 
(c1) = (1). 
(C2) There exist yo € [0, 1), y, € [0, +00), y2 € [0, 1), x9 € Q such that for each 
x,y,ZzEQ 
|Z +A(x,z) -AQ. DI <0, 


A (x,z) -AQy, Zl < alle — yIl 


_ 72 Ilxo — zll for xo Az 
IF (z) + A (0, 2) Go - DIS iP Cell forsee 
(c3) 

yo + V1 Ilxoll + y2 < 1 for y2 4 0, 


Yo +1 |lxoll < 1 for y2 = 0, 


2 
(rte) 2 Ol say, DO, 


V1 
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yo < | fory, =0 


and _ 
(c4) U (xo, r) € 2, where 


F i= 
IF Gol, 1 Got nol) 5 29, 
1—(% + 7% Ilxoll) V1 
F 
IF Goll <r fory = 0, 
1-0 


1—(% +7 Ilxoll) 
r< 
eal 


for z= x0, y1 4 0. 


Theorem 1.3 Suppose that the conditions (C) are satisfied. Then, for each n = 
0,1,2,... Eg.(1.1.5) is unique solvable. Moreover, if Ap € £L(Bo, B,), then 
Eq. (1.1.6) is also uniquely solvable for eachn = 0, 1, 2, ... 


Proof We base the proof on the contraction mapping principle. Let x, y € U (xo, r). 
Then, using (1.1.8) we have in turn by (c2) that 


|P. (x) — P, (y) I] = IU + A, 2) -AQ, 2) @ — y) -— (AG, 2) - AQ, 2)) ll 
<|Z+AGQ,z)-A(y, ZI lx —yIl +1 AG, 2) -AQ, ZI lz 


So llx — yll +n Cz — oll + IlxolD Ix — yl 


< p(llx — xoll) lla — yll, (1.2.15) 
where 
p= 0 iH . a uu ae — (1.2.16) 
Notice that y (tf) € [0, 1) fort € [0, r] by the choice of r in (cq). 
We also have that 
I| P: ©) — xoll < IP; @&) — Pz Go) ll + IP: Go) — xoll - (1.2.17) 
If z = xo in (1.2.17), then we get by (c3), (ca) and (1.2.15) that 
| Pro 2) — xo] < 9 (lx — xoll) Ile — xoll + IF Go) 
<(+N Ixoldr + IF Goll <r. (1.2.18) 


The existence of x; € U (xo, r) solving (1.1.5) for n = 0 is now established by the 
contraction mapping principle, (1.2.15) and (1.2.18). 
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Moreover, if z 4 xo, the last condition in (c3), (c3) , (c4) and (1.2.17) give instead 
of (1.2.18) that 


I| P; (x) — xoll < ¢ (lx — xoll) Ilx — xoll + 2 llx — xoll 
= (o+%1 Ilxol) + 42)7 = (1.2.19) 


Then, again by (1.2.15), (1.2.19) and the contraction mapping principle, we guarantee 
the unique solvability of Eq.(1.1.5) and the existence of a unique sequence {x,} 
for each n = 0, 1, 2, ... Finally, equation (1.1.6) is also uniquely solvable by the 
preceding proof and the condition A,! € L (Bo, Bi). a 


Remark 1.4 (a) The gamma conditions can be weakened, if 7; are replaced by func- 
tions 7; (t), i = 0,1, 2,3. Then, +; will appear as 7; (||x — xo||) and +4; (7) in the 
conditions (C). 

(b) Sect. 1.2 has an interest independent of Sect. 1.4. However, the results espe- 
cially of Theorem 1.1 can apply in Abstract Fractional Calculus as suggested in 
Sect. 1.4. As an example crucial condition (3) is satisfied in (1.4.8), if we choose 
a2 = a3 = Oandla, = ao where c is defined in (1.4.8). Similar choices can be given 
for the rest of the special cases of (h3) appearing in Sect. 1.4. 


1.3 Semi-local Convergence for Explicit Methods 


Theorem |.1 is general enough so it can be used to study the semi-local convergence of 
method (1.1.11), method (1.1.12), method (1.1.13) and method (1.1.14). In particular, 
for the study of method (1.1.12) (and consequently method (1.1.11)), we use the 
conditions (H' ) : 

(h\) F:QcC B, > B, is continuous and A (F) (x, x) € £(B,, Bo) for each 
xEQ. 

(h5,) There exist / > 0 and Qo C B, such that A (F) (x, x)! € £L (Bo, B)) and 


|A(F)(@, x) ' |] set. 


Set Q) = QN MQ. 
(h,) There exist real numbers 7, a2, y3 satisfying 


O<a.<y and0<73 
such that for each x, y € Qy 
lF@y—FW)=—4AG) OG -—yIS 


V1 
1 (Fle — yl + a2 ly — xoll +78) Ile — yl: 
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(h',) For each x, y € Q; and some y4 > 0, ys => 0 


|A (x,y) -AQ, y)Il Sly 


or 
|A (x,y) — A(y, y)Il < Lys Ilx — yll- 

Set ay = y, +5 and a3 = 73 + 4, if the second inequation holds or a, = 7, and 

a3 = 73 + y4, if the first inequation holds. Further, suppose 0 < a3 < 1. 


(h) There exist x9 € Qo and 7 > O such that A (F) (x0, xo)! € £ (Bo, B,) and 


|| A (F) (x0, x0)! F (xo) || < n- 


(he) = (he) 
and 
(h5) = (h7). 


Then, we can show the following semi-local convergence of method (1.1.12) using 
the conditions (H ’) and the preceding notation. 


Proposition 1.5 Suppose that the conditions (H') are satisfied. Then, sequence {x;} 
generated by method (1.1.12) starting at Xx) € Q is well defined in U (xo, t*), remains 
in U (xo, t*) for each n = 0, 1,2, ... and converges to a solution x* € U (xo, t*) 
of equation F (x) = 0. Moreover, if a; 4 0, the equation F (x) = 0 has a unique 
solution x* in ie where 


_ JU Go, t*) N Qo, ifh = 4 (1 — 03)" 


v= 
U (xo, 1°) N Qo, ifh < 5 (1 — a3) 


and, if a = 0, the solution x* is unique in U (xo, x); where t* and t** are given 


in Theorem 1.1. 


Proof Use in the proof of Theorem 1.1 instead of estimate (1.2.5) the analogous 
estimate 


IF xa) Il = WF Ore) — F rt) — A rei, Xk-1) OR = XK-D)IL = 
[LF Ge) — F Ge) — A Ges e-1) re — Xe-1)] + 
(A OK, Xk-1) — A K-21, Xk-1)) Xe = Xe 
<1 (FS tha, — seal + ae lxi1 — oll + 79) te — xe + 


A (xx, Xe—-1) — A (Xp—1, Xe—1) I MK — xe-1 1] < 
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Qa) 2 
l Z Mi — th—-1)" + 2 (te = te-1) te-1 + 03 (te — te-1)) 
where we used again that ||x, — xx_1|| < te — te-1, |lxx-1 — Xol] < t—1 and the con- 
dition (hi) . | 


Remark 1.6 Comments similar to Remark 1.2 (1)—(3) can follow but for method 
(1.1.11) and method (1.1.12) instead of method (1.1.5) and method (1.1.6), respec- 
tively. 
Similarly, for method (1.1.13) and method (1.1.14), we use the conditions (H”) : 
hi) = (Ay) 
hs) = (h2) 
(hy There exist real numbers a, a2, 3 satisfying 


O0<a. <a, and0< 73 
such that for each x, y € Q, 
lF@)-FY)-A@)&@W@e-yils 
Qa) 
1 (S lle — yll + a2 lly — xoll +8) x — yl. 
(hi) For each x, y, z € Q; and some 73 > 0 
|A (z, y) — A(y, x) I] < 153. 
Set a3 = 73 + 43 and further suppose 0 < a3 < 1. 
(hf) There exist x_1 € Q, x9 € Q and 7 > 0 such that A (F) Gaxay < 


L£ (B>, B,) and 
|A (F) (xo, x-1)7! F (x0) || < 2. 


(hg) = (he) 
and 
(hy) = (h7). 


Then, we can present the following semi-local convergence of method (1.1.14) 
using the conditions (H ") and the preceding notation. 


Proposition 1.7 Suppose that the conditions (H") are satisfied. Then, sequence {xp} 
generated by method (1.1.14) starting at x9 € Q is well defined in U (xo, t*), remains 
in U (xo, t*) for each n = 0, 1,2, ... and converges to a solution x* € U (xo, t*) 
of equation F (x) = 0. Moreover, if a; 4 0, the equation F (x) = 0 has a unique 
solution x* in U, where 


~ — [ U o, 0) N Qo, ifh = 5 A — a3)" 
| U Got) 1M, ifh < 4-03) 
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and, if a, = 0, the solution x™ is unique in U (xo, 5), where t* and t** are given 


in Theorem 1.1. 


Proof As in Proposition 1.5, use in the proof of Theorem 1.1 instead of estimate 
(1.2.5) the analogous estimate 
IF Goll = 


IF (xg) — F (xe-1) — A (Xx, Xe—1) (Xe = XK-1) 
+ (A (XK, Xk-1) — A K-15 Xe-2)) (Kk — XK! S 
IF (xg) — F (xe-1) — A (xe, Xe—1) (OR — Xe—1) I + 


A (xe, Xe—1) — A (Xp—1, Xe—2) I] Nove — Xe—1Ih 


Q) 
<1 (Flt — reall + a2 [lax — xoll +78) [lau — ial + 153 Ute — e- 


a 
<l (> (te — th-1)? + 02 (te — th-1) th-1 + 03 (te — n-1)) ‘ 


where we used again that ||, — x¢-1|| < te — te-1, l|Xx-1 — Xoll < te-1 and the con- 
dition (h{) . a 


Remark 1.8 Comments similar to Remark 1.2 (1)—(3) can follow but for method 
(1.1.13) and method (1.1.14) instead of method (1.1.5) and method (1.1.6), respec- 
tively. 


1.4 Applications to X-valued Fractional Calculus 


Here we deal with Banach space (X, ||-||) valued functions f of real domain [0, a], 
a > 0. All integrals here are of Bochner-type, see [8, 13]. The derivatives of f are 
defined similarly to numerical ones, see [16], pp. 83-86 and p. 93. 

In this section we apply our Newton like numerical methods to X-valued fractional 
calculus. 

We want to solve 


f(x) =0. (1.4.1) 


(1) Let 1 < v < 2,1.e. [v] = 2 ([-] ceiling of number); x, y € [0, a], a > 0, and 
f €C? ((0,a], X). 
We define the following left X-valued Caputo fractional derivatives (see [4]), 


Vv rs 1 : l-v gi 
(DY, f) @) =rary/ («—1)'” Ff" (t) dt, (1.4.2) 
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when x > y, and 


1 y 
(Dif) Q) = FOLD / (y—a'” fF" @ dt, (1.4.3) 


when y > x, where I is the gamma function. 
We define also the X-valued fractional linear operator 


f' () + (DY, f) @) SO a> y, 
(AoA) = 7 FO) + (DEF) OD Feu Ye (14.4) 
0, x=y. 


By X-valued left fractional Caputo Taylor’s formula (see [4]) we get that 


1 x 
f@-fM=fOG-y+ — | (x—t)"' Dy, f dt, forx > y, 
, (1.4.5) 


and 


1 pf 
fO-f@M=f @Mo-x*)+ — | (y—1t)’ | DY f @ dt, forx < y, 
‘ (1.4.6) 
equivalently, it holds 


1 Bf 
F@-FO=/M-yy--z ‘i Gn De pwede ey. 
wo I, 


(1.4.7) 
We would like to prove that 
Gay 
lf @&) — f(y) — (Ao (f)) @, y) @ — y)IL S 7 (1.4.8) 
for any x,y € [(0,a],O<c< 1. 
When x = y the last condition (1.4.8) is trivial. 
We assume x # y. We distinguish the cases: 
(1) x > y : We observe that 
lf ) — f &) — (Ao (f)) @ y) @ — y= (1.4.9) 


[Pore N+oo fe t)*" (Di, f) Qdt= 


! y (x = yy”! _ 
(10+ ON =P) a-»| = 
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v-1 Vv 
lol (x — 0" (Dif) dt — (Dy Os | - (1.4.10) 
rol (x —1)""! (DY, f )@dt- ar andy (De, f) atl = 
(1.4.11) 


(by [1], p. 426, Theorem 11.43) 


[eo (BNO (1) 9) at] = 


T (v) 
(by [8]) 
rw al @ =H (DLA) O- (BGA) Ola = ©, (1.4.12) 
(assume that 
| (DY, F) © — (DY, f) GI) < Arie— xP”, (1.4.13) 
for any t,x, y € [(0,a]:x >t > y, where A; < T'(v), Le. p) := _ <1), 
Therefore 
(gs) < rol C= Gon de (1.4.14) 
M @-yP @—yy? 
=sol (Ndr = PESO S15) 
We have proved that 
(x—y)? 
lf) -— Ff) -— 40) Gy) @—- ys ae (1.4.16) 
where 0 < p; < l,andx > y. 
(2) x < y : We observe that 
lf @) — Ff ) — (Ao (Ff) @, y) &@ - YI = (4.17) 


1 v—-1 Vv 
y) =f (y — t)’" Dy f (t) dt— 


j ‘ (y— x)! z 
(r @) + (DLS) Fe yp Jo »| = 
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1 y v-l Vy i (y = x)’ z 
-; wy | (y — 1)’ " DY. f (t) dt + (D%,.F) O) roaul= (1.4.18) 
| 1 f (y —t)”"' DY. f (dt — (DY. Ff) &) ar (1.4.19) 
DP (Y) Sx " me ret a 
1 y v—-l D” F 1 » sity i 7 
rv) | (y —f) wx Jd (t) -oo | (y —1f) ( on) (y) til = 
| 7 v-1 V Vv 
rv) | (y O°" (Dif O — Dix f ()) at} < (1.4.20) 


1 4 vy V V 
aon (y — 1)" [Def — DY. f | at 


(by assumption, 
|| Dy. f © — Dy, Ff (|| < Aalt- yl”, (1.4.21) 


foranyt,y,x €[0O,a]:y>t>x). 


1 
rw) 


= 


y 
(y <1)" dg |t — yl ade 


Aa [o-otto-par (1.4.22) 
(Vv) Jx 


r 


df? — ds @-y)? 
ay O~ D4 = Fay Zz 


r 


Assuming also p2 := < 1 (ie. Az < I (v)), we have proved that 


AD 
rv) 


(x — y) 
2 


Il fe) — f O) — (Ao (Ff) @, y) @& — y)II S 2 » forx<y. (1.4.23) 


Conclusion: choosing  := max (A), Az) and p := < 1, we have proved that 


AL 
T(v) 


2 
IFO -FO)- A @ N= s pL, forany x,y € [0,al. 
(1.4.24) 
This is a condition needed to solve numerically f (x) = 0. 
I) Letn-—1 <v<n,néN-—- {lI}, Le. [v] =n; x,y € [0,a],a>0, and f € 
C” ({0, a], X). 
We define the following X-valued right Caputo fractional derivatives (see [3]), 
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Dy_f QY)= ah (z—y)"” | f™ @dz, fory<x, (1.4.25) 


and 


Dy_f @)== = a ff (@—x)""! f(g) dz, forx<y. (1.4.26) 


By X-valued right Caputo fractional Taylor’s formula (see [3]) we have 


n—1 


fG) 1 fp? eddies 
fO-F=2D riage yt | (@— x)" (Dy_f) @ dz, 


(1.4.27) 
when x < y, and 


n—1 


ro-fo=y 


k=1 


YO, 


' 1 ° y— vy 
yay + 7 (—»)! (DY_A) @dz, 


(1.4.28) 
when x > y. 
We define also the fractional linear operator 


‘ia a re xr = (D’ DY’ of) (y) Tee, ee yy 
(Ao (f)) (®, y) = rt P80) Cy yt — (DY_f)@) SB y>s, 


k=l Twth ? 
0, x=y. 


(1.4.29) 
We would like to prove that 


|x — yl" 
lf @) — Ff &) — Ao) @, y) @ - YI S [= (1.4.30) 
for any x, y € [(0,a],O<c<l. 
When x = y the last condition (1.4.30) is trivial. 


We assume x # y. We distinguish the cases: 
(1) x > y : We observe that 


If @) — Ff O) — (Ao (fF) @, y) @ — y= (1.4.31) 
If () — f @)) — (Ao (f)) @ ¥) Y= XDI = 


n—-1 
f® (x) 1 % _ ; 
(= ao me Tr (v) / g-y"" (DLS) @ ic) - 
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n—-1 _ 
f®(@) = ; (x — yy’! 
(= (y — x) (DY_f) 0) eae 


k! 


i - (x —y)’! _ 
[e-» (Dy_f) (z) dz + (DY No) P= y-a| = 
(1.4.32) 


i * v-1 v 7 (x 7 yy’ _ 
rol (z—y) (Def) (z) dz — (D¥_f) OFo+h = 


lr 


yy! (DY_f) @dz- [ (— yy"! (DY_f) dz 


— yy’! ((DY_f) (@) — (DY_f) O)) dz (1.4.33) 
1 : v-1 Vy vy 
we | G-y""! (DLN) @ — (DLA) Ohaz 
(we assume that 
| (PY_f) @ — (DY_F) O)| < Arlz-yI"™”, (1.4.34) 
Xr; > 0, for all x, z, y € [0, a], with x > z > y) 
1 * = v—-l yyy = 
=D / (Z— yy @— yy dz = (1.4.35) 
= ° — j=l Ag Air @— y)" 
- worl Re oe Ta a 
(assume A, <I (v), ie. py = x <1) 
_ @—y)" 
= pj. 
n 
We have proved, when x > y, that 
cies = (1.4.36) 


If @) — f &) — (Ao (f)) @, y) & — WIS pt 


(2) y > x : We observe that 


lf @) — f(y) — Ao f)) @, y) @ - y= 
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n—-1 
f©O) 1p a 
(= a © tre | C=H) (0;.N eas) 


n—1 (k) 7 ea 
(Fe (x — yh! = (Dy_f) @) oY - y)}= 437) 
k=1 . 

7 v-l Vv i 
| = f 2 (DLA) @de- (BEA) Po 7 he or 


| : a (¢— x)" (DY_f) @dz— z (— x)" (DY_f) (#)az| = 
ry) x a= lr (v) x 7 


on a (= xy"! ((DY_f) @ — (DY_f) @)) dz (1.4.39) 
of = 9" (Dy) @ - (DEF) Col de 
(we assume that 
|(25_f) @ — (By_F) @|| s ele 2", (1.4.40) 
2 > 0, for all y, z, x € [0, a] with y > z > x) 
= a a (z— x)’ (z—x)""dz= (1.4.41) 
Tone a a — aa 
Assume now that \2 < T (v), that is py := atl <1. 
We have proved, for y > x, that 
If®—fO)- ANG @- I< m2 aL (1.4.42) 
Set A := max (\;, A2), and 
O<p:= 2 (1.4.43) 
Pr (v) 


Conclusion: We have proved that 
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If @) -— f &) — (Ao (f)) @, y) & - y)IS pao for any x, y € [0, a]. 
(1.4.44) 
In the special case of 1 < v < 2, we obtain that 
(x—y)? 
If &) — f () — (Ao (f)) @, y) @ — y)IS > (1.4.45) 


for any x,y € [0,a],O<p<l. 


R 


1 
2 


3. 


4. 


This is a condition needed to solve numerically f (x) = 0. 
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Chapter 2 
Convergence of Iterative Methods 
in Abstract Fractional Calculus 


We present a semi-local convergence analysis for a class of iterative methods under 
generalized conditions. Some applications are suggested including Banach space 
valued functions of fractional calculus, where all integrals are of Bochner-type. It 
follows [6]. 


2.1 Introduction 


Sections 2.1—2.3 are prerequisites for Sect. 2.4. 

Let B,, By stand for Banach spaces and let Q stand for an open subset of By. 
Let also U (z, p) := {u € By: ||u —z|| < p} and let U (z, p) stand for the closure of 
U (Zz, p). 

Many problems in Computational Sciences, Engineering, Mathematical Chem- 
istry, Mathematical Physics, Mathematical Economics and other disciplines can be 
brought in a form like 

F(x) =0 (2.1.1) 


using Mathematical Modeling [1-18], where F : 82 — Bp is a continuous operator. 
The solution x* of Eq. (2.1.1) is sought in closed form, but this is attainable only 
in special cases. That explains why most solution methods for such equations are 
usually iterative. There is a plethora of iterative methods for solving Eq. (2.1.1). We 
can divide these methods in two categories. 

Explicit Methods [7, 8, 12, 16, 17]: Newton’s method 


Kat = Xn — F’ (tn) * F Cn). (2.1.2) 

Secant method: : 
Xn4+1 = Xn — [%n—15 Xn3 Fl F (Xn) ; (2.1.3) 
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where [-, -; F'] denotes a divided difference of order one on Q x  [8, 16, 17]. 
Newton-like method: 
hi ae EF (2.1.4) 


where E,, = E (F) (x,) and E : Q — L (Bj, Bp) the space of bounded linear oper- 
ators from B, into Bz. Other explicit methods can be found in [8, 12, 16, 17] and 
the references there in. 

Implicit Methods [7, 10, 12, 17]: 


F (Xn) + An Qnt1 — Xn) = 0 (2.1.5) 
Xntl = Xn — Ay! F (Xn) ’ (2.1.6) 

where A, = A (Xn41, Xn) = A(F) (%n41, Xn) and A: Q x Q—> L(B,, Bo). 
There is a plethora on local as well as semi-local convergence results for explicit 


methods [1—9, 11-17]. However, the research on the convergence of implicit methods 
has received little attention. Authors, usually consider the fixed point problem 


P(x) =x, (2.1.7) 
where 
P, (x) =x + F(z) + A(x, z) (* — 2) (2.1.8) 
or 
P, (x) =z— A(x, 2)! F(z) (2.1.9) 


for methods (2.1.5) and (2.1.6), respectivelly, where z € Q is given. If P is a con- 
traction operator mapping a closed set into itself, then according to the contraction 
mapping principle [12, 16, 17], P, has a fixed point x* which can be found using the 
method of succesive substitutions or Picard’s method [17] defined for each fixed n 
by 
Vk+1,n = Py (Yen) » YOn = Xn, Xn4+1 = lim Ykn- (2.1.10) 
k—>+00 


Let us also consider the analogous explicit methods 


F (Xn) ae A ‘Can Xn) (Xn41 = Xn) — 0 (2.1.11) 
et =e A OG a) PD (2.1.12) 
F (Xn) + A (Xn, Xn—1) (Xn+1 — Xn) =0 (2.1.13) 


and 
Xnt1 =Xn — A (ns %n-1) | F (in). (2.1.14) 
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In this chapter in Sect.2.2, we present the semi-local convergence of method 
(2.1.5) and method (2.1.6). Section2.3 contains the semi-local convergence of 
method (2.1.11), method (2.1.12), method (2.1.13) and method (2.1.14). Some appli- 
cations to Abstract Fractional Calculus are suggested in Sect. 2.4 on a certain Banach 
space valued functions, where all the integrals are of Bochner-type [9]. 


2.2 Semi-local Convergence for Implicit Methods 


We present the semi-local convergence analysis of method (2.1.6) using conditions 
(S): 
(s}) F:8&C B, > By is continuous and A(x, y) € £(B), By) for each 
(x,y)EQxX Q, 
(sz) There exist @ > 0 and Qo C B, such that A (x, yy} € £L (Bo, B,) for each 
(x, y) € Qo x Qo and 
JAG») |=". 


Set Qy = 010) Qo. 


(s3) There exists a continuous and nondecreasing function w : [0, +00)? > 
[0, +00) such that for each x, y € Q 


lF @)-F(y)-AG, y)@-y)lls 
By (lx — yll, lle — xoll lly — xoll) llx — yl. 
(s4) For each x € Qo there exists y € Qo such that 
y=x-A(y, xy FP). 
(s5) For x9 € Qo and x; € Qo satisfying (s4) there exists 7 > 0 such that 
| A Ga, x0)? F @o) |] <7. 
(s6) Define g (t) := wv (n, t, t) for each t € [0, +00). Equation 
t(l—q(@t))-7n=0 


has positive solutions. Denote by s the smallest such solution. 
(s7) U (xo, 8) C , where 


sS= 4 and qo = W (7), 5S, 8). 
l= qo 
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Next, we present the semi-local convergence analysis for method (2.1.6) using 
the conditions (S) and the preceding notation. 


Theorem 2.1 Assume that the conditions (S) hold. Then, sequence {x,} generated by 
method (2.1.6) starting at x9 € Q is well defined in U (xo, s), remains in U (xo, 5) for 
eachn = 0, 1, 2, ... and converges to a solution x* € U (x0, 8) of equation F (x) = 
0. Moreover, suppose that there exists a continuous and nondecreasing function 
w, : [0, +00)* > [0, +00) such that for each x, y,z € Qy 


|F (x) — F(y) -Az, y) (x — y)|l S 
Bur x — yll, lx — xoll. lly — xoll, Ilz — xoll) Ilx — ylh 


and q, = # (, §, 5,8) <1. _ 
Then, x* is the unique solution of equation F (x) =O in U (xo, S). 


Proof By the definition of s and (s5), we have x; € U (xo, 5). The proof is based on 
mathematical induction on k. Suppose that |x, — x%—1|| < ae 4 and ||x, — xol| <s. 
We get by (2.1.6), (s2) — (55) in turn that 


Ilxeea — Xell = Ag! F Gx) |] = Ag! (F Ge) — F Gri) — Agi Ore — xe-1))| 
< Ag" || WF Ge) — F Ge-1) — Ani Oe — DI < 
A" Bap (xe — Xe—-il 5 ee—1 — xoll. Ive — Xoll) lee — xe] S 
W (n, 8,8) xe — Xe-111 = o lle — xe-1ll < 96 x1 — x0ll Sahn (2.2.1) 


and 
x41 — Xoll S llxee1 — Xe |] +. + llr — oll 


1— gk"! n 
n< = 5 
1 — qo 


<qint..tn= ; = 
— 4 


The induction is completed. Moreover, we have by (2.2.1) that for m = 0, 1, 2, ... 


=O 2 
0"). 


[xem — axl] < - 

1 — qo 

It follows from the preceding inequation that sequence {x;,} is complete in a Banach 

space B, and as such it converges to some x* € U (xo, s) (since U (xo, 8) is a closed 

ball). By letting k > +00 in (2.2.1) we get F (x*) = 0. To show the uniqueness 

part, let x** € U (xo, 5) be a solution of equation F (x) = 0. By using (2.1.6) and 
the hypothesis on w, we obtain in turn that 


\|x** - xe as ||x** — x + Ay F (xe) — Ay F (x) | = 
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|Ac' TF @*) — F Gx) — Ae @* — xx) s 


» Wxx—1 — Xoll, Ilxe — xoll . 


87" Bry (|[x** — xe 


x" — xo) |x" — xl] 


40k k+1 40k 
qu ||x** — xxl] < apt" |x — x0]. 
so lim x, = x**. We have shown that lim x, = x*,sox* =x™. |_| 
k—>+00 k>+00 


Remark 2.2 (1) The equation in (s¢) is used to determine the smallness of 77. It can 
be replaced by a stronger condition as follows. Choose ju € (0, 1). Denote by so the 
smallest positive solution of equation g (t) = py. Notice that if function q is strictly 
increasing, we can set so = q~! (). Then, we can suppose instead of (s¢) : 

(56) n<C—p)50 

which is a stronger condition than (56). 

However, we wanted to leave the equation in (s¢) as uncluttered and as weak as 
possible. 

(2) Condition (sz) can become part of condition (53) by considering 

(s3)’ There exists a continuous and nondecreasing function ¢ : [0, +00)3 > 
[0, +00) such that for each x, y € Q 


|A (x, y) | LF @) — Fy) — A(x, y) @, y)]] < 
y (lx — yl], lx — xoll, ly — xoll) Ix — yl. 


Notice that 
p (uy, U2, u3) < 1) (uy, U2, U3) 


for each u; > 0, uz > O and uz > 0. Similarly, a function y; can replace wv, for yhe 
uniqueness of the solution part. These replacements are of Mysovskii-type [7, 12, 
16] and influence the weaking of the convergence criterion in (s¢), error bounds and 
the precision of s. 
(3) Suppose that there exist 6 > 0, 6; > 0 and L € £(B,, By) with L7! e€ 
L (B, B,) such that 
Jz" ] se" 


|A@,y)-Lil sf 


and 


i= hpi <1, 
Then, it follows from the Banach lemma on invertible operators [12], and 


|e] AGN -LI< 8 B= <1 
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that A (x, y)~! € £ (Bo, B,). Let B = eae Then, under these replacements, condi- 
tion (s2) is implied, therefore it can be dropped from the conditions (5S). 

(4) Clearly method (2.1.5) converges under the conditions (S), since (2.1.6) 
implies (2.1.5). 

(5) We wanted to leave condition (s4) as uncluttered as possible, since in practice 
Eqs. (2.1.6) (or (2.1.5)) may be solvable in a way avoiding the already mentioned con- 
ditions of the contraction mapping principle. However, in what follows we examine 
the solvability of method (2.1.5) under a stronger version of the contraction mapping 
principle using the conditions (V) : 

(v1) = (1). 

(v2) There exist functions wy : [0,+00)* > [0,+00), wy»: [0,+00)* > 
[0, +-oo) continuous and nondecreasing such that for each x, y, z € Q 


IT + A(x, z) -—A(y, Zl S wi Cx — yl. lx — xoll, lly — xoll, lz — xoll) 
|A (x, z) —A(y, z)l] S we (lx — yl, lle — xoll, lly — xoll, Ilz — xoll Ila — yll 
and 


w (0, 0, 0, 0) = w2 (0, 0, 0, 0) = 0. 


w (2, t,t, t) + we (2t, t,t, t) (t+ Ilxoll), z F Xo 
w (2t, t,t, 0) + we (2t, t, t, 0) ||xol], z = Xo. 


A(t,t,t,t)= 
(v3) There exists T > 0 satisfying 


h(t,t,t,t) <1 


and 
h(t,t,0,t)t+ ||F Goll <t 
(v4) U (x0, T) © D. 


Theorem 2.3 Suppose that the conditions (V) are satisfied. Then, Eq. (2.1.5) is 
uniquely solvable for each n=0,1,2,.... Moreover, if Ay! € £L(Bo, Bi), the 
Eq. (2.1.6) is also uniquely solvable for eachn = 0, 1, 2, ... 


Proof The result is an application of the contraction mapping principle. Let x, y, z € 
U (xo, T). By the definition of operator P,, (v2) and (v3), we get in turn that 


IP: @) — P; (Yl = Il + A @, 2) — AQ, 2) @ — y) — (AGI -AVY, Dall 
< ||P+ AQ, 2) -A(y,2)1l lle — yl + IAG, 2) -AG, 2M Mell 


< [w1 (lx — yll, le — xoll. lly — oll, Iz — xoll) + 
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w2 ([|x — yl}, lx — xoll, Ily — xoll, Ilz — xoll) (lz — xoll + loll] Ila — yl 
<h(7,7,7,7) lx — yl 


and 
I| P. () — xoll < ||P: @&) — P; (xo) Il + IP, (0) — xoll 


<A (lx — xoll, lx — xoll 9, Iz — xoll) lx — xoll + IF @o)ll 
<h(7,7,0,7)7 + IF (o)ll ST. 


Remark 2.4 Sections 2.2 and 2.3 have an interest independent of Sect. 2.4. It is worth 
noticing that the results especially of Theorem 2.1 can apply in Abstract Fractional 
Calculus as illustrated in Sect.2.4. By specializing function 7), we can apply the 
results of say Theorem 2.1 in the examples suggested in Sect. 2.4. In particular for 


(2.4.1), we choose w (u1, U2, U3) = ae for u; > 0, uz > 0, uz > O and c, p are 
given in Sect. 2.4. Similar choices for the other examples of Sect. 2.4. It is also worth 
noticing that estimate (2.4.2) derived in Sect.2.4 is of independent interest but not 


needed in Theorem 2.1. 


2.3 Semi-local Convergence for Explicit Methods 


A specialization of Theorem 2.1 can be utilized to study the semi-local convergence 
of the explicit methods given in the introduction of this study. In particular, for 
the study of method (2.1.12) (and consequently of method (2.1.11)), we use the 
conditions (S’) : 

(s) F :Q Cc B, — B» iscontinuous and A (x, x) € £ (By, Bo) foreach x € Q. 

(s,) There exist 3 > 0 and Qo C B, such that A (x, x)! € £L (Bo, B)) for each 
x € Qo and 

Aways e 


Set Qy =QN Qo. 
(s,) There exist continuous and nondecreasing functions ~%p : [0, +00)? > 


[0, +00), we : [0, +00)? — [0, +00) with wo (0, 0, 0) = w2 (0, 0, 0) = 0 such that 
for each x, y € Qy 


IF (x)-F(y)-Aty, y)@—y)I 


Bro (lx — yll, le — xoll, lly — xoll) Ile — yl 


28 2 Convergence of Iterative Methods in Abstract Fractional Calculus 


and 
|A (x,y) —A(Y, Il S< By2 Cx — yIl, lle — xoll, lly — xoll) - 


Set b = Yo + Yr. 
(s,) There exist x9 € Qo and 7 > 0 such that A (xo, xo)! € £ (Bo, B,) and 


|| A (xo, x0)! F (xo) || < 0. 


(85) = (so) 

(56) = (s7). 

Next, we present the following semi-local convergence analysis of method (2.1.12) 
using the (S ‘) conditions and the preceding notation. 


Proposition 2.5 Suppose that the conditions (S’) are satisfied. Then, sequence {x;} 
generated by method (2.1.12) starting at x9 € (2 is well defined in U (xo, S), remains 
inU (xo, 8) foreachn = 0, 1, 2, ...and converges to a unique solution x* € U (x0, 8) 
of equation F (x) = 0. 


Proof We follow the proof of Theorem 2.1 but use instead the analogous estimate 
|F re) | = WF Orn) — F ei) — A Xe, Xe-1) Oe — Xk-)I] S 
IF xe) — F Ore-1) — A Ore, Xe-1) Ore — Xe-1) || + 


ICA (xp, Xe-1) — A (Xe=1, Xe—-1)) ORK = Xe-DI S 


[Wo (lx — xe—all, xx—1 — xoll , lla — xoll) + 
we (ae — xa). lei — xoll . lle — xoll)] Nee — xe-1 11 = 
W (lx — Xe-r111 5 lee—1 — oll, Ulex — Xoll) [lex — xK-111 - 


The rest of the proof is identical to the one in Theorem 2.1 until the uniqueness part 
for which we have the corresponding estimate 


[2 xn | = fa ac + APF OW — Ap ("I< 


JAc' TF @*) — F Gx) — Ae @* — xx) s 


, llxx—1 — xoll, Ilex — xoll) < 


G7" Babo (|[x** — xx 


k+1 ||x** 


q||x*—x| <4 — oll. 
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Remark 2.6 Comments similar to the ones given in Sect.2.2 can follows but for 
method (2.1.13) and method (2.1.14) instead of method (2.1.5) and method (2.1.6), 
respectively. 


2.4 Applications to Abstract Fractional Calculus 


Here we deal with Banach space (X, ||-||) valued functions f of real domain [a, D]. 
All integrals here are of Bochner-type, see [15]. The derivatives of f are defined 
similarly to numerical ones, see [18], pp. 83-86 and p. 93. 

In this section we apply the earlier numerical methods to X-valued fractional 
calculus for solving f (x) = 0. 

Here we would like to establish for [a, b] C R,a < b, f € C? (a, b], X),p EN, 
that 


|x — y|?*! 
lF O=F @—AG NY — 2s a ——___—; (2.4.1) 
pt+l 
Vx, y € [a, b], where c; > 0, and 
|A@,x) -—A(y, WIl<calx—yl, (2.4.2) 


with c. > 0, Vx, y € [a, b]. 

Above A stands for a X-valued differential operator to be defined and presented 
per case in the next, it will be denoted as A, (f), A_ (f) in the X-valued fractional 
cases, and Ag (f) in the X-valued ordinary case. 

We examine the following cases: 

(I) Here see [4, 5]. 

Let x, y € [a, b] suchthatx > y,v > 0,v ¢ N, such that p = [v], [-] the integral 
part,a=v—p(0<a<l1). 

Let f € C? ([a, b], X) and define 


y — 1 ° _ 27-1 
(PF) Qi ae al G47" 7 Od p27 =o, (2.4.3) 


the X-valued left generalized Riemann-Liouville fractional integral. 

Here I" stands for the gamma function. 

Clearly here it holds (J? f) (y) = 0. We define (J) f) (x) = 0 for x < y. By [4] 
(J? f ) (x) is a continuous function in x, for a fixed y. 

We define the subspace CY, ([a, b], X) of C? ([a, b], X): 


Cy, ((a, 6], X) = {f €C? ((a,b], X): JP_f” €C' (ly, b],X)}. (2.4.4) 


So let fec : , (a, b], X), we define the X-valued generalized »—fractional 
derivative of f over [y, b] as 
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DPS Ural”) (2.4.5) 


that is 


1 d [* 
WANW=—_—oef o-o 1 d, 046 


which exists for f € CY, ([a,b], X),fora<y<x <b. 

Here we consider f € C? ([a, b], X) such that f € CY, ([a, b], X), for every 
y € [a,b], which means also that f € Cy, ([a, b], X), for every x € [a, b] (i.e. 
exchange roles of x and y), we write that as f € CY ([a,b], X). 

That is 


1 d f” 
(LE Oy= Fd-ady Z (y —t)-* f™ (t) dt (47) 


exists for f € Cy, ([a,b], X),fora<x<y<b. 

We mention the following left generalized X-valued fractional Taylor formula 
(f € CY, (la, b], X), v > 1), see [5]. 

It holds 


p-1 


(k) 1 x 
6 F0= PAP wt ef we 9 (Dyn od, 


aa k! 
(2.4.8) 
allx, y € [a, b] with x > y. 
Similarly for f ¢ CY, (la, b], X) we have 
p-1 (k) 1 y 
fo-f@m= yO y-ntt 5 — f=" (DEF) at 
= (2.4.9) 


all x, y € [a, b] with y > x. 
So here we work with f € C? ([a, b], X), such that f € CY ([a, b], X). 
We define the X-valued left linear fractional operator 


ye 1 LO ete (D oy) (x) Tea, ey, 
(Ay (f)) &, y) = Y WPe 1 — x)! + (Df) oO) 


k=l a eee 
ye? ven x=y. 


(2.4.10) 
Notice that (see [13], p. 3) 


Ay (f)) x) — (Ag (A) I= PFET? — FPP OI] 4D 
< | fe —yl. Vx. y € [a5], 


so that condition (2.4.2) is fulfilled. 
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Next we will prove condition (2.4.1). It is trivially true if x = y. So we examine 
the case of x A y. 

We distinguish the subcases: 

(1) x > y : We observe that 


If O- F@)-A4N@ NO -*DI = 


(by (2.4.8), (2.4.10) 


If @) — f ) — As FP) @, ») & — DI 


p-l 
fOr) 1 oP ae, 

2a ee a. (eH (DSF) Odt— 2.4.12) 

p-1 
f (y) : (x _ py : 

a Bee (BNO TO 

| : v-l v 1 * i 

l/l (=2) (Dy f) oda | (x — ft) (D¥ f) (x) dt 


by [1], p. 426, Theorem 11.43 


— 1)! (Dy f) ) — (DSF) @) dt] < aa 


(by [9]) : 
cof Or MONO- (BSN colar 


(we assume that 


| (D7 F) © — (DXF) @I| <A.) It x1PA”, (2.4.14) 
forall x, y,t € [a, b] withx >t > y, with A; (vy) > Oand sup A, (vy) =: A; < ~, 
yela,b] 
also itisO < p+1—v< 1) 
vl y — pet’ dt = (2.4.15) 


p M(x — y)Ptt 
my fe -9rae= Tw) @+) 
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We have proved condition (2.4.1) 


A A @—y)Pt ; 
If) -— f @)- +N@NO-MlS FH Gap” orx > y. 
(2.4.16) 
(2) x < y : We observe that 
If) — f= A) Gy») = ay BE 
PO eh) i Pp 
2 ; 2 ome = | (yn (Def) @at— (2.4.17) 
p-l 
OD 9. . we = (pipe OD" ||. 
2. (ya - (Dif) O) 5 le 
: v-l v Vv (y = xy’ 
ol (y— 1)" (Dif) © dt — (D¢F) O) ean 
1 y eee 1 y ee - 
> / (y—1t) (D‘ f) @at-—— f (y —t) (DY f) (y) dt} = 


= 


—1)’' ((D¥ f) ® — (DY F) &)) dt 


of o-o MBEA) © - (Be) Orla 


(we assume here that 


forallx, y,t € [a, b] with y > t > x, with A2 (xv) > Oand sup 2 (x) 
xe[a,b] 


(DY f) () -— (DEF) &) | <2 @)Ie- Ph”, 


= 


2 _ v—-l _ p+l—v —_ 
eo. (yA yt) dt = 


de (y —x)P# 
_—1?’dt— ant 
oO 2 Pw @+D 


T (v) 


(2.4.18) 


(2.4.19) 


=: A2 < oo) 


(2.4.20) 


2.4 Applications to Abstract Fractional Calculus 33 


We have proved that 
dy (y—x)Pt! 
- -—(A F ; 2.4.21 
lf O) — f @)— A+) @ yO — XDI S ro @tD ( ) 
for all x, y € [a, b] such that y > x. 
Call X := max (A), A2). 
Conclusion We have proved condition (2.4.1), in detail that 
PO-'S-“ee nose "a5 ce 
—f(x)- x, x)|| < Xx, a,b]. 
Be $F YY lO) @+D y: 
(2.4.22) 


(ID) Here see [3] and [5]. 

Let x, y € [a, b] such that x < y,v > 0, v ¢ N, such that p = [v], a=v-— p 
(0<a< 1). 

Let f € C? ([a, b], X) and define 


(7 f) @) = —— =f 2" f@d,a<x<y, (2.423) 


r(v) 


the X-valued right generalized Riemann-Liouville fractional integral. 
Define the subspace of functions 


Cy (fa, b], X):= {f €C? ([a,b], X): J* fF €C! (la, y], X)}. (2.4.24) 
Define the X-valued right generalized »—fractional derivative of f over [a, y] 
as 


DY f= (1? (ref. (2.4.25) 


Notice that 


1 y 
l-a ¢(p) = —a ¢(p) 
Jy PE OO) = a | (g—x)* f (z) dz, (2.4.26) 


exists for f € C’_ ({a, b], X), and 


y- 


: _ (bet d Sp ees 
(Di_f) @) = Td a) ae =f (=a) "7" (dz. (2.4.27) 
Le. 
y = (-1)?"! v ¢(p) 
(Dy_f) @) = P= e= ae (z— x)?” f (z) dz. (2.4.28) 
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Here we consider f € C? ([a,b], X) such that f € oe ({a, b], X), for every 
y € [a,b], which means also that f € CY_ ([a, b], X), for every x € [a, b] (ie. 
exchange roles of x and y), we write that as f € CY ([a,b], X). 

That is 


Eye 


(Pr NO= Em v +1) dy 


an (z—y)?-” f™ (2) dz (2.4.29) 


exists for f € CY_ ([a,b], X),fora<y<x<b. 

We mention the following X-valued right generalized fractional Taylor formula 
(f € Cy_ ([a, b], X), v > 1), see [5]. 

It holds 


p-1 


f° oO) 1 ff St tess 
= ae | (@— x)! (DY_f) @dz, 
(2.4.30) 


all x, y € [a, b] with x < y. 
Similarly for f € CY_ ([a, b], X) we have 


p-1 


« ae ee i 
fO)—fO) = er Sy -ott of @- 9 (DLA Ode 


(2.4.31) 
allx, y € [a, b] with x > y. 
So here we work with f € C? ([a, b], X), such that f € C” ([a,b], X). 
We define the X-valued right linear fractional operator 
ft G2 - at - (Df) Oo) Se >, 
ANN Ga, erie Ba yy! — (DY f) @) G2 y >, 
= DG) p= y. 
(2.4.32) 


Condition (2.4.2) is fulfilled, the same as in (2.4.11), now for A_ (f) (x, x). 
We would like to prove that 


|x _ ylPtt 
If O-fOY-A-G)G@y@-yils<e- oor (2.4.33) 
for any x, y € [a, b], where c > 0. 
When x = y the last condition (2.4.33) is trivial. We assume x # y. 
We distinguish the subcases: 
(1) x > y : We observe that 
If @) — f O) — (A- (Ff) @ ) @ - I= (2.4.34) 


IFO) — f @) — A-G) @ )  — x1 = 
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p-1 . 
fQ) 1 ft pes 
(E52 0-94 [en wn e«)- 


k=1 


— v-l 
f® (x) 23 : a9) 7 
1 x pat Vy iy (x a adie 7 
| Tv) / @— wy" (Def) @dz+ (DEA) O FET O- »| 7 
1 a _ v—-1 v _ : Gey" 7 
| rw) i: (2 — yy’! (DY_f) @) dz — (D¥_f) 0) eel 


_l 7 v—-1l Vv " v-1 y 
TW / (z—y) (DY_f) (z) dz =) (z—y) (D7) (y) dz 


1 
lv) 


= (2.4.36) 


[ e-v (LN O- (PLANO) 
: 


1 
rv) 


/ @—y)"" |(DL_Ff) @ - (DLA) OD] dz 
(we assume that 


| (Df) @ — (DY_F) || < Arlz— yl’, (2.4.37) 


A; > 0, for all x, z, y € [a, b] with x > z> y) 


r x 
oo (z—y)’ '(z— yt’ dz= (2.4.38) 
ry) Jy 
Xr x r — yyptl — yyptl 
I [ e-va= 1 (&—y) sary, y) 
rw J, rw) pti p+l 
where p; := a > 0. 


We have proved, when x > y, that 


(x — y)Ptt 
If ~)- f@)- (A-)@,. y) &- IS ree (2.4.39) 


(2) y > x : We observe that 


If) — f(y) — (A- (PD) @ 9) @& = YI = 
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p-l 
f® (y) 1 y - F 
(= ge et woul (@—x)"! (D)_f) @ ts) 


pol v-1 
f©G) tl O=2) ) 
(x — y) (D¥_ f) @) (x—y)| = (2.4.40) 
(= k! ? rw+l) 
: — yyy-l v _ v (y = x)” 
| rv) i Ga (PT) ae, On Gay 
1 ? v—1 Vv 1 v—-1 Vy 
Pw | (g—x)”* (D_f) @dz—- aol (z—x)”"! (D¥_f) (x) dz] = 
; : (2.4.41) 
1 ¥ 
Pw / (g—x)’"! ((DS_f) @ — (D¥_f) (&)) dz] < (2.4.42) 
yf &- 9 OLA) @ = (DIA) GI de 
(we assume that 
|| (DS_f) @ — (D¥_F) @]] s Agle -—x1P”, (2.4.43) 
A2 > 0, for all y, z, x € [a, b] with y > z > x) 
<TD zl (z—x)"! (g—x)?t! "dz= (2.4.44) 
ag Go ayrh 
as)P des Se 
5 f @ oP BO ray ed 
We have proved, for y > x, that 
(y _ xyPtl 
lf @)- f£O)-GA-G)@, y) @—-y)Il <= ars ae (2.4.45) 
where p2 := ae > 0. 
Set A := max (\1, A2) and p := ry > 0. 
Conclusion We have proved (2.4.1) that 
oar" 
lf @)-fO)-A_-G)G.n@-wWIs a a (2.4.46) 


for any x, y € [a, b]. 
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(ID Let again f € C? ({a,b], X), pE N, x, y € [a, D]. 
By vector X-valued Taylor’s formula we have, see [3, 4, 18] 


f@—f 52 (y) (x yk | 1 a (x ry! (Fe (t) — f ( )) ar 
a Gait d, ie sa 
(2.4.47) 
Vx,y € [a,b]. 
We define the X-valued function 


D ro k-1 
(Ao (f)) 9) := { 2 ba. a ees (2.4.48) 


Then it holds, by [13], p. 3, 
(Ao (f)) @&, x) — (Ao (f)) &, I = FP? @- FPP OD) 24-49) 


< [fl], le -—yl,Vax,ye la, b], 


so that condition (2.4.2) is fulfilled. 
Next we observe that 


lf @) — fy) — Ao f)) @, y) @ - y= 


Pp 
f® (y) iT x : : 
> qt aa, GOP (FP OF" On at 
(2.4.50) 
pP (k) 
k=1 
— G1! [f (OP (FO -f% O)) dtl =O. (2.4.51) 
Here we assume that 
(2.4.52) 


lf? O- Ff O)| <clt-—yl, Vey €labl,c>0. 


(1) Subcase of x > y : We have that (by [9]) 


©<— ail (x — 1! | () — f O)|| ar < 


@-—t : DI i (1) TG yl dts (2.4.53) 
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‘ I: (~) T (2) (x yr! ys (p _ 1)! (x ale 
(p — I)!IT (p + 2) (p—)! (p+ 1)! 
_e@—y)t 
(p+ 1)! 
Hence oa 
(x — y) 
(Os are Te Ly, (2.4.54) 


(2) Subcase of y > x. 
We have that 


1 y pe 
©=o-m | | aay sP Oya fo) at 


< (2.4.55) 


1 y 
cof er lF? = 1” olla < 
oT | (=O dt = (2.4.56) 
c Tr (2)T (p) Pate c (p-1)! pre 
@-DIT@+) 9” -“G@-Di~ry?-” 
(y —x)?*? 
= c——___.. 
(p+ 1)! 
That is 
(y—x)P*! 
(os ce eh yw, (2.4.57) 
Therefore it holds 
|x — y|?*" 
(€) < Cit all x, y € [a, b] such that x # y. (2.4.58) 
We have found that 
_— ylptl 
If @) — f &) — (Ao (f)) @, y) @ — IS Pee a c>0, (2.4.59) 
(p+ 1)! 


forallx A y. 
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When x = y inequality (2.4.59) holds trivially, so (2.4.1) it is true for any x, y € 


[a, b]. 
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Chapter 3 
Equations for Banach Space Valued 
Functions in Fractional Vector Calculi 


The aim of this chapter is to solve equations on Banach space using iterative methods 
under generalized conditions. The differentiability of the operator involved is not 
assumed and its domain is not necessarily convex. Several applications are suggested 
including Banach space valued functions of abstract fractional calculus, where all 
integrals are of Bochner-type. It follows [5]. 


3.1 Introduction 


Sections 3.1—3.3 are prerequisites for Sect. 3.4. 
Let B,, Bz denote Banach spaces and let Q stand for an open subset of B,. Let also 
U (z, p) := {u € By: ||u —z|| < p} and let U (z, p) denote the closure of U (z, p). 
Many problems in Computational Sciences, Engineering, Mathematical Chem- 
istry, Mathematical Physics, Mathematical Economics and other disciplines can be 
brought in a form like 
F(x) =0 (3.1.1) 


using Mathematical Modeling [1-17], where F : 2 — By is a continuous operator. 
The solution x* of Eq. (3.1.1) is sought in closed form, but this can be achieved only 
in special cases. That is why most solution methods for such equations are usually 
iterative. There is a plethora of iterative methods for solving Eq. (3.1.1). We can 
divide these methods in two categories. 

Explicit Methods: Newton’s method [7, 8, 12, 16, 17] 


Kat = Xn — F’ (Xn) * F Cn). (3.1.2) 

Secant method: : 
Xn4+1 = Xn — [%n—15 Xn3 Fl F (Xn) ; (3.1.3) 
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where [-, -; F'] denotes a divided difference of order one on Q x  [8, 16, 17]. 
Newton-like method: 
hii Say — EF Ga) (3.1.4) 


where E,, = E (F) (x,) and E : Q — L (Bj, Bp) the space of bounded linear oper- 
ators from B, into Bz. Other explicit methods can be found in [8, 12, 16, 17] and 
the references there in. 

Implicit Methods: [7, 10, 12, 17]: 


F (Xn) + An (Xn41 = Xn) =0 (3.1.5) 
Xnt1 = Xn — Ay! F (Xn) ’ (3.1.6) 


where A, = A (Xn41, Xn) = A(F) On41,%) and A: Q x Q— L(B), Bo). We 
also write A (F) (x, x) = A (x, x) = A (x) foreach x € Q. 

There is a plethora on local as well as semi-local convergence results for explicit 
methods [1—9, 11-17]. However, the research on the convergence of implicit methods 
has received little attention. Authors, usually consider the fixed point problem 


Pz (x) =x, (3.1.7) 
where 
P, (x) =x + F(z) +A (x, z) (x — 2) (3.1.8) 
or 
P, (x) =z—A(x,z)' F(z) (3.1.9) 


for methods (3.1.5) and (3.1.6), respectively, where z € Q is given. If P is a con- 
traction operator mapping a closed set into itself, then according to the contraction 
mapping principle [12, 16, 17], P, has a fixed point x* which can be found using the 
method of successive substitutions or Picard’s method [17] defined for each fixed n 
by 
e+ an = Pe (Ye,n) » On = Xn, Xn4+1 = lim Ykn- (3.1.10) 
k—>+00 


Let us also consider the analogous explicit methods 


F (Xn) fe A (Xn, Xn) (Xn41 = Xn) =0 (3.1.11) 
tnt =F A etn FO) (3.1.12) 
F (Xn) + A (in, Xi) (Xn41 — Xn) =0 (3.1.13) 


and 
iu Sy Ae) FG): (3.1.14) 
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In Sect.3.2 of this chapter, we present the semi-local convergence of method 
(3.1.5) and method (3.1.6). Section3.3 contains the semi-local convergence of 
method (3.1.11), method (3.1.12), method (3.1.13) and method (3.1.14). Several 
applications to Abstract Fractional Calculus are suggested in Sect.3.4 on Banach 
space valued functions, where all the integrals are of Bochner-type [8, 14]. 


3.2 Semi-local Convergence for Implicit Methods 


We present the semi-local convergence analysis of method (3.1.6) using conditions 
(S): 
(s}) F : QC B, > By iscontinuous and A (x, y) € £(B), By) foreach (x, y) € 
Q x Q. 
(s2) There exist @ > 0 and Qo C B, such that A (x, yy} € £L (Bo, B,) for each 
(x, y) € Qo x Qo and 
JAG ye". 


Set Qy = 010) Qo. 


(s3) There exists a continuous and nondecreasing function w : [0, +00)? > 
[0, +00) such that for each x, y € Q 


lF @)-F(y)-AG, y)@-y)lls 
By (lx — yll, lle — xoll lly — xoll) llx — yl. 
(s4) For each x € Qo there exists y € Qo such that 
y=x-A(y, xy FP). 
(s5) For x9 € Qo and x; € Qo satisfying (s4) there exists 7 > 0 such that 
| A Ga, x0)? F @o) |] <7. 
(s6) Define g (t) := wv (n, t, t) for each t € [0, +00). Equation 
t(l—q(@t))-7n=0 


has positive solutions. Denote by s the smallest such solution. 
(s7) U (xo, 8) C , where 


sS= 4 and qo = W (7), 5S, 8). 
l= qo 
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Next, we present the semi-local convergence analysis for method (3.1.6) using 
the conditions (S) and the preceding notation. 


Theorem 3.1 Assume that the conditions (S) hold. Then, sequence {x,} generated by 
method (3.1.6) starting at x9 € Q is well defined in U (xo, s), remains in U (xo, 5) for 
eachn = 0, 1, 2, ... and converges to a solution x* € U (x0, 8) of equation F (x) = 
0. Moreover, suppose that there exists a continuous and nondecreasing function 
w, : [0, +00)* > [0, +00) such that for each x, y,z € Qy 


|F (x) — F(y) -Az, y) (x — y)|l S 
Bur x — yll, lx — xoll. lly — xoll, Ilz — xoll) Ilx — ylh 


and q, = # (, §, 5,8) <1. _ 
Then, x* is the unique solution of equation F (x) =O in U (xo, S). 


Proof By the definition of s and (s5), we have x; € U (xo, 5). The proof is based on 
mathematical induction on k. Suppose that |x, — x%—1|| < ae 4 and ||x, — xol| <s. 
We get by (3.1.6), (s2) — (55) in turn that 


IIxeaa — Xell = Ag! F Gx) |] = Ag! (F Ge) — F Gri) — Agi Ore — xe-1))| 
< Ag" || WF Ge) — F Ge-1) — Ani Oe — DI < 
A" Bap (xe — Xe—-il 5 ee—1 — xoll. Ive — Xoll) lee — xe] S 
W (n, 8,8) xe — Xe-111 = o lle — xe-1ll < 96 x1 — x0ll Sahn (3.2.1) 


and 
x41 — Xoll S llxee1 — Xe |] +. + llr — oll 


_ k+1 
<qgin+..+n= 40 n< i =S. 
1 — qo 1— qo 


The induction is completed. Moreover, we have by (3.2.1) that for m = 0, 1, 2, ... 


~4q 
2 ghn- 


I|Xe+m — Xx] < 

1 — qo 

It follows from the preceding inequation that sequence {x;,} is complete in a Banach 

space B, and as such it converges to some x* € U (xo, s) (since U (xo, 8) is a closed 

ball). By letting k > +00 in (3.2.1) we get F (x*) = 0. To show the uniqueness 

part, let x** € U (xo, s) be a solution of equation F (x) = 0. By using (3.1.6) and 
the hypothesis on w, we obtain in turn that 


\|x** - xe as ||x** — x + Ay F (xe) — Ay F (x) | = 
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|Ac' TF @*) — F Gx) — Ae @* — xx) s 


» Wxx—1 — Xoll, Ilxe — xoll . 


87" Bry (|[x** — xe 


x" — xo) |x" — xl] 


40k k+1 40k 
qu ||x** — xe|| < apt! |x — x0]. 
so lim x, = x**. We have shown that lim x, = x*,sox* = x. |_| 
k—>+00 k—>+00 


Remark 3.2 (1) The equation in (s¢) is used to determine the smallness of 77. It can 
be replaced by a stronger condition as follows. Choose ju € (0, 1). Denote by sg the 
smallest positive solution of equation qg (t) = py. Notice that if function q is strictly 
increasing, we can set so = q~! (). Then, we can suppose instead of (s¢) : 

(56) n<C—p)50 

which is a stronger condition than (56). 

However, we wanted to leave the equation in (s¢) as uncluttered and as weak as 
possible. 

(2) Condition (sz) can become part of condition (53) by considering 

(s3)’ There exists a continuous and nondecreasing function ¢ : [0, +00)3 > 
[0, +00) such that for each x, y € Q 


|A (x, y) | LF @) — Fy) — A(x, y) @, y)]] < 
y (lx — yl], lx — xoll, ly — xoll) Ix — yl. 


Notice that 
p (uy, u2, u3) < cH) (uy, U2, U3) 


for each u; > 0, u2 > 0 and uz > 0. Similarly, a function y; can replace y, for the 
uniqueness of the solution part. These replacements are of Mysovskii-type [7, 12, 
16] and influence the weaking of the convergence criterion in (s¢), error bounds and 
the precision of s. 
(3) Suppose that there exist 6 > 0, 6; > 0 and L € £(B,, B2) with L7! € 
L (B, B,) such that 
Jz" ] se" 


|A@,y)-Lil sf 


and 


i= hpi <1, 
Then, it follows from the Banach lemma on invertible operators [12], and 


|e] AGN -LI< 8 B= <1 
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that A (x, y)~! € £ (Bo, B,). Let 8 = a Then, under these replacements, condi- 
tion (s2) is implied, therefore it can be dropped from the conditions (S). 

(4) Clearly method (3.1.5) converges under the conditions (S), since (3.1.6) 
implies (3.1.5). 

(5) We wanted to leave condition (s4) as uncluttered as possible, since in practice 
Egs. (3.1.6) or (3.1.5) may be solvable in a way avoiding the already mentioned con- 
ditions of the contraction mapping principle. However, in what follows we examine 
the solvability of method (3.1.5) under a stronger version of the contraction mapping 
principle using the conditions (V): 

(v1) = (1). 

(v2) There exist functions wy : [0,+00)* > [0,+00), wy»: [0,+00)* > 
[0, +-oo) continuous and nondecreasing such that for each x, y, z € Q 


IT + A(x, z) -—A(y, Zl S wi Cx — yl. lx — xoll, lly — xoll, lz — xoll) 
|A (x, z) —A(y, z)l] S we (lx — yl, lle — xoll, lly — xoll, Ilz — xoll Ila — yll 
and 


w (0, 0, 0, 0) = w2 (0, 0, 0, 0) = 0. 


w, (2t, ft, t,t) + we (2t, t,t, t) (t + Ilxoll), z F x0 
w, (2t, t,t, 0) + we (28, t, t, 0) ||xoll, z = xo. 


A(t tt) = 
(v3) There exists T > 0 satisfying 


h(t,t,t,t) <1 


and 
h(t,t,0,t)t+ ||F Goll <t 
(v4) U (xo, 7) € D. 


Theorem 3.3 Suppose that the conditions (V) are satisfied. Then, Eq. (3.1.5) is 
uniquely solvable for each n=0,1,2,.... Moreover, if Ay! € £(Bo, Bi), the 
Eq. (3.1.6) is also uniquely solvable for eachn = 0, 1, 2, ... 


Proof The result is an application of the contraction mapping principle. Let x, y, z € 
U (xo, T). By the definition of operator P,, (v2) and (v3), we get in turn that 


|P, @) — Pe (Wl = + AG, 2) —AG, 2) @ —y) —- AG, 2 -AQ, DD Ell 
< ||P+ AQ, 2) -A(y,2)1l lle — yl + IAG, 2) -AG, 2M Mell 


< [w1 (lx — yll, le — xoll. lly — oll, Iz — xoll) + 
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w2 ([|x — yl}, lx — xoll, lly — xoll, lz — xoll) (lz — xoll + [ol] Ila — yl 
<h(7,7,7,7) lx — yl 


and 
I| P. () — xoll < ||P: @&) — P; (xo) Il + IP, (0) — xoll 


<A (lx — xoll, lx — xoll 9, Iz — xoll) lx — xoll + IF @o)ll 
<h(7,7,0,7)7 + IF (o)ll ST. 
| 


Remark 3.4 Sections 3.2 and 3.3 have an interest independent of Sect. 3.4. Itis worth 
noticing that the results especially of Theorem 3.1 can apply in Abstract Fractional 
Calculus as illustrated in Sect.3.4. By specializing function 7), we can apply the 


results of say Theorem 3.1 in the examples suggested in Sect. 3.4. In particular for 
p-l 
Uy 


(3.4.8), we choose ~ (uw, U2, U3) = Dp for u; > 0, u2 => 0, uz => 0 and c, p are 
given in Sect.3.4. Similar choices for the other examples of Sect. 3.4. 


3.3. Semi-local Convergence for Explicit Methods 


A specialization of Theorem 3.1 can be utilized to study the semi-local convergence 
of the explicit methods given in the introduction of this study. In particular, for 
the study of method (3.1.12) (and consequently of method (3.1.11)), we use the 
conditions (S’) : 

(s) F :Q Cc B, => Bis continuous and A (x, x) € £ (By, Bz) foreach x € Q. 

(s5) There exist 3 > 0 and Qo C B, such that A (x, x)! € £ (Bo, By) for each 
x € Qo and 

JAG. <a 


Set Q] = QA Qo. 
(s,) There exist continuous and nondecreasing functions wp : [0, +00)? > 


[0, +00), w : [0, +00)? — [0, +00) with wo (0, 0, 0) = y» (0, 0, 0) = 0 such that 
for each x, y € Qy 


IF (x)-F(y)-AQ,y)@—-y)Il< 
Bwo (lx — yll, lle — xoll. lly — xolD llx — yll 


and 
|A (x,y) ~A(y, WIS By2 Cx — yIl, lle — xoll, lly — xoll)- 
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Set w = Wo + Un. 
(s,) There exist x9 € Qo and 7 > 0 such that A (xo, xo)? € £ (Bo, B,) and 


|| A Go, x0)! F (xo) || < n- 


(85) = (56) 

(56) = (57). 

Next, we present the following semi-local convergence analysis of method (3.1.12) 
using the (S ‘) conditions and the preceding notation. 


Proposition 3.5 Suppose that the conditions (S’) are satisfied. Then, sequence {x,} 
generated by method (3.1.12) starting at x9 € (2 is well defined in U (xo, 5), remains 
inU (xo, 8) foreachn = 0, 1, 2, ...and converges to a unique solution x* € U (xo, 5) 
of equation F (x) = 0. 


Proof We follow the proof of Theorem 3.1 but use instead the analogous estimate 
| F re) Il = WF Gre) — F xe-1) — A (Xe-1, Xe-1) Oe — Xe-)II S 
IF x) — F Oe—-1) — A Ok, Xe—-1) Oe — Xk-1) || + 


ICA (x, Xe-1) — A (Xe=1, Xe—-1)) OR — X-DII S 


[Wo (llxe — xe—all 5 xe—1 — Xoll s xe — xoll) + 
> (xx — Xe—all » llxe—1 — xoll , lle — xoll)] lx — xe-111 = 
W (llx% — Xe-11l, llee—1 — xoll , llx% — xoll) llx% — xe—-1 11 « 


The rest of the proof is identical to the one in Theorem 3.1 until the uniqueness part 
for which we have the corresponding estimate 


[2° — seas = fo <n + AgTF G) - Ag (0) < 
[Az"I [F (e") — F Go — Ar @ — 8) < 


8" Brbo (|[x** — xx 


5 [Ixx-1 — oll, Ilex — oll) < 


a |x" — se] <a! fe" — xo n 


Remark 3.6 Comments similar to the ones given in Sect.3.2 can follows but for 
method (3.1.13) and method (3.1.14) instead of method (3.1.5) and method (3.1.6), 
respectively. 
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3.4 Applications to X-valued Fractional and Vector Calculi 


Here we deal with Banach space (X, ||-||) valued functions f of real domain [a, b]. 
All integrals here are of Bochner-type, see [14]. The derivatives of f are defined 
similarly to numerical ones, see [17], pp. 83-86 and p. 93. 

We want to solve numerically 


f(x) =0. (3.4.1) 


(I) Application to X-valued Fractional Calculus 

Let p € N — {1} such that p — 1 < v < p,wherev ¢ N,v > 0,ie. [v] = p([-] 
ceiling of the number), a < b, f € C? ([a, b], X). 

We define the following X-valued left Caputo fractional derivatives (see [3]) 


V os 1 ? —v-1 ¢(p) 
(Dy, f) @) =rorn/ (pai ge) Gad, (3.4.2) 


when x > y, and 


1 y 
(Dy. f) O) = Fon =f (y—1)?-"! f® (8) dt, (3.4.3) 


when y > x, where I’ is the gamma function. 
We define also the X-valued linear operator 


far Pe — y+ (DEF) @) Fete > y. 
ADVE) = | EEE GE yh + (DEF) 0) Sy > 
fon oa. 


(3.4.4) 
By X-valued left fractional Caputo Taylor’s formula (see [3]), we get that 
P=! ek) ( 1 x 
feo-fo=y OM a-y + [ 9 Digg nds, tors > y, 
= k! TY) Jy 
(3.4.5) 


and 


(k) y 
fO)-f@= 3 f — Gy=2)*+ <a f (=i)? DEF Wat, fore <y. 


k=l 
(3.4.6) 
Immediately, we observe that (by [12], p. 3) 


(Ai (f)) @ *) — (Ar (Ff) I = J FP7? @) -— FP7? ©) (3.4.7) 
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< | Ff]. lk — yl. Vx y € la, bl, 


We would like to prove that 


— y|P 
If @) -— f)- Ai) @ y) @ — IS a (3.4.8) 
for any x, y € [a, b] and some constant 0 < c < 1. 
When x = y, the last condition (3.4.8) is trivial. 
We assume x # y. We distinguish the cases: 
(1) x > y : We observe that 
If) — FY) — (Air Ff) @, ¥) & — WIL = (3.4.9) 
p-l 7 
7° (y) k 1 * v—-1 Vv 
a a ad) te w / (x —1)""! DY, f dt—- 
p-l k v 
f©O) k Vv (x — y) _ 
a a FPN Otean 


raf e-o"'@ a) (t) dt — rw ca gan! (DY, Jf) (x) dt 


(by [1], p. 426, Theorem 11.43) 


(Dy, f) () — (DY, f) () dt 


(by [8]) 
< ea / (x —1)""! | (DY f) © — (DY) @| at (3.4.10) 
(assume that 
| (DBA) © — (PSF) @| <Avit- xP”, (3.4.11) 
for any t,x, y € [a,b]: x >t > y, where A; <I (v), ie. p, = ny <1) 
ye Ge OP" d= (3.4.12) 
* ts pei, oe es eg)” 
Pw / (x —1)P'dt= [oe pl . (3.4.13) 
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We have proved that 


(x — y)? 
If @) -— f O) — (Ai (f)) @, y) & — yy) S oe (3.4.14) 


where 0 < p; < l,andx > y. 
(2) x < y : We observe that 


If @) -— Ff) — (AiG) @, y») @ - WI = (3.4.15) 
If -— F@) — (A171) OO — x)= 


3 f@) ( 


ay ro — fo (y —t)’"! (DY. f) Oat— 


k!} 
k=1 
p-1 
FPR k y G= x)" | 
d nm 7 *) (PANO FEaD 


ao i (y— a"! (DY, A\iar~—— (y — 1)" (DY, f) Oy at| = 
ro) One vx 


(3.4.16) 
—1)’' ((D¥, f) @ — (DY. Ff) ()) atl < 
a is = v 
co (y—"" [(Def) © — (Dif) On| dt G.4.17) 
(we assume that 
(Dy. f) (t) — (Dit) (y)|| < do |t- y|Po” (3.4.18) 
for any t, y,x € [a,b]: y>t>x) 
g : _ 7\y-1 _ n—v = 
= TOW (yy — PY dt = 
_ 1 Fe ear 
wate a daa ea (3.4.19) 
Assuming also 
r2 
(= =e (3.4.20) 
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(i.e. Az < I (v)), we have proved that 


— Pp 
iveoa7o<AateiG=-ai< no, forx <y. (3.4.21) 


Conclusion: Choosing \ := max (A;, A2) and p := < 1, we have proved that 


aA 
Tv) 


|x — y|? 
Pp 


Ifo) -— f)—- (Ar (f)) @, y) @— IS p , for any x, y € [a,b]. 


(3.4.22) 
(II) Application from Banach space Mathematical Analysis 
In [4], we proved the following general X-valued Taylor’s formula: 


Theorem 3.7 Let p € Nand f € C? ({A, B], X), where [A, B] C Rand (X, ||-||) 
is a Banach space. Let g € C! ({A, B]), strictly increasing, such that, g~'! € C? 
(Lg (A) g (B)]). Let any a,b € [A, B]. Then 


ied 


-1) 
f(b) = Fa+y (FoF) OO) Gay g(a) + Rp (a,b), (3.4.23) 


where 


1 b 
Rp (a,b) = 7 | (9) — 918)" (Fog) Gis) g' (s)ds (3.4.24) 


1 ” p-1 -1\()) 
=o |, GO-Or Ver)” war 


Theorem 3.7 will be applied next for g (x) = e*. One can give similar applications 
for g = sin, cos, tan, etc., over suitable intervals. 


Proposition 3.8 Let f ¢ C? ([A, B], X), p € N. Then 


p-l ol (k) 7a 
PO=f@+y ew O! 


k=1 


(ch —e)" +R, (a,b), (3.4.25) 


where ; 
b 


1 pes 
Ry (a,b) = oof (e — 1)! "(f oln)™ (ny dt 


b 
7 a, (cP —e")"" (folny” (e") e'ds, Va, b ELA, BI. 


We will use the following variant. 


3.4 Applications to X-valued Fractional and Vector Calculi a3 
Theorem 3.9 Let all as in Theorem 3.7. Then 


Pp fo -1)\&) 
POces> me) a) 


k=1 


(9 (8) — g (a))* + R% (a, B), (3.4.26) 


where 


R* (a, 8) = 


! p—1 -1)(?) —1)(P) , 
(=! / (g (8) — g(s))? (fog ) (g(s)) — (fog ) (g(a))) g (s) ds 
= (3.4.27) 


1 (8) 
— (DI 0 GH =O" (Fog) - (Fog) lay) dr, 
io 


Va, €[A, B]. 
Proof Easy. | 


Remark 3.10 Calll = f og!.Thenl, I’, ..., 1 are continuous from [g (A) , g (B)] 
into f ([A, B)). 


Next we estimate as (a, 2) : We assume that 
(fog)? O-(Fog)? GO| <Kit-g@l, 6.4.28) 
Vt,g(a) € [g (A), g(B)], where K > 0. 
We distinguish the cases: 


(i) if g (9) > g (a), then 
|X, (@,8)| s 


1 ne = ~1\(?) (p) 
(p — D! [ (g (8) — 1)? (Fog ) = (fog y (g (a))| dt < 
Lies 


K g(B) 
pH ea. (9 (8) — 1)?! (@- g(@)y" dt = (3.4.29) 


kK r'@mre 
(p— 1)! P'(p +2) 


(g (8) — g(a)?" = 


(3.4.30) 


K (p- ptt ~ G(A—g(@)P" 
eerie oS OE OFT 
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We have proved that 


= p+i 
eeaizec=- ron 


(p+ 1)! 


when g (8) > g(a). 
(ii) if g (a) > g (8), then 


|, (a, 8)|| = 
i ae | -1\() —1\(P) 
ool, Csr (For) ©- (Fos) Gay) at} = 
7 Gf 
1 a p-l -1\() —1)\(P) 
aol, 9B)" |(Fo97)” © - (Fog) Ga)| at = 
7 (3.4.32) 
K wn 2-1 p—1 
owl. (g(a) —t)' (t—g(B))? dt = 
! Jo 
K TQ? 
=) ar fee (g(a) — 9 (2)? = (3.4.33) 
K (p-))! pti » Gla) — 9 (8))?*" 
Ganley Oe = Gam 
We have proved that 
, (g(a) — 9 (3)? 
|R5 A < K a a (3.4.35) 
whenever g (a) > g (2). 
Conclusion: It holds 
_ p+1 
|| R= (a, 8)|| < c 99 OI (3.4.36) 


(p+ 1)! 


Va, €[A, B]. 
Both sides of (3.4.36) equal zero when a = £3. 
We define the following X-valued linear operator: 


(As (f)) (ty) = | Uke (fos YP (g (x) — g (1, when g (x) #9), 
fe D (x), -— y, 
(3.4.37) 


for any x, y € [A, B]. 
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Easily, we see that ({12], p. 3) 
IAs (f)) 2) = (43 (FO Y= [FPP GC) - FPP OD G.-4.38) 
< | Ff], le -— yl. Vx. y € IA, BI. 
Next we observe that (case of g (x) # g (y)) 


If @) -— fO) — 43) @& 9 G@) -g ODI = 


P og7! (k) 
— d = (g (x) — 9 (y)) + RS, x) (3.4.39) 
k=1 a 
Pp g” (k) 
(pes ) 2) G(x) oy"@0) 90) = 
k=1 


3.4.36) — |g (x) — g(y)|?*1 


R* (yy, < K 3.4.40 
Vx,y €[A, Bl: g() #g(y). 
We have proved that 
Ig (x) —gQ)|?*! 
If O-fO)- 43 (9) @ y@@)-gO)I S es 
(p+D! 
(3.4.41) 
Vx,y €[A, B] 
(the case x = y is trivial). 
We apply the above theory as follows: 
(II, ) We define 
(em 5 ) 75x yyk-1 
A ysl met Sey ie 3.4.42 
(Ast (A) 9) = | AIH Gy (3.4.42) 
for any x, y € [A, B]. 
Furthermore it holds 
jer — er P+! 
_ —(A ; * — e?) < K, ———__, 3.4.43 
| f &) — f &) — (Asi (FY) @, y) (e* - e”) | TES (3.4.43) 
V x,y € [A, B], where we assumed that 
(fo In) (1) — (f o ln) (e”) || < Ki |t—e’], (3.4.44) 


Vt,e” €[e4, 8], A < B, with K; > 0. 
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(II) Next let f ¢ C? ([—-3 +e, 4 —e]), p | N,e > O small. 
Here we define that 


Lie} Ge) en (sinx — sin y)*~!, when x # 
(A32 (f)) @, y) = =1, ra yy y, 
fe = %, 
(3.4.45) 
for any x,y € [—3 +e, 5 —e]. 
We assume that 
| (f osin”) () — (f osiné') (sin y) | < Ky|t —sinyl, (3.4.46) 


Vt,siny € [sin (-§ + €) , sin (F = e)], where K> > 0. 
It holds 
: : |sinx — sin y|?*! 
If @&) — f &) — (Asa (f)) (&, y) (sin x — sin y)|| < K2 ——~—_-—_, 
(p+ 1)! 
(3.4.47) 
Vx,ye[-Zte, $e]. 
(II;) Next let f € C? (fe, 7-—€]), p€ N,e > O small. 
Here we define 


Besecl 
(focos x (cos y) (cos x — cos sy! , when x # y, 


Pp 
(Azz (f)) (*, Y) = {= Ns ee y, 


(3.4.48) 
for any x, y € [e,m7—e]. 
We assume that 
| (f co) cos7!)” (t) — (f co) cos!) (cos y) | < K3|t—cosy|, (3.4.49) 


V t, cos y € [cose, cos (7 — €)], where K3 > 0. 
It holds 


|cos x — cos y|?*! 
(p+ 1)! 


’ 


Il f @&) — f(y) — (Az3 (f)) @, y) (cosx — cos y)|| < K3 


(3.4.50) 
Vx,y ele,a7—e]. 
Finally we give: 
(IL) Let f e C? ([-$ +e, $ —e]), p € Ne > O small. 
We define 


1) Gan y) 
(fotan~ 2 (tan y) (tan x tan y)*"!, when x # y, 


(Aza (f)) (% y) = {= 1) =, 
(3.4.51) 


for any x,y € [—3 +e, 5 —e]. 
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We assume that 


| (f otan~!) (1) — (fo tan!) (tan y) | <K,|t—tanyl, (3.4.52) 


Vt,tany € [tan (-3 + é) , tan (5 = e)], where Ky > 0. 
It holds that 


|tanx — tan y|?*! 


Il f @&) — f(y) — (Asa (F)) @, y) (tanx — tan y)|| < Kg (p+)! 


] 


(3.4.53) 
Vx,ye [-3+e¢,5—-e]. 
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Chapter 4 
Iterative Methods in Abstract Fractional 
Calculus 


The goal of this chapter is to present a semi-local convergence analysis for some 
iterative methods under generalized conditions. The operator is only assumed to 
be continuous and its domain is open. Applications are suggested including Banach 
space valued functions of fractional calculus, where all integrals are of Bochner-type. 
It follows [5]. 


4.1 Introduction 


Sections 4.1-4.3 are prerequisites for Sect. 4.4. 

Let B,, By stand for Banach spaces and let Q stand for an open subset of By. 
Let also U (z, p) := {u € By: ||u —z|| < p} and let U (z, p) stand for the closure of 
U (z, p). 

Many problems in Computational Sciences, Engineering, Mathematical Chem- 
istry, Mathematical Physics, Mathematical Economics and other disciplines can be 
brought in a form like 

F (x) =0 (4.1.1) 


using Mathematical Modeling [1-17], where F : 2 — By is a continuous operator. 
The solution x* of Eq. (4.1.1) is sought in closed form. However, this is attainable 
only in special cases, which explains why most solution methods for such equations 
are usually iterative. There is a plethora of iterative methods for solving Eq. (4.1.1). 
We can divide these methods in two categories. 

Explicit Methods [7, 8, 12, 16, 17]: Newton’s method 


Kat = Xn — F! (tn) * F Gn). (4.1.2) 

Secant method: ‘ 
Xn+1 = Xn — me Xn F | F (Xn) , (4.1.3) 
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where [-, -; F'] denotes a divided difference of order one on Q x Q [8, 16, 17]. 
Newton-like method: 
Ant = Xn — Ey F (Xm), (4.1.4) 


where E,, = E (F) (x,) and E : Q — L (Bj, Bp) the space of bounded linear oper- 
ators from B, into Bz. Other explicit methods can be found in [8, 12, 16, 17] and 
the references there in. 

Implicit Methods [7, 10, 12, 17]: 


F (Xn) = An (Xn+1 _ Xn) =0 (4.1.5) 
Xntl = Xn — Ay F (Xn) ’ (4.1.6) 


where A, = A (Xn41, Xn) = A(F) 41,2) and A: Qx Q—> L(B), Bo). We 
denote A (F) (x, x) = A (x, x) = A (x) foreach x € Q. 

There is a plethora on local as well as semi-local convergence results for explicit 
methods [1—9, 11-17]. However, the research on the convergence of implicit methods 
has received little attention. Authors, usually consider the fixed point problem 


P,(x) =X, (4.1.7) 
where 
P, (x) =x + F (z)+A(x,z) @—2) (4.1.8) 
or 
P, (x) =z— A(x, 2)! F(z) (4.1.9) 


for methods (4.1.5) and (4.1.6), respectivelly, where z € Q is given. If P is a con- 
traction operator mapping a closed set into itself, then according to the contraction 
mapping principle [12, 13, 16, 17], P, has a fixed point x* which can be found using 
the method of succesive substitutions or Picard’s method [17] defined for each fixed 
n by 
e+ a = Pe (Ye,n) » YOn = Xn, Xn+1 = lim Yk.n- (4.1.10) 
k—> +00 


Let us also consider the analogous explicit methods 


F (Xn) +A (Xn, Xn) (Xn+1 = Xn) =0 (4.1.11) 
tnt =F A etn FO) (4.1.12) 
F (Xn) + A (in, Xi) (Xn41 — Xn) =0 (4.1.13) 


and 
Xnt = Xm — A Wns n=1) | F Gp). (4.1.14) 
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In this chapter in Sect.4.2, we present the semi-local convergence of method 
(4.1.5) and method (4.1.6). Section4.3 contains the semi-local convergence of 
method (4.1.11), method (4.1.12), method (4.1.13) and method (4.1.14). Some appli- 
cations to Abstract Fractional Calculus are suggested in Sect. 4.4 on a certain Banach 
space valued functions, where all the integrals are of Bochner-type [8, 14]. 


4.2 Semi-local Convergence for Implicit Methods 


We present the semi-local convergence analysis of method (4.1.6) using conditions 
(S): 
(s}) F : QC B, > By iscontinuous and A (x, y) € £(B), By) foreach (x, y) € 
Q x Q. 
(s2) There exist @ > 0 and Qo C B, such that A (x, yy} € £L (Bo, B,) for each 
(x, y) € Qo x Qo and 
JAG ye". 


Set Qy = 010) Qo. 


(s3) There exists a continuous and nondecreasing function w : [0, +00)? > 
[0, +00) such that for each x, y € Q 


lF @)-F(y)-AG, y)@-y)lls 
By (lx — yll, lle — xoll lly — xoll) llx — yl. 
(s4) For each x € Qo there exists y € Qo such that 
y=x-A(y, xy FP). 
(s5) For x9 € Qo and x; € Qo satisfying (s4) there exists 7 > 0 such that 
| A Ga, x0)? F @o) |] <7. 
(s6) Define g (t) := wv (n, t, t) for each t € [0, +00). Equation 
t(l—q(@t))-7n=0 


has positive solutions. Denote by s the smallest such solution. 
(s7) U (xo, 8) C , where 


sS= 4 and qo = W (7), 5S, 8). 
l= qo 
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Next, we present the semi-local convergence analysis for method (4.1.6) using 
the conditions (S) and the preceding notation. 


Theorem 4.1 Assume that the conditions (S) hold. Then, sequence {x,} generated by 
method (4.1.6) starting at x9 € Q is well defined in U (xo, s), remains in U (xo, 5) for 
eachn = 0, 1, 2, ... and converges to a solution x* € U (x0, 8) of equation F (x) = 
0. Moreover, suppose that there exists a continuous and nondecreasing function 
w, : [0, +00)* > [0, +00) such that for each x, y,z € Qy 


|F (x) — F(y) -Az, y) (x — y)|l S 
Bur x — yll, lx — xoll. lly — xoll, Ilz — xoll) Ilx — ylh 


and q, = # (, §, 5,8) <1. _ 
Then, x* is the unique solution of equation F (x) = 0 in U (xo, S). 


Proof By the definition of s and (s5), we have x; € U (xo, 5). The proof is based on 
mathematical induction on k. Suppose that ||, — x%—1|| < ae 4 and ||x, — xol| <s. 
We get by (4.1.6), (s2) — (55) in turn that 


Ica — Xell = Ag! F Gx) |] = Ag! (F Ge) — F Gr) — Agi Ore — xe-1))| 
< Ag" || WF Ge) — F Ge-1) — Ani Oe — DI < 
A" Bap (xe — xxi) 5 ee—1 — xoll. Ive — Xoll) ee — xe] S 
W (n, 8,8) xe — Xe-1l1 = go lle — xe-1ll < 96 x1 — x0ll Sahn (4.2.1) 


and 
x41 — Xoll S llxee1 — Xe |] +... + llr — xoll 


1— gk"! n 
n< = 5 
1 — qo 


<qint..tn= ; = 
— 4 


The induction is completed. Moreover, we have by (4.2.1) that for m = 0, 1, 2, ... 


=O <p 
0"). 


[xem — axl] < - 

1 — qo 

It follows from the preceding inequation that sequence {x;,} is complete in a Banach 

space B, and as such it converges to some x* € U (xo, 5) (since U (xg, 8) is a closed 

ball). By letting k > +00 in (4.2.1) we get F (x*) = 0. To show the uniqueness 

part, let x** € U (xo, 5) be a solution of equation F (x) = 0. By using (4.1.6) and 
the hypothesis on w, we obtain in turn that 


\|x** - xe = \|x** — x + Ay F (xe) — Ay F (x) | = 
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|Ac' TF @*) — F Gx) — Ae @* — xx) s 


87" Bry (|[x** — xe 


» Wxx—1 — Xoll, Ilxe — xoll . 


x" — xo) |x" — xl] 


40k k+1 40k 
qu ||x** — xe|| < apt! |x — x0]. 
so lim x, = x**. We have shown that lim x, = x*,sox* = x*™. |_| 
k—>+00 k—>+00 


Remark 4.2 (1) The equation in (s¢) is used to determine the smallness of 77. It can 
be replaced by a stronger condition as follows. Choose ju € (0, 1). Denote by so the 
smallest positive solution of equation qg (t) = py. Notice that if function q is strictly 
increasing, we can set so = q~! (4). Then, we can suppose instead of (s¢) : 

(56) n<C—p)50 

which is a stronger condition than (56). 

However, we wanted to leave the equation in (s¢) as uncluttered and as weak as 
possible. 

(2) Condition (sz) can become part of condition (53) by considering 

(s3)’ There exists a continuous and nondecreasing function ¢ : [0, +00)3 > 
[0, +00) such that for each x, y € Q 


|A (x, y) | LF @) — Fy) — A(x, y) @, y)]] < 
y (lx — yl], lx — xoll, ly — xoll) Ix — yl. 


Notice that 
p (uy, u2, u3) < cH) (uy, U2, U3) 


for each u; > 0, u2 > 0 and uz > 0. Similarly, a function y; can replace y, for the 
uniqueness of the solution part. These replacements are of Mysovskii-type [7, 12, 
16] and influence the weaking of the convergence criterion in (s¢), error bounds and 
the precision of s. 
(3) Suppose that there exist 6 > 0, 6; > 0 and L € £(B,, B2) with L7! € 
L (B, B,) such that 
Jz" ] se" 


|A@,y)-Lil sf 


and 


i= hpi <1, 
Then, it follows from the Banach lemma on invertible operators [12], and 


|e] AGN -LI< 8 B= <1 
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that A (x, y)~! € £ (Bo, B,). Let 8 = a Then, under these replacements, condi- 
tion (s2) is implied, therefore it can be dropped from the conditions (5S). 

(4) Clearly method (4.1.5) converges under the conditions (S), since (4.1.6) 
implies (4.1.5). 

(5) We wanted to leave condition (s4) as uncluttered as possible, since in practice 
Eqs. (4.1.6) or (4.1.5) may be solvable in a way avoiding the already mentioned con- 
ditions of the contraction mapping principle. However, in what follows we examine 
the solvability of method (4.1.5) under a stronger version of the contraction mapping 
principle using the conditions (V): 

(v1) = (81). 

(v2) There exist functions wy :[0,+00)* — [0,+00), wy»: [0,+00)* > 
[0, +-oo) continuous and nondecreasing such that for each x, y, z € Q 


IT + A(x, z) -—A(y, Zl S wi Cx — yl. lx — xoll, lly — xoll, lz — xoll) 
|A (x, z) —A(y, z)l] S we (lx — yl, lle — xoll, lly — xoll, Ilz — xoll Ila — yll 
and 


w (0, 0, 0, 0) = w2 (0, 0, 0, 0) = 0. 


w, (2t, ft, t,t) + we (2t, t,t, t) (t + Ilxoll), z F x0 
w, (2t, t,t, 0) + we (28, t, t, 0) ||xoll, z = xo. 


A(t tt) = 
(v3) There exists T > 0 satisfying 


h(t,t,t,t) <1 


and 
h(t,t,0,t)t+ ||F Goll <t 
(wa) -U (%,.7) SD. 


Theorem 4.3. Suppose that the conditions (V) are satisfied. Then, Eq. (4.1.5) is 
uniquely solvable for each n=0,1,2,.... Moreover, if Ay! € £(Bo, Bi), the 
Eq. (4.1.6) is also uniquely solvable for eachn = 0, 1, 2, ... 


Proof The result is an application of the contraction mapping principle. Let x, y, z € 
U (xo, T). By the definition of operator P,, (v2) and (v3), we get in turn that 


|P, @) — Pe (Wl = + AG, 2) —AG, 2) @ —y) —- AG, 2 -AQ, DD Ell 
< ||P+ AQ, 2) -A(y,2)1l lle — yl + IAG, 2) -AG, 2M Mell 


< [w1 (lx — yll, le — xoll. lly — oll, Iz — xoll) + 
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w2 (l|x — yll, lx — xoll, lly — xoll, lz — xoll) (lz — xoll + [ol] Ila — yl 
<h(7,7,7,7) lx — yl 


and 
I| P. () — xoll < ||P: @&) — P; (xo) Il + IP, (0) — xoll 


<A (lx — xoll, lx — xoll 9, Iz — xoll) lx — xoll + IF @o)ll 
<h(7,7,0,7)7 + IF (o)ll ST. 


Remark 4.4 Sections 4.2 and 4.3 have an interest independent of Sect. 4.4. Itis worth 
noticing that the results especially of Theorem 4.1 can apply in Abstract Fractional 
Calculus as illustrated in Sect.4.4. By specializing function 7), we can apply the 
results of say Theorem 4.1 in the examples suggested in Sect. 4.4. In particular for 


yy ba 
(4.4.21), we choose (uy, U2, U3) = BP GntDaathath for u, > 0,u2 > 0, u3 > 0 
and A, a are given in Sect.4.4. Similar choices for the other examples of Sect. 4.4. 
It is also worth noticing that estimate (4.4.2) derived in Sect.4.4 is of independent 


interest but not needed in Theorem 4.1. 


4.3 Semi-local Convergence for Explicit Methods 


A specialization of Theorem 4.1 can be utilized to study the semi-local convergence 
of the explicit methods given in the introduction of this study. In particular, for 
the study of method (4.1.12) (and consequently of method (4.1.11)), we use the 
conditions (S’): 

(s) F :Q Cc B, — B» iscontinuous and A (x, x) € £(B,, Bo) foreach x € Q. 

(s,) There exist 3 > 0 and Qo C B, such that A (x, x)! € £L (Bo, B)) for each 
x € Qo and 

Aways e 


Set Qy =QN Qo. 
(s,) There exist continuous and nondecreasing functions ~%p : [0, +00)? > 


[0, +00), we : [0, +00)? — [0, +00) with wo (0, 0, 0) = w2 (0, 0, 0) = 0 such that 
for each x, y € Qy 


IF (x)-F(y)-Aty, y)@—y)I 


Bro (lx — yll, le — xoll, lly — xoll) Ile — yl 
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and 
|A (x, y) —A(y, y)Il S By2 Cx — yIl, lle — xoll, lly — xoll) - 


Set b = Yo + Yr. 
(s,) There exist x9 € Qo and 7 > 0 such that A (xo, xo)! € £ (Bo, B,) and 


|| A (xo, x0)! F (xo) || < 7 
(55) = (56) 
(86) = (57). 


Next, we present the following semi-local convergence analysis of method (4.1.12) 
using the (S ‘) conditions and the preceding notation. 


Proposition 4.5 Suppose that the conditions (S’) are satisfied. Then, sequence {x;} 
generated by method (4.1.12) starting at x9 € Q is well defined in U (xo, 5), remains 
inU (xo, 5) foreachn = 0, 1, 2, ...and converges to a unique solution x* € U (xo, 5) 
of equation F (x) = 0. 


Proof We follow the proof of Theorem 4.1 but use instead the analogous estimate 
IF Gre) Il = WF rn) — F e-1) — A Xe, Xe-1) Oe — Xe-I] S 
IF K) — F e-1) — A Ok, Xe—-1) Oe — Xe-1) || + 
(A ee, Xe-1) — A Xe-1, Xe-1)) (Xk — XK S 
[wo (lle — xxe—a11, x1 — xoll. le — xoll) + 
Wo (xe — xxl]. xe—1 — xoll . lle — xoll)] Ile — xe—-11] = 
W (Wee — Xe-r 1, lxe—1 — oll, loc — xoll) lle% — xe-111 - 


The rest of the proof is identical to the one in Theorem 4.1 until the uniqueness part 
for which we have the corresponding estimate 


Js" — seul = [x — e+ Ay! Go — Ag! (")] 
JAF LF Ge") — F Go — Ae (e* 2) < 


B~* Bao (]x** — xt 


» [lee—1 — xoll , Ilex — xoll) < 


k+1 ||x** 


q||x* — || <4 ~ oll - a 


Remark 4.6 Comments similar to the ones given in Sect.4.2 can follows but for 
method (4.1.13) and method (4.1.14) instead of method (4.1.5) and method (4.1.6), 
respectively. 
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4.4 Applications to X-valued Fractional Calculus 


Here we deal with Banach space (X, ||-||) valued functions f of real domain [a, b]. 
All integrals are of Bochner-type, see [14]. The derivatives of f are defined similarly 
to numerical ones, see [17], pp. 83-86 and p. 93. 

Let f : [a,b] > X such that f™ € Ly ({a, b], X), the X-valued left Caputo 
fractional derivative of order a ¢ N, a> 0, m = [a] ([-] ceiling) is defined as 
follows (see [3]): 


Qa i 1 e m—a-—1 ¢(m) 
(Dif) =p | (x —1) f™ Odt, (4.4.1) 


where I’ is the gamma function, V x € [a, b]. 
We observe that 


a I 7 m—-Aa— m 
lo Nols aos f e-em Ls” ol ar 
< [ee ll. ([ (x- jet ar) = |e loc (x — ay"° 
~ T(m—a) Va T'(m—a) (m—-a) 
__ |r le. = 
> Fina (x —a) ; (4.4.2) 
We have proved that 
a | hae | CO m—-Q ce CO m—-Q 
[PO aoa” =e OO 
(4.4.3) 


Clearly then (D? f) (a) = 0. 
Let n € N we denote D’° = D? D?... Di (n-times). 
Let us assume now that 


fec' (a,b),X), Df eC' (a, b),X), k= 1,..40: 


Dee ¢ € C ([a,b],X),n€N,O<aK<1. (4.4.4) 
By [4], we have 
{@O= 3 cca ae (Di? f) (a) + (4.4.5) 
a Tl (ia + 1) 


- — 50 ll (x joe (pre Fy (t) dt, Vx €[a, bl. 
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Under our assumption and conclusion, see (4.4.4), Taylor’s formula (4.4.5) becomes 


n 


f@-f@=>, 


i=2 


=a)" 


Fuotd (Di? f) (@ + 


1 x 
TaD | (x — 1)@tDer! (pete f) @) dt, Vx € [a,b], for0 <a <1. 
n Q) Ja 


(4.4.6) 
Here we are going to operate more generally. Again we assume 0 < a < 1, and 
f : [a,b] — X, such that f’ € C ({a, b] , X). We define the following X-valued left 
Caputo fractional derivatives: 


a 3 I . —a / 
(DS f) (x) = rd— >/ (x —t) ° f' (t) dt, (4.4.7) 


for any x > y;x, y € [a, b], and 


1 y 
(Dif) O) = aol (yt) f’ dt, (4.4.8) 


for any y > x;x,y € [a, Dd]. 

Notice D) f = f’, Di f = f’ by convention. 

Clearly here (D¢ f) , (Df) are continuous functions over [a, b], see [3]. We 
also make the convention that (Df) (x) = 0, for x < y, and (De F) (y) = 0, for 
y<x. 

Here we assume that Des: pe eC? (a, bl, XK) KS 1, vane 83 DPN y, 
DO*he f EC ({a,b], X),n EN; Vx, y € [a,b]. 
Then by (4.4.6) we obtain 


n 


a= yo Oe 
i=2 
- a 5 - i ze pores (oe) (t) dt, (4.4.9) 


Vx>y;x,y €[a,b], for0 <a <1, 
and also it holds 


n 


fO-F@=>, 


i=2 


(y — xi? 


Pune (Di? f) (~) + 


rn : Ia) [ (y = erent (DEN f) (dt, (4.4.10) 
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Vy>x;x,y € [a,b], forO <a <1. 
We define the following X-valued linear operator 


(A(f) Gy) = 


jie=l 


n (x—-y i (n+1)c (x—y)@tDer! 
Di=2 Tear (Dif) iar (Dy “f ()) Taitbath?* > ¥> 


ia-1 


— oe ie . yy y(t la-1 
Dio TIES (ey) (x) + (DP a 3 (y)) to Darth? 2 7% a 


f' (x), when x = y, 


Vx,yela,b],0<a<l. 
We may assume that 


I(A (f)) @ x) — (ACP) OY I= [Lf - £' OD (4.4.12) 
<@|x—y|,Vx,y € [a,b], with ® > 0, 


see also ([12], p. 3). 
We estimate and have: 
i) caseofx >y: 


If) -FfO-AAMYG@y@-wWIl= 


| Tr (( : 1 7 eae (Dae (4.4.13) 
n (ay y ) 
(x 23 yet ba 
= (n+l)a ay 
(DS FQ) Piet pet 


(by [1], p. 426. Theorem 11.43) 


il 
~ T(n+l)a) 


[ é _ pperiie-" (ores) (t) = oy ae 8. (x)) dt 


(by [8)) 


1 . (n+l)a-1 (n+l)a _ (nt+la 
meu) oe JDerPe Fo — (DYP*F) CI] de 


(we assume here that 


| Det Pe F@ — DEFY? F || <ALIt— x1, (4.4.14) 
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Vt,x,y € [a,b]: x >t > y, where \; > 0) 


Xj 7 (n+lja-1 _ 
<peapal (x —f) (x —t)dt= 


a Gea yee 


Ai i “Oe = yet at = _ (4.4.15) 
l((an+1)a) J, P(@ + 1a) (+ lat) 
We have proved that 
M (x _ cd iiaia ll 
lf @)-— fO) -— (ACF) @, y) & — y)Il Ss Pathe (@+tbath’ 
(4.4.16) 
forany x,y € [a,b]: x >y,0<a< 1. 
(ii) case of x < y: 
If @) — f ) -— (ACS) @, ») @ — yl = 
lf) -— fF @)- (AS) @ yO - x)= 
1 y 
(n+la-1 (n+l)a 
| Famed = | (y —t) (D¢ f) (dt (4.4.17) 
(y- ayers 
= (n+l)a _— 
(DFO) T(n+hat+h 
I 2 _ 4\(nt+Da-1 (ntla = (ntla 
oy |f oat! (memes) @ - (Det) On) at} 


1 y 
F(a + la) ] (y = OPPO | DEP A) (= (DEP F) | at 
(we assume that 


| (DETP*F) @ = (DE*P?F) OD] < Ax It yl, (4.4.18) 


Vt, y,x € [a,b]: y >t > x, where 2 > 0) 


X2 (ntla-1l ey _ 
<oatas | (y-1) (y-1dt = 


r2 a (y jo dt- A2 (y- Fo aie adie 
(n+ 1)a) J, T(ntDa)(n+tDatl) 


(4.4.19) 
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We have proved that 
NS (y _ Eo ke ay 
x)— —A x, x < ‘ 
lf o=— FO) (f) (x, y) @ — y)Il Pa@+ba(a+batD 
(4.4.20) 
Vx,yela,bl:y>x,0<a<1. 
Conclusion Let \ := max (j, A2) . It holds 

r |x = yi er tert 


If) -FO-ANMGE—-IS PE De ws cave 


1) 
Vx,y € [a,b], where0O <a <1,neEN. 

One may assume that Taba << IL, 

(Above notice that (4.4.21) is trivial when x = y.) 

Now based on (4.4.12) and (4.4.21), we can apply our numerical methods pre- 
sented in this chapter, to solve f (x) = 0. 


To have (n+ 1)a+ 1 > 2, we need to take 1 > a > 1, where n € N. 
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Chapter 5 
Semi-local Convergence in Right Abstract 
Fractional Calculus 


We provide a semi-local convergence analysis for a class of iterative methods under 
generalized conditions in order to solve equations in a Banach space setting. Some 
applications are suggested including Banach space valued functions of right frac- 
tional calculus, where all integrals are of Bochner-type. It follows [5]. 


5.1 Introduction 


Sections 5.1—5.3 are prerequisites for Sect. 5.4. 

Let B,, Bz stand for Banach spaces and let Q stand for an open subset of By. 
Let also U (z, €) := {u € By : |lu — z|| < €} and let U (z, €) stand for the closure of 
U (z, &). 

Many problems in Computational Sciences, Engineering, Mathematical Chem- 
istry, Mathematical Physics, Mathematical Economics and other disciplines can be 
brought in a form like 

F(x) =0 (5.1.1) 


using Mathematical Modeling [1-18], where F : 82 — Bp is a continuous operator. 
The solution x* of Eq.(5.1.1) is needed in closed form. This is possible only in 
special cases, which explains why most solution methods for such equations are 
usually iterative. There is a plethora of iterative methods for solving Eq. (5.1.1). We 
can divide these methods in two categories. 

Explicit Methods [7, 8, 12, 16, 17]: Newton’s method 


Kat = Xn — F’ (tn)? F Gn). (5.1.2) 

Secant method: : 
Xn4+1 = Xn — [%n—15 Xn3 Fl F (Xn) ; (5.1.3) 
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where [-, -; F'] denotes a divided difference of order one on Q x  [8, 16, 17]. 
Newton-like method: 
hii Say — EF Ga) (5.1.4) 


where E,, = E (F) (x,) and E : Q — L (Bj, Bp) the space of bounded linear oper- 
ators from B, into Bz. Other explicit methods can be found in [8, 12, 16, 17] and 
the references there in. 

Implicit Methods [7, 10, 12, 17]: 


F (Xn) + An (Xn41 = Xn) =0 (5.1.5) 
Xnt1 = Xn — Ay F (Xn) ’ (5.1.6) 


where A, = A (%n41,%n) = A (F) (Xn41, Xn) and A: Q x Q — L(B,, Bz). We let 
A(F) (x, x) = A(x, x) = A (x) for each x € Q. 

There is a plethora on local as well as semi-local convergence results for explicit 
methods [1—9, 11-17]. However, the research on the convergence of implicit methods 
has received little attention. Authors, usually consider the fixed point problem 


Pz (x) =x, (5.1.7) 
where 
P, (x) =x + F(z) +A (x, z) (x — 2) (5.1.8) 
or 
P, (x) =z—A(x,z) | F(z) (5.1.9) 


for methods (5.1.5) and (5.1.6), respectivelly, where z € Q is given. If P is a con- 
traction operator mapping a closed set into itself, then according to the contraction 
mapping principle [12, 16, 17], P, has a fixed point x* which can be found using the 
method of succesive substitutions or Picard’s method [17] defined for each fixed n 
by 
e+ an = Pe (Ye,n) » On = Xn, Xn+1 = lim Ykn- (5.1.10) 
k—++00 


Let us also consider the analogous explicit methods 


F (Xn) fe A (Xn, Xn) (Xn41 = Xn) =0 (5.1.11) 
tnt =F A etn FO) 1.12) 
F (Xn) + A (in, Xi) (Xn41 — Xn) =0 (5.1.13) 


and 
iu Sy Ae) FG): (5.1.14) 


5.1 Introduction is) 


In this chapter in Sect.5.2, we present the semi-local convergence of method 
(5.1.5) and method (5.1.6). Section5.3 contains the semi-local convergence of 
method (5.1.11), method (5.1.12), method (5.1.13) and method (5.1.14). Some appli- 
cations to Abstract Fractional Calculus are suggested in Sect. 5.4 on a certain Banach 
space valued functions, where all the integrals are of Bochner-type [9, 15]. 


5.2 Semi-local Convergence for Implicit Methods 


We present the semi-local convergence analysis of method (5.1.6) using conditions 
(S): 
(s}) F : QC B, > By iscontinuous and A (x, y) € £ (By), By) foreach (x, y) € 
Q x Q. 
(s2) There exist @ > 0 and Qo C B, such that A (x, yy} € £L (Bo, B,) for each 
(x, y) € Qo x Qo and 
JAG ye". 


Set Qy = 010) Qo. 


(s3) There exists a continuous and nondecreasing function w : [0, +00)? > 
[0, +00) such that for each x, y € Q 


lF @)-F(y)-AG, y)@-y)lls 
By (lx — yll, lle — xoll lly — xoll) llx — yl. 
(s4) For each x € Qo there exists y € Qo such that 
y=x-A(y, xy FP). 
(s5) For x9 € Qo and x; € Qo satisfying (s4) there exists 7 > 0 such that 
| A Ga, x0)? F @o) |] <7. 
(s6) Define g (t) := wv (n, t, t) for each t € [0, +00). Equation 
t(l—q(@t))-7n=0 


has positive solutions. Denote by s the smallest such solution. 
(s7) U (xo, 8) C , where 


sS= 4 and qo = W (7), 5S, 8). 
l= qo 
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Next, we present the semi-local convergence analysis for method (5.1.6) using 
the conditions (S) and the preceding notation. 


Theorem 5.1 Assume that the conditions (S) hold. Then, sequence {x,} generated by 
method (5.1.6) starting at x9 € Q is well defined in U (xo, s), remains in U (xo, S) for 
eachn = 0, 1,2, ... and converges to a solution x* € U (x, 8) of equation F (x) = 
0. Moreover, suppose that there exists a continuous and nondecreasing function 
w1 : [0, +00)* > [0, +00) such that for each x, y,z € Qy 


|F (x) -—F(y)-Atz,y)@—y)Il < 
Bu (lx — yll. lle — xoll. lly — xoll, lz — xoll Ix — yIl 


and q, = % (, §, 5,8) <1. _ 
Then, x* is the unique solution of equation F (x) = Oin U (xo, 5). 


Proof By the definition of s and (s5), we have x; € U (xo, 5). The proof is based on 
mathematical induction on k. Suppose that ||, — x%—1|| < qn and ||x, — xo|| <s. 
We get by (5.1.6), (s2) — (55) in turn that 


Ilxe41 — xel] = | Ag’ F @)|| = | Ag! F Gad — F G1) — A-1 Gre — X4-1))| 
< Ag! | WF Gx) — F Grn) — Acai Ore — xD < 
A" Brp (xe — xxi) 5 ee—1 — Xolls Ive — Xoll) lee — xe] S 
b(n, 8, 8) [Xx — Xe-1l] = go lxe — Xe-1ll S 9h lei — x0ll Sahn 2.1) 


and 
x41 — Xoll < Wxep1 — xe |] +--+ + llr — xoll 


i—ge n 
<qint-:-+n=—2-n< =e 
1— qo 1—4qo 


The induction is completed. Moreover, we have by (5.2.1) that for m = 0, 1, 2, ... 


=O 2p 
0"). 


[xem — axll < - 

1 — qo 

It follows from the preceding inequation that sequence {x;,} is complete in a Banach 

space B, and as such it converges to some x* € U (xo, 5) (since U (xg, 8) is a closed 

ball). By letting k — +00 in (5.2.1) we get F (x*) = 0. To show the uniqueness 

part, let x** € U (xo, 5) be a solution of equation F (x) = 0. By using (5.1.6) and 
the hypothesis on w, we obtain in turn that 


\|x** - xe as ||x** — x + Ay F (xe) — Ay F (x) | = 
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|Ac' TF @*) — F Gx) — Ae @* — xx) s 


» Wxx—1 — Xoll, Ilxe — xoll . 


87" Bry (|[x** — xe 


x" — xo) |x" — xl] 


40k k+1 40k 
qu ||x** — xe|| < apt! |x — x0]. 
so lim x, = x**. We have shown that lim x, = x*,sox* = x*™. |_| 
k—>+00 k—>+00 


Remark 5.2 (1) The equation in (s¢) is used to determine the smallness of 77. It can 
be replaced by a stronger condition as follows. Choose ju € (0, 1). Denote by sq the 
smallest positive solution of equation g (t) = py. Notice that if function q is strictly 
increasing, we can set sy = q~! (4). Then, we can suppose instead of (s¢) : 

(56) n<C—p)50 

which is a stronger condition than (56). 

However, we wanted to leave the equation in (s¢) as uncluttered and as weak as 
possible. 

(2) Condition (sz) can become part of condition (53) by considering 

(s3)’ There exists a continuous and nondecreasing function ¢ : [0, +00)3 > 
[0, +00) such that for each x, y € Q 


|A (x, y) | LF @) — Fy) — A(x, y) @, y)]] < 
y (lx — yl], lx — xoll, ly — xoll) Ix — yl. 


Notice that 
p (uy, u2, u3) < cH) (uy, U2, U3) 


for each u; > 0, u2 > 0 and uz > 0. Similarly, a function y; can replace y, for the 
uniqueness of the solution part. These replacements are of Mysovskii-type [7, 12, 
16] and influence the weaking of the convergence criterion in (s¢), error bounds and 
the precision of s. 
(3) Suppose that there exist 6 > 0, 6; > 0 and L € £(B,, B2) with L7! € 
L (B, B,) such that 
Jz" ] se" 


|A@,y)-Lil sf 


and 


i= hpi <1, 
Then, it follows from the Banach lemma on invertible operators [12], and 


|e] AGN -LI< 8 B= <1 
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that A (x, y)~! € £ (Bo, B,). Let B = - Then, under these replacements, condi- 
tion (s2) is implied, therefore it can be dropped from the conditions (S). 

(4) Clearly method (5.1.5) converges under the conditions (S), since (5.1.6) 
implies (5.1.5). 

(5) We wanted to leave condition (s4) as uncluttered as possible, since in practice 
equations (5.1.6) (or (5.1.5)) may be solvable in a way avoiding the already men- 
tioned conditions of the contraction mapping principle. However, in what follows we 
examine the solvability of method (5.1.5) under a stronger version of the contraction 
mapping principle using the conditions (V) : 

(v1) = (1). 

(v2) There exist functions wy : [0,+00)* > [0,+00), wy»: [0,+00)* > 
[0, +-oo) continuous and nondecreasing such that for each x, y, z € Q 


IT + A(x, z) -—A(y, Zl S wi Cx — yl. lx — xoll, lly — xoll, lz — xoll) 
|A (x, z) —A(y, z)l] S we (lx — yl, lle — xoll, lly — xoll, Ilz — xoll Ila — yll 
and 


w (0, 0, 0, 0) = w2 (0, 0, 0, 0) = 0. 


w, (2t, ft, t,t) + we (2t, t,t, t) (t + Ilxoll), z F x0 
w, (2t, t,t, 0) + we (28, t, t, 0) ||xoll, z = xo. 


A(t tt) = 
(v3) There exists T > 0 satisfying 


h(t,t,t,t) <1 


and 
h(t,t,0,t)t+ ||F Goll <t 
(v4) U (xo, 7) € D. 


Theorem 5.3 Suppose that the conditions (V) are satisfied. Then, equation (5.1.5) 
is uniquely solvable for each n = 0, 1,2, .... Moreover, if Ay! € £L(B, B), the 
Eq. (5.1.6) is also uniquely solvable for eachn = 0, 1, 2, ... 


Proof The result is an application of the contraction mapping principle. Let x, y, z € 
U (xo, T). By the definition of operator P,, (v2) and (v3), we get in turn that 


IP: @) — P; (Yl = Il + A @, 2) — AQ, 2) @ — y) — (AGI -AVY, Dall 
< ||P+ AQ, 2) -A(y,2)1l lle — yl + IAG, 2) -AG, 2M Mell 


< [w1 (lx — yll, le — xoll. lly — oll, Iz — xoll) + 
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w2 (l|x — yl}, lx — xoll, lly — xoll, Ilz — xoll) (lz — xoll + [oll] Ila — yl 
<h(7,7,7,7) lx — yl 


and 
I| P. () — xoll < ||P: @&) — P; (xo) Il + IP, (0) — xoll 


<A (x — xoll, le — xoll 9, Iz — xoll) lle — xoll + IF @odll 


<h(t,7,0,7) T+ ||F (x0) || < T. |_| 


Remark 5.4 Sections 5.2 and 5.3 have an interest independent of Sect. 5.4. Itis worth 
noticing that the results especially of Theorem 5.1 can apply in Abstract Fractional 
Calculus as illustrated in Sect.5.4. By specializing function 7), we can apply the 
results of say Theorem 5.1 in the examples suggested in Sect.5.4. In particular for 


dere 
(5.4.28), we choose (Uy, U2, u3) = a tDa@thath foru, > 0,u2 > 0, uz > 0 
and A, a are given in Sect.5.4. Similar choices for the other examples of Sect. 5.4. 
It is also worth noticing that estimate (5.4.2) derived in Sect.5.4 is of independent 


interest but not needed in Theorem 5.1. 


5.3. Semi-local Convergence for Explicit Methods 


A specialization of Theorem 5.1 can be utilized to study the semi-local convergence 
of the explicit methods given in the introduction of this study. In particular, for 
the study of method (5.1.12) (and consequently of method (5.1.11)), we use the 
conditions (S’) : 

(s}) F :Q Cc B, — Bis continuous and A (x, x) € £ (By, Bz) foreach x € Q. 

(s}) There exist G > 0 and Qo C B, such that A (x, x)! € £L (Bo, B)) for each 
x € Qo and 

JAG aso 


Set Q] = QN Qo. 
(s,) There exist continuous and nondecreasing functions %p : [0, +00)? > 


[0, +00), Wy : [0, +00)? — [0, +00) with wo (0, 0, 0) = y» (0, 0, 0) = 0 such that 
for each x, y € Qy 


IF (x)-F(y)-AQ,y)@-y)Il< 


Brbo (lx — yl. le — xoll. ly — xolD) lle — Il 


and 
|A (x,y) ~A (Cy, IL S By2 Cx — yl, lle — xoll, lly — xoll)- 


Set Y = Yo + Yr. 
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(s,) There exist x9 € Qo and 7 > 0 such that A (xo, xo)! € £ (Bo, B,) and 
|| A (xo, x0)! F (x0) || < n- 
(85) = (56) 
(56) = (87). 


Next, we present the following semi-local convergence analysis of method (5.1.12) 
using the (S ‘) conditions and the preceding notation. 


Proposition 5.5 Suppose that the conditions (S') are satisfied. Then, sequence {xy} 
generated by method (5.1.12) starting at x9 € Q is well defined in U (xo, 5), remains 
inU (xo, 5) foreachn = 0, 1, 2, ...and converges to a unique solution x* € U (xo, s) 
of equation F (x) = 0. 
Proof We follow the proof of Theorem 5.1 but use instead the analogous estimate 
IF xe) = WF a) — F tei) — A ri, Xe-1) OR = Xe-1) | S 
IF xn) — F e-1) — A Xe, Xe-1) Oe = Xe-1) | + 
I|(A (a, Xk-1) — A e-1, Xe-1)) (Oe — Xe-1) | S 
[Wo (lle — xe—all. x1 — xoll. lle — xoll) + 
We (Wx — Xe—al) 5 lee—1 — xoll s le — xoll)] lle — axa = 
W (Wee — Xe-r11), lXe—1 — oll, loc — xoll) Ile% — xe-rl - 


The rest of the proof is identical to the one in Theorem 5.1 until the uniqueness part 
for which we have the corresponding estimate 


\|x** — Xk41 | = ||x** —x+ Ap (xx) — Ae (x**) | < 
[Ac |F @*) — F Ge) — Ac (2 — 1%) S 


, llee—-1 — Xoll, lle — xoll) < 


B~* Bybo (]x** — x 


q |x" — xx] sa" Ja — x0]. = 


Remark 5.6 Comments similar to the ones given in Sect.5.2 can follows but for 
method (5.1.13) and method (5.1.14) instead of method (5.1.5) and method (5.1.6), 
respectively. 
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5.4 Applications to X-valued Right Fractional Calculus 


Here we deal with Banach space (X, ||-||) valued functions f of real domain [a, b]. 
All integrals here are of Bochner-type, see [15]. The derivatives of f are defined 
similarly to numerical ones, see [18], pp. 83-86 and p. 93. 

Let f : [a,b] > X such that f © L. ({a, b], X). The X-valued right Caputo 
fractional derivative of order a ¢ N, a > 0, m = [a] ([-] ceiling), is defined as 
follows (see [3]): 


(- 1)” 
(m — a) 


b 
(DLA W= 5 / (c= xy"! £™ (dz, (5.4.1) 
Vx € [a,b], with D™_f (x) = (—1)" f™ (x), D9_f := f, where I is the gamma 
function. 

We observe that 


1 b 
lo Nols os f e-2r LF ola 


(n) b (m) _ y\ti—a@ 
_ If ls ( 0" az) = [Flo @=2) (5.4.2) 


~ T(m — a) CT(m-a) m-a 


[PO Leo @ = x" 


T(im—a+tl) 
We have proved that 
| (25) GI s [Pel Oo < | PhO ar (5.4.3) 
Clearly here (D¢_f) (b) =0,0 <a dN. 
Let n € N. We denote 
Dl := D®_D®...D2_ (n- times). (5.4.4) 


The X-valued right Riemann-Liouville fractional integral of order a, is defined as 
follows: 


1 b 
Lj) es =a (zg —x)*"! f (2) dz, (5.4.5) 


Vx € [a, b], a := I (the identity operator). 
We denote also 
[n° := Tf Tp... (n - times). (5.4.6) 
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From now on we assume 0 < a < 1, that is m = 1. 
In [4], we proved the following X-valued right generalized fractional Taylor’s 
formula: 


Theorem 5.7 Suppose that f € C! ({a, b], X) and Dey € C! ({a, b], X), fork = 
1,..,.n EN; Dt?" f € C ({a, b], X), where 0 < a < 1. Then 


b— 


n 


= (b = ia ia 
f(x) = 2, TGa+D (D;" f) (b) + (5.4.7) 


b 
1. (= xP! (DIM F) @dz, Vx € [a,b], 


We make 


Remark 5.8 In particular, when f’ € L~ ([a,b], X), 0 <a <1, we have that 
D$_f (b) = 0, also (Dj_f) (x) = —f’ (x), and 


_] b 
(Dp ) (x) = Td-@ ‘| (z—x)° f' (z) dz, Vx € [a,b]. (5.4.8) 


Thus, from (5.4.7) we derive 


n 


f@-fo=>) 


i=2 


(b _ xo 


rGaLD (Dj° f) (b) + (5.4.9) 


1 b 
rin sal (zg —x)@tDer! (os) (<)dz, Vx € [a,b];0<a<1. 


Here we are going to operate more generally. Again we assume 0 < a < 1, and 
f : [a,b] > X, such that f’ € C ([a, b], X). We define the following X-valued 
right Caputo fractional derivatives: 


a a eal : -—a ¢/ 
(Dy_f) @) = saa (t-— x) f’ dt, (5.4.10) 


for any x < y; x, y € [a, b], and 


a _ =1 * —a ¢/ 
(PNW = / Gy" Od, G41) 


for any y < x;x,y € [a, D]. 
Notice D}_ f = —f’, Di_ f = —f’, by convention. 
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Clearly here Dy f, D¢_f are continuous functions over [a, b], see [3]. We also 
make the convention that (DS_f) (x) = 0, for x > y, and (D¢_f) (y) = 0, for y > 
x. 

Here we assume that 


DIX ft, DEF eC" (a,b, X); (5.4.12) 


k=0,1,....,2 € Nand DOTY £¢, DEM? F & C ((a, b], X); Vx, y € [a, b]; and 
O0<a<l. 
By (5.4.9) we derive the X-valued formulae: 


= (y—2)/* 
x)= = ————_——_ + 
fa@-f) DT Ga Se )O) 
1 (n+l)a—1 (n+l1)a 
pee (DSF) @az, (5.4.13) 
Vx <y; x,y € [a,b]; 0 < a < 1, and also it holds 
n (x — yyia 4 
- ————~ (D 
fO)-f@= Lrg apy (Pf) @)+ 
1 et y— n a 
(n+ la) / aes (pi f) @)dz, Ba 
y 
Vy <x;x,y€f[a,b];0<a<1. 
We define the following X-valued linear operator 
(A(f)) @, y) = 
 (y—x)8"! (ig (nt a (y—aymtder! 
4 , Tard (D\* f) = (oe f@) TCCESICE SG x<y; 


a 
Q yiam ‘ (ntla ( —y)y@tder! (5.4.15) 
5 Sains (Die f) @) — (DEF Oy) SAAT. ey, 


7 (x), when x = y, 


Vx,y €[a,b];0<a< 1. 
We may assume that 


(A (f)) x) — (ACA) OWI = Lf’ G) — f’ O) | < @lx — yl, Vx, y € Ia, 8], 
(5.4.16) 
(see [13], p. 3), with ® > 0. 
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We estimate and have: 
(i) casex < y: 


lf @) — FO) — (ACA) Gy) @ — y= 
If O- Ff@)-AAY@yY- oll = (5.4.17) 


1 y 
(n+la-1 (n+l)a = 
| r(+ha) l on} (De? 4) @az 


pee (y = al 
(0) f°) aes 


(by [1], p. 426, Theorem 11.43) 


y 
i (z—- xyarbe-l (pees (2) - pry (x)) de 
(5.4.18) 


~ T(n+)a) 


(by [9)) 


a — yy(ntbDa-1 n+l)a (n+l)a 
Sent (/ @—x) | Dy Ff (2) — Dy- f0| a) 


(we assume here that 
[oe ¢@ — DPF | sArle-al, (5.4.19) 


Vz,x,y€ [a,b]: y>z>x;, > 0) 


AI . (ntla-ly, _ 
2naspol (z—x) (z—x)dz= 


y M (y _ xyerdetl 
(n+l)a = 
harp | SE Gace WEL). 


(5.4.20) 


We have proved that 


NM (y _ xyerdetl 
lS Taha @t aD 


If @) — f£O)- ACA) Gy) & 
21) 


forany x,y € [a,b]: x <y;0<a<l. 
(ii) Case of x > y : We have 


If O-Ff@)-AAYGy&Y- xl = (5.4.22) 
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: : (ntla-l (pinta _ 
cent! @— yer (DEM 7) @de 


(x _ Daa 
C(n+ lat) 


(DF 09) 


= 


[ e-nerr (net) @ - (oF) 0») a, 
. (5.4.23) 


mae / “@— yer! |(DeP2F)  — (DE*F) On] az 


(we assume that 


1 
P((n+ Ia) | 


| (ay) (z) — (ores) (y) | <d|z- yl, (5.4.24) 


Vz,y,x € [a,b]: x => z= y; A2 > 0) 


d2 7 (ntla-1 7, = 
< aceanel (z—y) (Z-—y)dz= 


d i r = (n+l)at+l 
= / (z — yy@tD@ dz = = ey) . (5.4.25) 
T(a+ la) Jy P(t la) (a+ Dat) 
We have proved that 
ro (x = yeni 


If O-FO-AMMGy@a—-yIls Pita) (oF DetD 


26) 
forany x,y € [a,b]: x >y;0<a< 1. 
Conclusion 5.9 Let \ = max (1, Az). Then 

d |x 7 yer tert 


yils r ; 
(n+ l)a) (a+ la+]) 
(5.4.2 


If @) — f(y) — ACP) Gy) & 
7) 
Vx,y € [a,b]; whereOQ<a<I1neN. 
One may assume that 
a 


ei (5.4.28) 
T((n+ l)a) 


Above notice that (5.4.27) is trivial when x = y. 
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Now based on (5.4.16) and (5.4.27), we can apply our numerical methods pre- 


sented in this chapter to solve f (x) = 0. 


To have (n+ 1)a+ 1 > 2, we need to take 1 > a > ii Where n € N. 
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Chapter 6 
Algorithmic Convergence in Abstract 
g-Fractional Calculus 


The novelty of this chapter is the design of suitable algorithms for solving equations 
on Banach spaces. Some applications of the semi-local convergence are suggested 
including Banach space valued functions of fractional calculus, where all integrals 
are of Bochner-type. It follows [6]. 


6.1 Introduction 


Sections 6.1—6.2 are prerequisites for Sect. 6.3. 

Let B,, By stand for Banach spaces and let Q stand for an open subset of By. 
Let also U (z, p) := {u € By : ||u — zl] < p} and let U (z, p) stand for the closure of 
U (Zz, p). 

Many problems in Computational Sciences, Engineering, Mathematical Chem- 
istry, Mathematical Physics, Mathematical Economics and other disciplines can be 
brought in a form like 

F(x) =0 (6.1.1) 


using Mathematical Modeling [1-18], where F : {2 — Bp is a continuous operator. 
The solution x* of equation (6.1.1) is sought in closed form, but this is attainable 
only in special cases. That explains why most solution methods for such equations 
are usually iterative. There is a plethora of iterative methods for solving equation 
(6.1.1), more the [2, 7, 8, 10-14, 16, 17]. 

Newton’s method [7, 8, 12, 16, 17]: 


Kati = Xn — F! tq)’ F On). (6.1.2) 

Secant method: : 
Xn4+1 = Xn — [%n—15 Xn3 Fl F (Xn) ; (6.1.3) 
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where [-, -; F'] denotes a divided difference of order one on Q x  [8, 16, 17]. 
Newton-like method: 
Xnt1 = Xn — Ey F (x), (6.1.4) 


where E, = E(F)(x,) and E : Q — L(B,, Bo) the space of bounded linear 
operators from B, into By. Other methods can be found in [8, 12, 16, 17] and the 
references therein. 
In the present study we consider the new method defined for each n = 0, 1, 2, ... 
by 
Xn+1 = G (Xn) 


G (Xn41) = G (tn) — Ay F Gn), (6.1.5) 


where xo € Q is an initial point, G : Bz; — Q (B3 a Banach space), A, = 
A(F) (n41,%n) = A (Xn41,%) and A: Q x Q > L(B), Bo). Method (6.1.5) 
generates a sequence which we shall show converges to x* under some Lipschitz- 
type conditions (to be precised in Sect. 6.2). Although method (6.1.5) (and Sect. 6.2) 
is of independent interest, it is nevertheless designed especially to be used in g- 
Abstract Fractional Calculus (to be precised in Sect.6.3). As far as we know such 
iterative methods have not yet appeared in connection to solve equations in Abstract 
Fractional Calculus. 

In this chapter we present the semi-local convergence of method (6.1.5) in 
Sect.6.2. Some applications to Abstract g-Fractional Calculus are suggested in 
Sect.6.3 on a certain Banach space valued functions, where all the integrals are 
of Bochner-type [9, 15]. 


6.2 Semi-local Convergence Analysis 


We present the semi-local convergence analysis of method (6.1.5) using conditions 
(M): 
(m,) F : Q Cc B, — Bp» is continuous, G : B; — Q is continuous and 
A(x, y) € £(B), Bo) for each (x, y) € Q x Q. 
(m) There exist 3 > 0 and Qo C B; such that A (x, y)~! € L (Bo, B;) for each 
(x, y) € Qo x Qo and 
JAG yy" <8". 


Set Qy =QN Qo. 


(m3) There exists a continuous and nondecreasing function w : [0, +00)? > 
[0, +00) such that for each x, y € Q) 


IF @)-FO)-AG, y)(G@)-GO))I < 


Bw (lx — yl. lle — xoll ly — xoll) IG @) -— GOI. 
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(m4) There exists a continuous and nondecreasing function wo : [0, +00) > 
[0, +oo) such that for each x € Qy 


|G (x) — G (ro)Il_ S Yo (lx — xoll) lx — xoll - 
(ms) For x9 € Qo and x; = G (xq) € Qo there exists 7 > 0 such that 
||A G1, x0)! F (xo) || < n- 
(m6) There exists s > 0 such that 
v(m, s,s) <1, 
Yo (8) < 1 


and 


(Go) =x <s2—"— 
— 40 


where go = w (7), 5, 8). 

(m7) U (x0, 5) CQ. 

Next, we present the semi-local convergence analysis for method (6.1.5) using 
the conditions (/) and the preceding notation. 


Theorem 6.1 Assume that the conditions (M) hold. Then, sequence {x,} generated 
by method (6.1.5) starting at xy € Q2 is well defined in U (xq, 5), remains in U (xo, S) 
for each n = 0,1,2,... and converges to a solution x* € U (xo, 5) of equation 
F (x) = 0. The limit point x* is the unique solution of equation F (x) = 0 in 
U (x0, S) : 


Proof By the definition of s and (ms), we have x; € U (xo, s). The proof is based on 
mathematical induction on k. Suppose that ||x, — xx—-1|| < qn and ||x, — xol| <s. 
We get by (6.1.5), (m2) — (ms) in turn that 
IG Gs) — G oll = Ag’ F Go|] = 
| Ac! F Ge) — F Ge-v) — Aci (G On) — G G1) | 
< Ag! | WF Ge) — F Gre) — Aci (G Oe) — G OI S 
B7" Bw (xe — xxi) 5 lxe—1 — Xoll. lye — xoll) IG (xe) — G @%-VIl < 


Wn, 8,5) IG (xe) — G Xe) Il = G0 IG (xe) — G (ea) II < 96 Ile1 — xoll < 967 
(6.2.1) 
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and by (m6) 
Ilxx+1 — Xoll = IG (xe) — Xoll S IG Gx) — G Co) Il + IG (0) — oll 
S Yo (lle — xoll) Ilex — xoll + IG Go) — xoll 
< Yo (s) 5 + [IG (xo) — xoll < ». 


The induction is completed. Moreover, we have by (6.2.1) that for m = 0, 1, 2, ... 


[een — mall S "gb, 

— 90 

It follows from the preceding inequation that sequence {G (x;,)} is complete in a 

Banach space B, and as such it converges to some x* € U (xo, 5) (since U (xo, 5) 

is a closed ball). By letting k — +00 in (6.2.1) we get F (x*) = 0. We also get 

by (6.1.5) that G (x*) = x*. To show the uniqueness part, let x** € U (xo, 5) bea 

solution of equation F (x) = 0 and G (x**) = x**. By using (6.1.5), we obtain in 
turn that 


\|x** — G Gea) || = |x — G Ge) + Ap! F Gx) — Ap F (x) | < 
4c" |F (e) — F Gx) — At (6) - Gen)|| s 


B~* Bao (]x** — xe 


» Mxez1 — xoll llxe — xoll) |G (x) — G @x)| < 


qo |G (x) — G Gx) S40" x — x0 


’ 
so lim x, = x**. We have shown that lim x, = x*,sox* = x. | 
k—>-+00 k—>+00 


Remark 6.2 (1) Condition (m2) can become part of condition (m3) by considering 
(m3)’ There exists a continuous and nondecreasing function ¢ : [0, +00)? > 
[0, +00) such that for each x, y € Q 


|A (x, y) | LF @&) — FQ) — A(x, ») (E(x) — G(y))]] < 
y (Ix — yl], lla — xoll, ly — xl) IG @) — GO)IL. 


Notice that 
p (uy, U2, u3) < 1) (uy, U2, u3) 


for each u, > 0, u2 > O and uz > O. Similarly, a function y; can replace y, for the 
uniqueness of the solution part. These replacements are of Mysovskii-type [7, 12, 
16] and influence the weaking of the convergence criterion in (m6), error bounds and 
the precision of s. 


6.2 Semi-local Convergence Analysis 91 


(2) Suppose that there exist @ > 0, 3; > O and L € L(B,, Bo) with L~! € 
L£ (B>, B,) such that 
Je] se" 


|A@,y)-Lll sf 


and 


ao", = 1. 
Then, it follows from the Banach lemma on invertible operators [12], and 
|e" AG») -LI< 6" =h <1 


that A (x, y)~! € £ (Bo, B,). Let B = = Then, under these replacements, condi- 


tion (mz) is implied, therefore it can be dropped from the conditions (M). 

Remark 6.3 Section 6.2 has an interest independent of Sect. 6.3. It is worth noticing 
that the results especially of Theorem 6.1 can apply in Abstract g-Fractional Calculus 
as illustrated in Sect. 6.3. By specializing function 7, we can apply the results of say 


Theorem 6.1 in the examples suggested in Sect.6.3. In particular for (6.3.4), we 
choose for u; > 0, u2 > 0, u3 > 0 


OT al 
BT (n+ Day(n+Dat+’ 


W (4, U2, U3) = 


if |g (x) — 9 (y)| < tm for each x, y € [a,b]; 


TT Saale 
BT (m+ Day(n+Dat+l)’ 


W (U4, U2, U3) = 


if |g (x) — 9g (y)| < & |x — y|| for each x, y € [a, b] and pz = & |b — al; 


ala 
BT (n+ Day(n+Dat+l)’ 


W (U1, U2, U3) = 


if |g (x)| < & for each x, y € [a, b] and paz = 2&3. 

Other choices of function ~ are also possible. 

Notice that with these choices of function 7 and f = F and g = G, crucial 
condition (m3) is satisfied, which justifies our definition of method (6.1.5). We can 
provide similar choices for the other examples of Sect. 6.3. 
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6.3 Applications to X-valued Modified g-Fractional 
Calculus 


Here we deal with Banach space (X, ||-||) valued functions f of real domain [a, D]. 
All integrals here are of Bochner-type, see [15]. The derivatives of f are defined 
similarly to numerical ones, see [18], pp. 83-86 and p. 93. 

Let0 <a < 1,m= [a] =1({-] ceiling of number), g is strictly increasing and 
g € C! (fa,b]), 9"! € C ([g (a), g (b)]). Assume that f € C! ({a, b], X). In both 
backgrounds here we follow [5]. 

(1) The X-valued right generalized g-fractional derivative of f of order a is defined 
as follows: 


y = Z -—a) —1)/ 
(Dp gf) @) = moa! (g(t) —g@)) ° 9 O(fog') G@)dt, 


(6.3.1) 
a<x<b. 
If 0 < a < 1, then (Deaf) € C ([a, b], X), see [4]. 
Also we define 
(Dh. f) @) = -((fog"') og) @), (6.3.2) 
(D}_., f) («) := f (x), Vx € [a,b]. 
When g = id, then 
Def OSD Gi @ SLT: (6.3.3) 
the usual X-valued right Caputo fractional derivative, see [3]. 
Denote by 
Dye g = Dy .gDp_.g---Do_zg (n times), n € N. (6.3.4) 


We consider the X-valued right generalized fractional Riemann-Liouville integral 


b 
/ g@)—g@)*' OF OMdt,asx<b. (635) 


a 1 
(Up_.gf) @) = Te S, 


Also denote by 


Lae = I; :glb- gulb- 1g (n times). 62-9) 


We will be using the following X-valued modified g-right generalized Taylor’s for- 
mula 


Theorem 6.4 ([5]) Let f € C!({a,b], X), g € C! ({a,b]), strictly increasing, 
such that g~' € C' ([g (a), g (b)]). Suppose that Fy, := De at, k=1,...,n, fulfill 
F, € C! ({a, b], X), and Fis; € C ([a, b], X), where0 <a <1,n EN. Then 
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3 (g (b) — g (x))'* 


fa)-fOo= Pdat) 


(Die. f) )+ (6.3.7) 
i=l 


awa! (g(t) —g (x) re 9 (1) (pes) (t) dt, 


Vx eéla,b]. 


Here we are going to operate more generally. We consider f € C! ({a, b], X). 
We define the following X-valued right generalized g-fractional derivative: 


(Desf) = ol (g(t) —gQ))* g O (fog) g@)dt, 


(6.3.8) 
ala<x<y; ye [a,b], 


(DI_,f)@:= -((fog)'og) OO, ¥xe lal. (6.3.9) 
Similarly we define: 


D>) Oe a. (g(t}-gQ))* g (fog!) G@)adt, 


(6.3.10) 
alla <y<x;x€[a,b], 


(DL f)O):=- ((F og ')'o 4) (y), Vy € [a,b]. (6.3.11) 


When 0 < a < 1, De a f and D®.. j f are continuous functions on [a, b] , see [5]. 
Note here that by convention we have that 


Vg 


and (6.3.12) 
(eee) (y) =0, for y >x 


(ve : f) (x) =0, forx > y 


Denote by 
PSD ee De Va ela): (6.3.13) 
We assume that 
F?, Fg € C' (a, b], X), and F?,,, FX,, € C (a, b], X), (6.3.14) 


k=1,...,n,Vx,y € [a,b];0<a< 1. 
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We also observe that (0 < a < 1) ({9]) 


| (Dh if Nwolseqrs [i Gg) -g@)*g O|(Fog7) Ge] a = 
(6.3.15) 

(roa 'yes| 

ld-a) 


b 
costa P| / (g(t)-—g(x))*' gg (dt = 


| (/ i v') i Al an (g (b) — g(x))'~° _ 
r(l—-a) l-a 


loreal, 
Tl (2-a) 


1°l (g (b) — g (x))'*, Vx € [a,b]. 


We have proved that 


=q\" 
a | (f fies ) ve oer l—a 
|(De_.,f) Ol < Poa (g (b) — 9 (x)) (6.3.16) 


\ron'ea, 
= lr Q-a) 


laitl (oe) = (a))-*, Var ye la Bl: 


Clearly here we have 
(D¢_,,f) ®) =0, 0<a<1. (6.3.17) 


In particular it holds 
(DEF ) (x) = (Deaf) (vy) =0, Vx,y€ [a,b]; O0<a<l. (6.3.18) 
By (6.3.7) we derive 


i=2 


(Di*., f) O)+ (6.3.19) 


ae ra eer — (n+l)a-1 7 ogee 
rin a a a (g(t) — g(x) g' ) (DMP L) Oar, 


Vx < y; x,y € [a,b]; 0 < a < 1, and also it holds: 
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foy- f= eee 


i=2 


(Dit, ,f) (@)+ (6.3.20) 
a _ (n+lja-1 7 9 al 
Tint orl (g(t) — g(y)) g (1) (D oo (t) dt, 


Vy <x;x,y€f[a,b];0<a<1. 
We define also the following X-valued linear operator 


(Ai (f)) @, y) = 
( (x)) ic (n+l)a (g(y)—g(x)) Der! 
ps ee - (Die. f) (y) ~ (De) “ft (x)) se ee eee x<y, 


(g@)-9))""! ( pia (ntl)a (g@)-ggy) "te! 
= Tat] - (Die, iS) (x) — (oe ca (»)) “Ta@thbath ? ~ 7 9: 


7" (x), when x = y, 


(6.3.21) 
Vx,y €f[a,b];0<a< 1. 
We may assume that (see [13], p. 3) 


(Ai (f)) @, x) — (Ar (f)) Oy) = If’) — OD) 
=|(fog')@@)-(fog') GO) < Pl9@—- gO), 3.22) 
Vx,y € [a,b]; with ® > 0. 
We estimate and have 


(i) casex <y: 


If O-fO)- AiG) @ y@G@-gONI = 


1 y 
| sae! (g(t) -—g@yerPe"! 7 (1) (a7) (t) dt— 


= (nt+la 
(ery (x)) (g(y) — 9 @)) (6.3.23) 


T(m+latl) 


(by [1], p. 426, Theorem 11.43) 


1 
Tas Do). 


| 7 G1) = g PO () (DEEP L) © — (DMP F) @) at 
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(by [9)) 
1 


<F@tDa) 
| “a — gyn gf) | (DEP) © — (DEPP) CO] ae 6.3.24) 
(we assume that 


(0 pide f)o- (D' pete PA) co <M lg —g(@)I, (6.3.25) 


Vt,x,y€ [a,b]: y>t>x; 1, >0) 


ee I (g(t) — g(x)! oD) (g(t) — g(x) dt = (6.3.26) 
~T(a+la)J, * 


pa Ae = (nt+l)a 
= 1 =| (g(t) — 9g ()) g (t)dt = 


1 (9 (y) — g (x) Ore? 


3.27 
Pantha) (a+hatl iSeet 
We have proved that 
Ifo -fO)- (Ai (f) & wWg@)-gO)I S 
x1 OG) eat 064) Waaaaele 
Pim+ha (a+Da+h (62-28) 
forany x,y € [a,b]: x <y;0<a< 1. 
(li) casex > y: 
If @)- f O)- (Ai (f) @. y) 9G @) -— 9g O)I = 
Ifo - f @)- (Ai) @ GO) -—g@)Il= (6.3.29) 
aed ar _ (n+l)a-1 (n+la _ 
| ; aa ar . (g(t) — 9) a (1) (De) f) (1) dt 
(nt-Ia (g (x) — g (y))@*D* 
(Pes ‘) "T+ Dat) ere 


1 
~ T+ ta) 
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= 


| i “Gaon gf (1) (DEP F) = (DELP F) On) at 


—— aT [ (9 (8) = g MPO! g! () [DEP F () — DELP” F | at 
(we assume that 
[oer s co — DEP FO] < 29 -gOI1, (63.31) 


Vt, y,x €[a,b]:x >t > y; 2. >0) 


— : _ (nt+Da-1 _ _ 
= T+ a) / @@— 90) J) (gt) —g(y))dt = (6.3.32) 


—— = (n+l)a 7 
mca (g(t) -9gQ)) g (t)dt = 


2 (g(x) —gQ))Orret! 
T(nt+la (a+Dath) ~ 


We have proved that 
If O-FfO-A)@ y@@- gO S 


do G@j=9gQ)yrr 
Tq(a+la) (a+ljatl) 


(6.3.33) 


Vx,yela,bl:x>y0<a<1. 


Conclusion 6.5 Set \ = max (1, 2). We have proved that 


If O-FO- AG) y@G@)—- IOI S 


r Ig (x) — g (yy [Order 


T(n+la) (n+la+l) (6.3.34) 


Vx,yé[a,b];0<a<LneNn. 


(Notice that (6.3.34) is trivially true when x = y.) 
One may assume that 


A 


———— <] (6.3.35) 
(n+ 1a) 
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Now based on (6.3.22) and (6.3.34), we can apply our numerical methods presented 
in this chapter to solve f (x) = 0. 

To have (n+ 1)a+ 1 > 2, we need to take 1 > a > <7 Where n EN. 

Some examples of g follow: 


g(x) =e*, x E€ [a,b] CR, 

g(x) = sinx, 

g (x) = tanx, 

where x € [-3 +e, z _ ae € > O small. 


(6.3.36) 


(I) The X-valued left generalized g-fractional derivative of f of order a is defined 
as follows (see [5]): 


1 * , -~i\! 
Ping) =p | @@-go"d O (Fog) GOnat, 


(6.3.37) 
Vxela,b]. 
If0 <a < 1, then (02.7) € C ({a, b] , X) (see [5]). 
Also, we define 
Das G)= (Gj og ')'o 4) (x), (6.3.38) 
Dis.gf @) =f), Vx € [a,b]. 
When g = id, then 
Divigt = Diciad = DES: 
the usual X-valued left Caputo fractional derivative (see [4]). 
Denote by 
Dit ig = DatigPatig:Patg (1 times), n EN. (6.3.39) 


We consider the X-valued left generalized fractional Riemann-Liouville integral 


(see [5]) 
Cth) @) = Pia ol ga)-gM)"'g OF Odt,asx<b. (6.3.40) 


Also denote by 


na ya a a 
Lae ie ere 2 ee 


(n times). (6.3.41) 


We will be using the following X-valued modified g-left generalized Taylor’s for- 
mula. 


Theorem 6.6 ([5]) Letp0 <a <1,neEN, f € C! ({a,b], X), g € C' ({a, b)), 


strictly increasing, such that g~'! € C! ({g(a),g(b)]). Let Fy : =F ass —— 


1, ...,n, that fulfill F, € C! ({a, b], X), and F,4; € C ({a, b], X). Then 
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fa-f@=-> eee 


i=] 


(Dig.,f) @+ (6.3.42) 


Feaha | ose (Dis) de 


Vx € [a,b]. 


Here we are going to operate more generally. We consider f € C! ({a,b], X). 
We define the following X-valued left generalized g-fractional derivative: 


De ae. (g@)-g9@)* 9 (fog) GM)dt, 


(6.3.43) 
forany y<x <b; x,y e€[a,b], 


(D Hed, (x) = (fog) g (x)), Vx € [a,b]. (6.3.44) 


Similarly, we define 


(Deeg f) O) = ool (g(y)-g@)° gd (fog) G@)adt, 


(6.3.45) 
foranyx <y <b; x,y €[a,b], 


(Di..,f) 0) = (fog) GQ), Vy €la, bl. (6.3.46) 


When 0 < a < l, Dis a and De ad are continuous functions on [a, b], see [5]. 
Note here that by convention, we have that 


(2537) (x) = 0, when x < y, 


and (6.3.47) 
(D2, ro (y) =0, when y <x. 


Denote by 
Go =DE fF GiaHDy Tf. Vayeiae: (6.3.48) 


We assume that 


GG, eC’ ((a,8) :X),. andl Cg, Gg 


€ C (fa, b], X), (6.3.49) 


n+l? 


k=1,...,n,Vx,y € [a,b];0<a< 1. 
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We also observe that (0 < a < 1) (by [9]) 
[ON Os f o@-soso|(for') Goolas 


loo ed 
rd-a) 


a / GO -gOy gd dt = 


| (fog!) 1 (g@) = g@)" = 
Tr (1—a) —@ 


(6.3.50) 


(ror oa), 
Pr (2—a) 


etl g(x) — g(a). 


We have proved that 


(fog) og 
| (Deo) I] s | ry tls (g(x) - g(a)!" 


(roa ea 
7 rQ2-a) 


vole"! (g (b) — g(a), Vx € [a, DI. 6350) 


In particular it holds 


(Dds 


fjwst, 02a<1, (6.3.52) 
and 

(D°,.,f) (0) = (D&,.,f) @) =0, Vx,y € [a,b]; 0<a<1. (63.53) 
By (6.3.42) we derive 


f@)- f= POs” eee 


i=2 


(Dig f) OD + 


1 “ 
rae / (G(x) — 9 QVM gf (DEL) Odt, 6.3.54) 


for any x > y: x, y € [a, b];0 < a < 1, also it holds 
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fM-fO=>5 oS (Ded) G+ (6.3.55) 


i=2 
1 y 
maapa | (g(y)-g a) ea cae oa (t) (oC) (t) dt, 


foranyy >x:x,y€[a,b];0<a<l. 
We define also the following X-valued linear operator 


(Az (f))  y) = 


n 5 

(g(x)—g(y))'*! ia (n+1)a (gx)-gQ) "De! 
ae TGatl) (Di, jf) ) + (DY fs fe x) F@tbary ? x > ys 
i= 
n 3 

(gQ)—g(a))'*! ia (n+l (G@)-ga)y "re! 
>. Tdatl) (oe jf) @) + (DI +39 fo )“Taaperp —? Yo 4: 
i= 
f' (x), when x = y, 

(6.3.56) 


Vx,y €[a,b];0<a <1. 
We may assume that (see [13], p. 3) 


II(Az (f)) @ x) — (42 (f)) DIL = [FG - FO) (6.3.57) 
< O* |g(x)—g(y)|. Vx, y € [a, DI); 
with ®* > 0. 
We estimate and have 
(i) case of x > y: 
If @) -— f &) — (42 f/f) & y) (gO) -—g OI = (6.3.58) 


= eee eee = (nt+l)a-1 _ (n+ Da _ 
pcesrone (9) — 9) g (t) (a ) (t) dt 


(nt la (g(x) —g Qtr 
jo f) (x) rine Dat D 


(by [1], p. 426, Theorem 11.43) 


f 
~ T(n+Da) 


fa tds ((pet>*f) @= (pi) (x)) dt | 
(6.3.59) 
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(by [9)) 
1 
———— 
r(a+t+ 1a) 


| “ala — gy?" g © | (DEEPA) © (DELPF) Col at 


(we assume here that 


[(oeersyin (ose 


Meer) CO] <erlg@-g@l, 


x >t >y; pi > 0) 


(6.3.60) 
Vt,x,y € [a,b]: 


= eee (9 (x) — g t)) +P"! gf (D) (g ®) — g(t) dt = 


SEE a EERE ED = (n+la 7s 
CEST sf (g (*) — 9) g (t)dt = 


p1 (9 (x) — g (yyy @t DoH? (6.3.61) 
T(mt+ la) (m+Dat+l © 7 


We have proved that 


Ifo) -—fO)-A2)@ VG) -I ODI S 


pi (g(x) —g Qyyerdor! (6.3.62) 
ri@+)e) (@+herI) ~* 7 


Vx,yela,bl:x>y0<a< 1. 
(ii) case of y>x: 


If @) -— fy) -— A2) @ 9G @) - 9g ODI = (6.3.63) 
If) — f@) — (Arf) @ GO) - 9g @)DI = 


1 y 
| maape giv) —g@) er"! 79 () (oes) (t) dt— 


(n+l)a g (y) = 9 ayer 
(Des f)( ) Tint Dat) 


1 
= Tne 
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| | “0 — 9 OP GD (Deere r) = (DEEP*F) 0) at 


a 
~T(n+1)a) 


; 
[ oor-genm a @ | (OEP) © = (BAL) On| ar 3.64 
(we assume here that 
(OA O- (BEA) O]| <mlgO@-9O)1, 63.65) 
Vt,y,x € [a,b]: y>t>x; po > 0) 
P2 / _ (nt+la-1 71 = 
=Tmabo (+a) | (9 (y) — 9 @)) I HG) —g@) dt 


p2 (@G)-—g@)or 


(n+ 1a) (n+latl) (6.3.66) 


We have proved that 


If @) — fy) — (427) @ ) G@) - 9g ODI S 


pr (9 (y) — g(x) Order! 


T(a+la) (a+lat+l) 


’ 


Vx,yela,bl:y>x0<a< 1. 


Conclusion 6.7 Set p = max (1, p2). Then 


If) — f(y) — (A2(f)) & yg) -gO)I S 
p Ig (x) — g(Q)|OrPot 


; (6.3.67) 
Ane AVES Ege a 


Vx,yéla,b]l;0<a<l. 


(Notice (6.3.67) is trivially true when x = y.) 


One may assume that 
p 


ee (6.3.68) 
C((n-+ 1a) 


Now based on (6.3.57) and (6.3.67), we can apply our numerical methods presented 
in this chapter to solve f (x) = 0. 
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Chapter 7 
Iterative Procedures for Solving Equations 
in Abstract Fractional Calculus 


The objective in this study is to use generalized iterative procedures in order to 
approximate solutions of an equation on a Banach space setting. In particular, we 
present a semi-local convergence analysis for these methods. Some applications are 
suggested including Banach space valued functions of fractional calculus, where all 
integrals are of Bochner-type. It follows [5]. 


7.1 Introduction 


Sections 7.1—7.3 are prerequisites for Sect. 7.4. 

Let B,, Bz stand for Banach spaces and let Q stand for an open subset of By. 
Let also U (z, p) := {u € By: ||u —z|| < p} and let U (z, p) stand for the closure of 
U (z, p). 

Numerous problems in Computational Sciences, Engineering, Mathematical 
Chemistry, Mathematical Physics, Mathematical Economics and other disciplines 
can be brought in a form like 

F (x) =0 (7.1.1) 


using Mathematical Modeling [1-17], where F : 2 — By is a continuous operator. 
The solution x* of Eq.(7.1.1) is sought in closed form, but this is attainable only 
in special cases. That explains why most solution methods for such equations are 
usually iterative. There is a plethora of iterative methods for solving Eq. (7.1.1). We 
can divide these methods in two categories. 

Explicit Methods [6, 7, 11, 15, 16]: Newton’s method 


Kat =n — F’ (tn) ' F Gn). (7.1.2) 

Secant method: ‘ 
Xn+1 = Xn — eee Xn F | F (Xn) , (7.1.3) 
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where [-, -; F'] denotes a divided difference of order one on Q x  [7, 14, 15]. 
Newton-like method: 
Xnp1 = Xn — Ey F (on), (7.1.4) 


where E, = E(F)(x,) and E : Q + L(B,, Bo) the space of bounded linear 
operators from B, into By. Other explicit methods can be found in [7, 11, 14, 15] 
and the references there in. 

Implicit Methods [6, 9, 11, 16]: 


F (Xn) =F An (X%n+1 _ Xn) =0 (7.1.5) 
Xnt1 = Xn — Ay F (Xn) ’ (7.1.6) 


where A, = A (%n41,%n) = A(F) (X41, %) and A: Q x Q > L(B, Bo). We 
also denote A (F’) (u, v) = A (F) (vu) = A (u), if u = v for eachu,v € Q. 

There is a plethora on local as well as semi-local convergence results for explicit 
methods [1—7, 9-16]. However, the research on the convergence of implicit methods 
has received little attention. Authors, usually consider the fixed point problem 


P, (x) =x, (7.1.7) 
where 
P, (x) =x + F (z)+A(x,z) @—2) (7.1.8) 
or 
P, (x) =z— A(x, 2)! F(z) (7.1.9) 


for methods (7.1.5) and (7.1.6), respectively, where z € Q is given. If P is a con- 
traction operator mapping a closed set into itself, then according to the contraction 
mapping principle [11, 15, 16], P, has a fixed point x* which can be found using the 
method of successive substitutions or Picard’s method [16] defined for each fixed n 
by 
Vk+1,n = Py, (Yn) > YOn = Xn, Xn41 = lim Ykyn- (7.1.10) 
k—+00 


Let us also consider the analogous explicit methods 


F (Xn) +A (Xn, Xn) (Xn+1 = Xn) =0 (7.1.11) 
ics, AGE) Fy (7.1.12) 
F (Xn) + A (in, Xi) (Xn41 — Xn) =0 (7.1.13) 


and 
Xnt = Xm — A Wns n=1) | F Gp). (7.1.14) 
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In this chapter in Sect. 7.2, we study the semi-local convergence of method (7.1.5) 
and method (7.1.6). Section7.3 contains the semi-local convergence of method 
(7.1.11), method (7.1.12), method (7.1.13) and method (7.1.14). Some applications 


to Abstract Fractional Calculus are suggested in Sect. 7.4 on a certain Banach space 
valued functions, where all the integrals are of Bochner-type [8, 14]. 


7.2 Semi-local Convergence for Implicit Methods 


The semi-local convergence analysis of method (7.1.6) that follows is based on the 
conditions (#7): 
(h;) F: QC B, > By is continuous and A (F) (x, y) € £(B), Bo) for each 
(x, y)EQxXQ, 
(hy) There exist ] > 0 and Qo C B, such that A (F) (x, yy} € £L(Bo, B,) for 
each (x, y) € Qo X Qo and 
JA) @, yy" se. 


Set Q]) = QA Qo. 
(h3) There exist real numbers a1, a2, a3 satisfying 


0 < a <q; and0 < a3 <1 
such that for each x, y € QQ; 
IF @) — F(y)-A(F)@, y)@—y)Il S 
1 (Sle = yl + a9 ly — oll +08) lle — yl. 
(h4) For each x € Qo there exists y € Qo such that 
y=x—Aly,x) F(a). 
(hs) For xo € Qo and x; € QQ satisfying (h4) there exists 7 > 0 such that 
||A (F) G1, x0)! F (x0) || < 2. 


(tg) h = aun < 3 (1— 03)’. 
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and 
(hz) U (x0, t*) C Qo, where 


Then, we can show the following semi-local convergence result for method (7.1.6) 
under the preceding notation and conditions (#7). 


Theorem 7.1 Suppose that the conditions (H) are satisfied. Then, sequence {xn} 
generated by method (7.1.6) starting at xy € Q is well defined in U (xo, t*), remains 
in U (xo, t*) for each n = 0, 1,2, ... and converges to a solution x* € U (xo, t*) 
of equation F (x) = 0. Moreover, provided that (h3) holds with A(F)(z,y) replacing 
A(F)(x,y) for each z € Q, , if a; # 0, the equation F (x) = 0 has a unique solution 
x*inU , where 
~ [U0 f)NQo, fh = 5 (1 — 03)" 
~ [| U Go, ON Qo, fh < 4-93) 


and, if a, = 0, the solution x* is unique in U (xo, ul ), where t** = 


l-a3 
1—a3+4/(1—a3)?—2h 


ay 


Proof Case a; # 0. Let us define scalar function g on R by g(t) = wf? _ 
(1 — a3) t + 7 and majorizing sequence {t,} by 


to = 0, te = t%&-1 + g (t%&e-1) foreachk = 1,2,.... (7.2.1) 
It follows from (h¢) that function g has two positive roots t* and t**, t* < r**, and 
ty < th4 . That is, sequence {t,} converges to ft”. 
(a) Using mathematical induction on k, we shall show that 


Xx41 — Xell S ee — te. (7.2.2) 


Estimate (7.2.2) holds for k = 0 by (hs) and (7.2.1), since ||x; — xo|| < 7 = t — fo. 
Suppose that for 1 <m<k 


Xm — Xm—11l S tn — tm—1. (7.2.3) 


Them, we get ||x, — xol| < t¢ — to = te < t* and A (xg, xx_1) is invertible by (h2). 
We can write by method (7.1.6) 


Xep1 — Xe = —Ag! (F (xn) — F (e-1) — Ant (te — Xe-1))- (7.2.4) 


In view of the induction hypothesis (7.2.3), (h2), (43), (ha), (7.2.1) and (7.2.4), we 
get in turn that 
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-l -l1 
lse-1 — xell = Ag’ F Gx) |] = Ag! OF Ge) — F @e-1) — Ani Ox — 2-1) 


< Ag! MF Ge) — F Gee) — Acai Ore — xe-)I) < 


= ay 

I (> IlX~¢ — Xe—1l| + 2 |lXe-1 — Xoll + as) llxz — Xx_11| < (7.2.5) 
QA) 2 
> (th — tei)" + O2 (te — th-1) h-1 + 03 (he — te-1) = 


a 
i (te — th-1)? + O12 (te — te—-1) te + 03 (te — the-1) — (te — th-1) + 8 (e-1) = 


& (tk) — (1 — O22) (te — te-1) ke-1 < 
8 (tk) = thei — th, (7.2.6) 


which completes the induction for estimate (7.2.2). 
That is, we have for any k 


X41 — Xell S they — tk (7.2.7) 


and 
Ixx —xoll <& < 0%. (7.2.8) 


It follows by (7.2.7) and (7.2.8) that {x;} is a complete sequence in a Banach space 


B, and as such it converges to some x* € U (xo, t*) (since U (xo, t*) is a closed 
set). By letting k — +00, using (h;) and (ho), we get 1 Jim || F (xz) || = 0, so 
—> +00 


F(x*)=0. | 
Let x** € U be such that F (x**) = 0. We shall show by induction that 


\|x** — x4|| < ¢* — % for each k = 0, 1,2,.... (7.2.9) 


Estimate (7.2.9) holds for k = 0 by the definition of x** and U. Suppose that 
l|x** — xg|| < t* — t. Then, as in (7.2.5), we obtain in turn that 


||x** - Xe = ||x** — x, + Ay F (x) — Ay F (x) | = 


Ac! (Ae — ax) + F Ga) — F (x) 


IA 


Jar" LF (0) - Fo) — Ax (0** =) 


IA 


(5 \|x** — xg|] + a2 lax — xoll + as) \|x** — x,|| < 
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(> (t* _ tk) + ot, + a3) (t* = tk) = 


= (r*) + 5 (th) — ay tyt® + a9 (t* — t) t +03 (t* — &) = 


—n+ (1 — a3) t+ st —aytt*® + artt*® — art? + a3t* — azty 
Spa dey, (7.2.10) 


which completes the induction for estimate (7.2.9). Hence, A lim x, = x**. But we 
—>+00 


showed that lim x, = x*, so x** = x*. 


k—>+00 


Case a; = 0. Then, we have by (43) that a2 = 0 and estimate (7.2.5) gives 
lxea1 — Kell S03 Xe — Xe-ill S ++ S04 [lar — xoll < 4 (7.2.11) 


and 
Wxea1 — Xoll < [ea — xell + lle — xg-1 |] +... + [ler — xl 
1—a; 
z 3 1) 


n< : (7.2.12) 
1- a3 1- a3 


Then, as in the previous case it follows from (7.2.11) and (7.2.12) that 


i 
1—a, 


Ilc4i — Xx] < akn, (7.2.13) 


= l-a; 


so sequence {x;,} 1s complete and x* solves equation F (x) = 0. Finally, the unique- 
ness part follows from (7.2.10) for a; = az = 0, since 


Je = xers | <a x — ae] <8" |x — x0] < o5"" 5 J. (244) 
= fis 
which shows again that lim x, = x**. a 
k->+00 


Remark 7.2 (1) Condition (2) can be incorporated in (3) as follows 
(hs) There exist real numbers @, @2, @3 satisfying 0 < @ < a and0 < a3 < | 
such that for each x, y € Q 
|A@.») "7 LF @)- FO)- AG.) @-y)]| < 


(1/2) lx = yl + lly — xoll +) x yl. 
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Then, (43) can replace (2) and (h3) in Theorem 7.1 for a) = @, @2 = a», 
Q@3 = a3 and Qo = Q. Moreover, notice that @| < a,j, Q@ < a, and a3 < a3, 
which play a role in the sufficient convergence criterion (hg), error bounds and the 
precision of t* and t**. Condition (h3) is of Mysowksii-type [11]. 
(2) Suppose that there exist Jo > 0, a4 > O and L € £(B), By) with L~! 
L (Bp, By) such that || L~!|| < 19! 


||A (F) (x, y) — L|| < a4 for each x, y € Q 


and 
as = Iau <i. 


Then, it follows from the Banach lemma on invertible operators [6, 10, 11, 15, 16] 
and 
[2° 1A) @, y) - Lil s Iga = a5 <1 


that A (F) (x, y)~! € £L (Bo, B,). Seti“! = 
condition (42) is implied, so it can be ek fon the condiGons (A). 

(3) Clearly method (7.1.5) converges under the conditions (1), since (7.1.6) 
implies (7.1.5). 

(4) Let R > O and define Ro = sup {t € [0, R): U (x, Ro) € D}. Set Qo = 
U (xo, Ro). Condition (h3) can be extended, if the additional term ap ||x — xo|l is 
inserted inside the parenthesis at the right hand side for some az > O . Then, the 
conclusions of Theorem 7.1 hold in this more general setting, provided that a3; = 
az Ro + a3 replaces a3 in conditions (h6) and (7). 

(5) Concerning the solvability of Eqs. (7.1.6) (or (7.1.5)), we wanted to leave 
condition (44) as uncluttered as possible in conditions (H). We did this because in 
practice these equations may be solvable in a way other than using the contraction 
mapping principle already mentioned earlier. 


Next, we show the solvability of method (7.1.5) using a stronger version of the 
contraction mapping principle and based on the conditions (C) : 
(c1) = (11). 
(c2) There exist yo € [0, 1), y1 € [0, +00), y2 € [0, 1), xo € Q such that for each 
x,y,ZzEQ 
I+ A@,z)-A(Q, ZI < 10, 


A @, 2) -A(y, DI < villa — yIl 


= 2 Xo — zl| for xo Az 
| F(z) + A (x0, z) (xo — ZI < Le Gell tory a 
(c3) 

yo +1 |lxoll + y2 < 1 for y2 40, 
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Yo + V1 IIlxoll < 1 for 7 = 0, 


= 2 
IF Goll < (— (70+ % Ilxoll)) foray 20, 


v1 


yo < | fory, =0 


and 
(c4) U (xo, r) € &, where 


F 1- 
IF (xo) I 23 (Yo + V1 IIxoll) foray eu 
1— (+7 Ilxoll) a 
F 
IF (xo) | Seep, 
1—% 


1—(%+% Ilxoll) 
r< 
onl 


for z= x0, y1 #0. 


Theorem 7.3 Suppose that the conditions (C) are satisfied. Then, for each n = 
0,1,2,... Eg.(7.1.5) is unique solvable. Moreover, if A;! € L(B), B), then 
Eq. (7.1.6) is also uniquely solvable for eachn = 0, 1, 2, ... 


Proof We base the proof on the contraction mapping principle. Let x, y € U (xo, r). 
Then, using (7.1.8) we have in turn by (c2) that 


I|P: ) — Pz (Yl = Il + A, 2) — AQ, 2) @ — y) — AG, I-A, Dall 
< ||P+AQ,2)-AQV, 2) lle — y+ IAG, 2) -— AG, ZI Mell 


< yo llx — yll +1 Cz — xoll + Ilxol) Ile — yIl 


< p(llx — xoll) Il — yll, (7.2.15) 
where 
p=} i ; < a uy nh aa (7.2.16) 
Notice that y (tf) € [0, 1) fort € [0, r] by the choice of r in (c4). 
We also have that 
I| Pz ©) — xoll < IP; @&) — Pz Go) Il + Il Pz 0) — xoll - (7.2.17) 


If z = xo in (7.2.17), then we get by (c3) , (ca) and (7.2.15) that 


|| Px) (&) — xol] < v (lx — xoll) lx — xoll + IF Go) 
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< (+7 Ilxoll) 7 + IF Goll <r. (7.2.18) 


The existence of x, € U (xo, 1r) solving (7.1.5) for n = 0 is now established by the 
contraction mapping principle, (7.2.15) and (7.2.18). 

Moreover, if z # xo, the last condition in (c3), (c3) , (ca) and (7.2.17) give instead 
of (7.2.18) that 


||P; (x) — xoll < (lle — xoll) Ix — xoll +72 Ix — oll 
< (+1 Ilxoll +72) r <r. (7.2.19) 


Then, again by (7.2.15), (7.2.19) and the contraction mapping principle, we guarantee 
the unique solvability of Eq. (7.1.5) and the existence of a unique sequence {x,} for 
each n = 0, 1, 2,.... Finally, Eq. (7.1.6) is also uniquely solvable by the preceding 
proof and the condition A;! € £ (Bo, Bi). | 


Remark 7.4 (a) The gamma conditions can be weakened, if 7; are replaced by func- 
tions 7; (t), i = 0,1, 2,3. Then, 7; will appear as 4; (||x — xo||) and 7; (7) in the 
conditions (C) . 

(b) Sections 7.2 and 7.3 have an interest independent of Sect. 7.4. However, the 
results especially of Theorem 7.1 can apply in Abstract Fractional Calculus as sug- 
gested in Sect. 7.4. 


7.3 Semi-local Convergence for Explicit Methods 


Theorem 7.1 is general enough so it can be used to study the semi-local convergence of 
method (7.1.11), method (7.1.12), method (7.1.13) and method (7.1.14). In particular, 
for the study of method (7.1.12) (and consequently method (7.1.11)), we use the 
conditions (H’) ; 

(h\) F : QC B, — Bp» is continuous and A (F) (x, y) € £(B,, Bo) for each 
xEQ. 

(h,) There exist / > 0 and Qo C B, such that A (F) (x, x)! € £L (Bo, B)) and 


|A CF) @,x) Tt] <r. 


Set Qy = 2M Qo. 
(h}) There exist real numbers 71, a2, 73 satisfying 


O0<a. < y and0 < 473 


114 7 Iterative Procedures for Solving Equations in Abstract Fractional Calculus 
such that for each x, y € Q 
IF @)-FO)-AP)(,y) @— ys 
V1 
1 (} le yll + a2 lly — oll +78) [lx — Il. 
(h,) For each x, y € Q; and some 74 > 0, ys => 0 


JA (x,y) —A(y, WIL Sly 


or 

|A (x,y) — A(y, y)Il Sls Ilx — yl. 
Set aj = ¥ +5 and a3 = 73 + 4, if the second inequation holds or a, = 7, and 
a3 = 73 + ya, if the first inequation holds. Further, suppose 0 < a3 < 1. 


(h) There exist x9 € Qo and 7 > O such that A (F) (x0, xo)! € £ (Bo, Bi) and 


|| A (F) (x0, x0)! F (xo) || < n- 


(he) = (he) 
and 
(h5) = (h7). 


Then, we can show the following semi-local convergence of method (7.1.12) using 
the conditions (H ’) and the preceding notation. 


Proposition 7.5 Suppose that the conditions (H') are satisfied. Then, sequence {x;} 
generated by method (7.1.12) starting at Xx) € Q is well defined in U (xo, t*), remains 
in U (xo, t*) for each n = 0, 1,2, ... and converges to a solution x* € U (x0, t*) 
of equation F (x) = 0. Moreover, if a, # 0, the equation F (x) = 0 has a unique 
solution x* in Uy, where 


_ [U Got") Qo, ifh = 3 (1-3)? 


Y= LU Go, t) DQ, fh <4 0-03) 


and, if a = 0, the solution x* is unique in U (xo, 7) where t* and t** are given 


in Theorem 7.1. 


Proof Use in the proof of Theorem 7.1 instead of estimate (7.2.5) the analogous 
estimate 


IF (xg) |] = WF (xe) — F e-1) — A (e-1, Xe—1) (Xe — Xe) I = 


[LF Ga) — F Cre) — A (re, e-1) Ce — xe-1)] + 
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(A (xq, Xe—-1) — A (Xp 1, Xe-1)) (Xk = XK-1) I 


onl 
<1 (2 fixe — seal + ac lxi—1 — oll + 79) lx — X41 + 
A (cK, Xe—-1) — A (K-15 Xe) NR — Xe-11] < 
Qa] 2 
l = (th — th-1)” + Q2 (te — th—-1) th-1 + 03 (th — Te-1)) 5 
where we used again that ||x, — xx-1|| < te — te-1, |lxe—-1 — Xoll < t%{-1 and the 


condition (h',) . 


Remark 7.6 Comments similar to Remark 7.2 (1)—(3) can follow but for method 
(7.1.11) and method (7.1.12) instead of method (7.1.5) and method (7.1.6), respec- 
tively. 


Similarly, for method (7.1.13) and method (7.1.14), we use the conditions (H ” ) : 
(hi) = Gi) 
hy) = (hz) 
hs) There exist real numbers a1, a2, 73 satisfying 
0 <a) < a, and0 < 73 
such that for each x, y € Q) 
IP @Q=F ()=—AGI Gy) @=yIs 
Qa) 
1 (S ll — yll + a2 lly — xoll +8) x — yl. 
(h{) For each x, y, z € Qy and some 73 > 0 
|A (zy) —-A Gy, x)I] < 163. 
Set a3 = 73 + 63 and further suppose 0 < a3 < 1. 
(hg) There exist x_; € Q, x9 € Q and 7 > O such that A (F) (x9, x1)! €E 


L (B2, B,) and 
A (F) (xo, x-1)7! F (x0) || < 2. 


(hg) = (he) 
and 
(ht) = (hy). 


Then, we can present the following semi-local convergence of method (7.1.14) 
using the conditions (H 7 ) and the preceding notation. 


Proposition 7.7 Suppose that the conditions (H") are satisfied. Then, sequence {x;} 
generated by method (7.1.14) starting at Xx) € Q is well defined in U (xo, t*), remains 
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in U (xo, t*) for each n = 0,1, 2, ... and converges to a solution x* € U (xo, t*) 
of equation F (x) = 0. Moreover, if a, # 0, the equation F (x) = 0 has a unique 
solution x* in U, where 


ga | U Got) MQ, ifh = 7 — 03)" 
~ [| U o, t*) NQo, ifh < $1 — 03) 


and, if a, = 0, the solution x* is unique in U (x0, 4), where t* and t** are given 


in Theorem 7.1. 


Proof As in Proposition 7.5, use in the proof of Theorem 7.1 instead of estimate 
(7.2.5) the analogous estimate 


IF (xe) Il = 
I|F (xk) — F Gre-1) — A (tk, Xe-1) (Xk — Xe-1) 
+ (A (Ce, Xk-1) — A Cre-1, Xk-2)) Oe — Xe-DII S 
I|F e) — F Oe—-1) — A Ok, Xe—-1) Oe — Xe-1) | + 


A (xx, Xe-1) — A (Xp—1, Xe—2) I] Nove — Xe—1Ih 


Qa) 
<1 (FS Mx — xual] + a tit — oll +79) tk — xe-ill +15 Ulxe — 24-1 


QA 2 
<l (> (te — th_1)7 + Ol (te = th—-1) te-1 + 03 (th — n-1)) ; 
where we used again that ||x, — xx-1|| < te — te-1, |lae—-1 — Xoll < t%{-1 and the 
condition (11) . 


Remark 7.8 Comments similar to Remark 7.2 (1)—(3) can follow but for method 
(7.1.13) and method (7.1.14) instead of method (7.1.5) and method (7.1.6), respec- 
tively. 


In particular, the results of Sect.7.4 connect to the crucial conditions (hs) and 
(hi) as follows: According to the definition of method (7.1.12) and using (7.4.1), 
(7.4.8) and (7.4.15), we can write 

If OI = IF O-FO-AO)@-yIl= 


I[Fo-fO-fMe-y]4+[f &)-AM] a -»y| 


IA 


Ifm-fM-fMOe@-ylt]/ O&-40)] Ie -yI 


[Fle 
2 


IA 


1 
2 


y+ 2 IF" loo Ly = 40l Ix — yl. 
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It follows that (h) and (h/,) hold for 


_ I'l. 
1 ’ 


IP" Te. 


Pa Soa ai 
and 73 = 74 = Ys = 0. Notice also that 0 < 73 and0 < az < 1, since l—a < 2—a 
and 73 = 0. 


7.4 Applications to Abstract Fractional Calculus 


Here we deal with Banach space (X, ||-||) valued functions f of real domain [c, d]. 
All integrals here are of Bochner-type, see [14]. The derivatives of f are defined 
similarly to numerical ones, see [17], pp. 83-86 and p. 93. 

We want to solve numerically f (x) = 0. 

Let 0 < a < 1, hence [a] = 1, where [-] is the ceiling of the number. Let also 
c<a<b<d,and f € C? ({c,d], X), with f” ¢ 0. Clearly we have (see [12], p 
3) 

If @-Ff Os [P14 — yl Vey € led]. (7.4.1) 


(I) The X-valued left Caputo fractional derivative (see [4]) of f of order a € (0, 1), 
anchored at a, is defined as follows: 


a _ 1 7 -—a ¢l 
(Def) (x) = aaa | (x —t) “ f' @dt, Vx €[a,d], (7.4.2) 


while (De ust ) (x) = 0, for c < x <a, where I is the gamma function. 
Next we consider a < a* < b, and x € [a*, b], also xo € (c, a). 
We define the function 


lr (2-a) ( 
(x =p)? 


Notice that A; (a) is undefined. 
We see that 


Ay (x) = D2, f) (x), Vx € [a*, bd]. (7.4.3) 


lr (2-—a) 


— gyice Paef) @)— £) 


Ar @ - f'@)| = |= (7.4.4) 


pee) 7 pF 2=D @ a | 
(x —a)' @T(1— a: ee Oe ea rool = 


re 


(x —a)'~ 


1 . —a yl 1 * —a zl 
saa « t) “ f @adt aaa « t) “ f (x)dt 
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(by [1] p. 246, Theorem 11.43) 


= ee [ (x—1)*(f'@) — f' (x) dt (7.4.5) 
(x — a) . a 
(by [8]) 
<i mc =f xn)" foO-f (x)|) dr "< 
C=O bes [op - 0) [7 "le ” 
@ —a)=4 | (x —t) “(x-ft)dt= aye =f (x —t) dt = 
d-o)|f’|,,@-@** _G-a)y., 
(x =<)? 2—a ~ (2Q- 2 If ll ¢ (7.4.6) 
We have proved that 


_ 
|A1@)- fs (5— ) IF" @-9 s (5) I" ln 0-9. 
2-—a 
(7.4.7) 
Vx € [a*, b]. 
In particular, it holds that 


i= 
|41.@)-f'@| < (; — “) IF" loo & — 0). (7.4.8) 


where x9 € (c, a), V x € [a*, b]. 
(I) The X-valued right Caputo fractional derivative (see [3]) of f 
of order a € (0, 1), anchored at b, is defined as follows: 


=| b 
(De_f) (x) = Pda) / (t—x)° f' (dt, Vx €[c, b], (7.4.9) 


while (D¢_ f) (x) = 0, ford > x > b. 
Next considera < b* <b,andx € [a, b*], also xo € (b, d). 
We define the function 


lr (2-—a) a : 
A2 (x) := ~@-nre (Dj_f) (x), Vxe [a, b | : (7.4.10) 
Notice that A> (b) is undefined. 
We see that 
r(2 
| A2@) — f’@]| = | DAF @ 
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Pr (2-a) 1 ee 24 
Ie / (¢—x)" f @)dt— f' @&) (7.4.11) 
Pr (2—a) _ rQ-a) b-x'* 
(b— x) ae sf C—x" i Odt— 7 ara Pea we 
r(2—a) : 
(b — x)!- Qa rda= a: (t x) fot i (t) dt a i (t x) f' (x) dt 
(7.4.12) 
r2—a) = 
=z)" are a) If @—xy“ ae Qa (x)) dt] < 
, (7.4.1) 
i Gx O- sf war "S 
—x) x 
l= " aa 
“orl " grees oy =f x) dt = 
x 
(7.4.13) 
1=— u _ 2-a —_ " 
(=o) 17h. @=9" _A-OlF"hs gy, 
(b—x)° 2-—a (2— a) 
We have proved that 
|| Az @&) — f’ @| <(= alae o-9<(— alae (b—a), 
(7.4.14) 
Vx ela, b*]. 
In particular, it holds that 
|A2@) — f’@I} < (= alae (xo — x), (7.4.15) 


where xo € (b, d), V x € [a, b*]. 
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Chapter 8 
Approximate Solutions of Equations in 
Abstract g-Fractional Calculus 


The novelty of this chapter is the design of suitable iterative methods for generating 
a sequence approximating solutions of equations on Banach spaces. Applications of 
the semi-local convergence are suggested including Banach space valued functions 
of fractional calculus, where all integrals are of Bochner-type. It follows [6]. 


8.1 Introduction 


Sections 8.1—8.2 are prerequisites for Sect. 8.3. 

Let B,, By stand for Banach spaces and let Q stand for an open subset of By. 
Let also U (z, €) := {u € By: ||u — z|| < €} and let U (z, &) stand for the closure of 
U (z, &). 

Many problems in Computational Sciences, Engineering, Mathematical Chem- 
istry, Mathematical Physics, Mathematical Economics and other disciplines can writ- 
ten like 

F(x) =0 (8.1.1) 


using Mathematical Modeling [1-18], where F : 82 — Bp is a continuous operator. 
The solution x* of Eq. (8.1.1) is needed in closed form. However, this is achieved 
only in special cases. That explains why most solution methods for such equations 
are usually iterative. There is a plethora of iterative methods for solving Eq. (8.1.1), 
more the [2, 7, 8, 10-14, 16, 17]. 

Newton’s method [7, 8, 12, 16, 17]: 


te =a I Gay P Ge) (8.1.2) 

Secant method: : 
Xn4+1 = Xn — [%n—15 Xn3 Fl F (Xn) ; (8.1.3) 
© Springer International Publishing AG 2018 121 
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where [-, -; F'] denotes a divided difference of order one on Q x  [8, 16, 17]. 
Newton-like method: 
Xnt1 = Xn — Ey F (x), (8.1.4) 


where E, = E(F)(x,) and E : Q — L(B,, Bo) the space of bounded linear 
operators from B, into Bz. Other methods can be found in [8, 12, 16, 17] and the 
references therein. 
In the present study we consider the new method defined for each n = 0, 1, 2, ... 
by 
Xn+1 = G (Xn) 


G (Xn41) = G (tn) — Ay F Gn), (8.1.5) 


where xo € Q is an initial point, G : Bz; — Q (B3 a Banach space), A, = 
A(F) (n41,%n) = A (Xn41,%) and A: Q x Q > L(B), Bo). Method (8.1.5) 
generates a sequence which we shall show converges to x* under some Lipschitz- 
type conditions (to be precised in Sect. 8.2). Although method (8.1.5) (and Sect. 8.2) 
is of independent interest, it is nevertheless designed especially to be used in g- 
Abstract Fractional Calculus (to be precised in Sect.8.3). As far as we know such 
iterative methods have not yet appeared in connection to solve equations in Abstract 
Fractional Calculus. 

In this chapter we present the semi-local convergence of method (8.1.5) in 
Sect.8.2. Some applications to Abstract g-Fractional Calculus are suggested in 
Sect.8.3 on a certain Banach space valued functions, where all the integrals are 
of Bochner-type [9, 15]. 


8.2 Semi-local Convergence Analysis 


We present the semi-local convergence analysis of method (8.1.5) using conditions 
(M): 
(m,) F : Q Cc By, — Bp» is continuous, G : B3 — Q is continuous and 
A(x, y) € £(B), Bo) for each (x, y) € Q x Q. 
(m) There exist 3 > 0 and Qo C B; such that A (x, y)~! € L (Bo, B;) for each 
(x, y) € Qo x Qo and 
JAG yy" <8". 


Set Qy =QN Qo. 


(m3) There exists a continuous and nondecreasing function w : [0, +00)? > 
[0, +00) such that for each x, y € Q) 


IF @)-FO)-AG, y)(G@)-GO))I < 


Bw (lx — yl. lle — xoll ly — xoll) IG @) -— GOI. 
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(m4) There exists a continuous and nondecreasing function wo : [0,-+oo) > 
[0, +00) such that for each x € Qy 


|G (x) — G (ro)Il_ S Yo (lx — xoll) lx — xoll - 
(ms) For x9 € Qo and x; = G (xq) € Qo there exists 7 > 0 such that 
||A G1, x0)! F (xo) || < n- 
(m6) There exists s > 0 such that 
v(m, s,s) <1, 
Yo (8) < 1 


and 


(Go) =x <s2—"— 
— 40 


where go = w (7), 5, 8). 

(m7) U (x0, 5) CQ. 

Next, we present the semi-local convergence analysis for method (8.1.5) using 
the conditions () and the preceding notation. 


Theorem 8.1 Assume that the conditions (M) hold. Then, sequence {x,} generated 
by method (8.1.5) starting at xy € (2 is well defined in U (xq, 5), remains in U (xo, 5) 
for each n = 0,1,2,... and converges to a solution x* € U (xo, 5) of equation 
F (x) = 0. The limit point x* is the unique solution of equation F (x) = 0 in 
U (xo, S) : 


Proof By the definition of s and (ms), we have x; € U (xo, 5). The proof is based on 
mathematical induction on k. Suppose that ||x, — xx—-1|| < qn and ||x, — xol| <s. 

We get by (8.1.5), (72) — (ms) in turn that 

IG Gs) — G oll = Ag’ F Go|] = 
| Ac! F Ge) — F Ge-v) — Aci (G On) — G G1) | 
< Ag! | WF Ge) — F Gre) — Aci (G Oe) — G OI S 
B-' Bw (xe — xxi) 5 xe—1 — Xolls Ive — xoll) IG (xe) — G Il < 
W (n, 8,8) IG (xe) — G (Xe) = 90 IG Gx) — G Ox-DIl < 96 lla — xoll < 


(8.2.1) 
gon 
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and by (m6) 
Ilxx+1 — Xoll = IG (xe) — Xoll S IG Gx) — G Co) Il + IG (0) — oll 
S Yo (lle — xoll) Ilex — xoll + IG Go) — xoll 
< Yo (s) 5 + [IG (xo) — xoll < ». 


The induction is completed. Moreover, we have by (8.2.1) that for m = 0, 1, 2, ... 


[een — all S gb, 

— 90 

It follows from the preceding inequation that sequence {G (x;)} is complete in a 

Banach space B, and as such it converges to some x* € U (xo, 5) (since U (xo, 5) 

is a closed ball). By letting k — +00 in (8.2.1) we get F (x*) = 0. We also get 

by (8.1.5) that G (x*) = x*. To show the uniqueness part, let x** € U (xo, 5) bea 

solution of equation F (x) = 0 and G (x**) = x**. By using (8.1.5), we obtain in 
turn that 


\|x** — G Gea) || = |x — G Ge) + Ap! F Ox) — Ap F (x) | < 
4c" |F (e) — F Gx) — At (6) - Gen)|| s 


B" Babo (x* — xe]. lea — oll. Ilex — xoll) |G (x) — G @)I| < 


’ 


qo |G (x) — G Gx) S40" x — x0 


so lim x, = x**. We have shown that lim x, = x*,sox* = x**. | 
k—>+00 k—>+00 


Remark 8.2 (1) Condition (m2) can become part of condition (m3) by considering 
(m3)’ There exists a continuous and nondecreasing function ¢ : [0, +00)? > 
[0, +00) such that for each x, y € Q 


|A (x, y) | LF @&) — FQ) — A(x, ») (E(x) — G(y))]] < 
y (Ix — yl], lla — xoll, ly — xl) IG @) — GO)IL. 


Notice that 
p (uy, U2, u3) < 1) (uy, U2, u3) 


for each u, > 0, u2 > O and uz > O. Similarly, a function y; can replace y, for the 
uniqueness of the solution part. These replacements are of Mysovskii-type [7, 12, 
16] and influence the weaking of the convergence criterion in (m6), error bounds and 
the precision of s. 
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(2) Suppose that there exist @ > 0, 8; > O and L € L(B,, Bo) with L7! € 
£ (B2, By) such that 
Je] se" 


|A@,y)-Lll sf 


and 


hao. <1 


Then, it follows from the Banach lemma on invertible operators [12], and 
|e" AG») -LI< 6" =h <1 


that A (x, y)~! € £ (Bo, B,). Let B = = Then, under these replacements, condi- 


tion (m2) is implied, therefore it can be dropped from the conditions (M). 


Remark 8.3 Section 8.2 has an interest independent of Sect. 8.3. It is worth noticing 
that the results especially of Theorem 8.1 can apply in Abstract g-Fractional Calculus 
as illustrated in Sect. 8.3. By specializing function 7, we can apply the results of say 
Theorem 8.1 in the examples suggested in Sect. 8.3. In particular for (8.3.33), we 
choose for u; > 0, u2 > 0, u3 > O 


Aut 


W (uy, U2, U3) = AT (@) (a+ 1)’ 


if |g (x) — g (y)| < 1 for each x, y € [a, b]; 


Mes 


W (uy, U2, U3) = BY (a) (a+ 1)’ 


if |g (x) — 9 (y)| S & llx — yl for each x, y € [a, b] and pun = & |b — al; 


AUS 


WY (uy, U2, U3) = BT (a) (a +1)’ 


if |g (x)| < & for each x, y € [a, b] and 3 = 2&3, where a, \ and TI are given in 
Sect. 8.3. 

Other choices of function ~ are also possible. 

Notice that with these choices of function 7 and f = F and g = G, crucial 
condition (m3) is satisfied, which justifies our definition of method (8.1.5). We can 
provide similar choices for the other examples of Sect. 8.3. 
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8.3. Applications to X-valued g-Fractional Calculus 


Here we deal with Banach space (X, ||-||) valued functions f of real domain [a, b]. 
All integrals here are of Bochner-type, see [15]. The derivatives of f are defined 
similarly to numerical ones, see [18], pp. 83-86 and p. 93. 

Here both backgrounds needed come from [5]. 

(1) We need 


Definition 8.4 ({5]) Let a > 0, [a] =n, [-] the ceiling of the number. Let f € 
C" ({a, b], X), where [a, b] C R,and (X, ||-||) isaBanach space. Let g € C! ({a, b]), 
strictly increasing, such that g~! € C" ([g (a), 9 (b)]). 

We define the left generalized g-fractional derivative X-valued of f of order a as 
follows: 


(D¢,.,f) @) = wool (g(x) —g ("0 (fog) GO) dt, 


(8.3.1) 
V x € [a, b]. The last integral is of Bochner type. 


If a ¢ N, by [5], we have that (ps fl) € C ({a, b], X). 


We see that a 
Ces ((F og) o s)) (x) = (D%,,,.f) @) Nx € [a, B). (8.3.2) 
We set 
Dt f (x)= ((F og!)"o 4) (x) € C (la, b], Xn EN, (8.3.3) 


Do. f @) =f), Vxe [a,b]. 


When g = id, then 
Disigl = Disial = Deeds 
the usual left X-valued Caputo fractional derivative, see [4]. 
We need the X-valued left general fractional Taylor’s formula. 


Theorem 8.5 ((5]) Leta > 0,n = [a], and f € C” ([a, b], X), where [a,b] CR 
and (X, ||-||) is a Banach space. Let g € C! ({a, b]), strictly increasing, such that 
“he C" (g(a), g(b))), a < x <b. Then 


— (g(x) = 9 (@))' é 
FO=f@+tD (fog) Gay t 


i=1 


: a-l ws a _ 
es) [ (g(x) — gO)" (0) (Diag f) Oat = 
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n—1 i ; 
e+ > LOO 59-9 G@+ 


i=1 


r (a) (a) 


The remainder of (8.3.4) is a continuous function in x € [a, b]. 


1 gr) 
=a! (g (x) —z)*" ((D2,.,f) og’) @ dz. 
g 
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(8.3.4) 


Here we are going to operate more generally. We consider f € C” ([a, b], X). We 
define the following X-valued left g-fractional derivative of f of order a as follows: 


1 x 
(Dyjigf) @) = Taco / ga-gM" "J @ (fog 


foranya<y<x<b; 


Disgf = (Fog!) 09) @).Vx.y € [a,b], 


and 
Dit = f(x) Vx € [a,b]. 


For a > 0, a €N, by convention we set that 


(Ds 


vugd) @) =O, form < y, Vx,y € [a,b]. 


Similarly, we define 


(De, 7/or= wo!/ G@-g@y"* 19 O (fog 
foranya<x<y<b; 


Disgf = (Fog!) 09) O).V x.y € La, BI, 


and 
Df =f OM y € [a4]. 


For a > 0, a € N, by convention we set that 


(Ds 


Cuigf) (y) = 0, fory < x, Vx, y € [a,b]. 


N gqa))dt, 


(8.3.5) 
(8.3.6) 


(8.3.7) 


(8.3.8) 


“© (g(t) dt, 
(8.3.9) 


(8.3.10) 


(8.3.11) 


(8.3.12) 
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We get that (see [9]) 


a 1 ° n—a-1 ~1)\(@ 
| ols mas [ e@-sore ty o|(tes)” een] at 
(8.3.13) 
(reg) oa, 7 -—a-1l 
= Cee Gaal | (g(x)-g@)"° f Odt = 
| (f . a)” og Pe n-a 
Paoatl (g@@)-g@)y""s 
| (f gy” ae — n-a 
Cra (g (b) — g(a))" °,Vx € [a,b]. (8.3.14) 
That is 
(DF, f) aQ=0, (8.3.15) 
and 
(Dy.gf) &) = (Dei. f) @) = 0,¥ x, y € [a, 5). (8.3.16) 


Thus when a > 0, a ¢ N, both Dy igh Digt € C ([a, b], X), (see [5]). 


Hence by Theorem 8.5 we obtain 


n-1 og (k) 
f@-f-y g y (9 0) 


k=1 


(g(x) — gy) + 


— : a-l ys a 
ra | (9) — gO)" 9' 0) (Dig f) dt. Vx ely, bl], (8.3.17) 


and 


a ee a 
fa-fe=y g y (9 (x)) 


k=1 


(g(y) — g(x) + 


1 : a— / a 
—a (gv) — 91)" o' © (Dig f) Odt.V y € [x, 5], (8.3.18) 
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We define also the following X-valued linear operator 


(Ai (f)) @ Y) = 


n—-l og7! (k) . x a 
> (fog ) OO) (9 (x) — g (yy + + (De, a) (x) fe) a forx > y, 


oq7 x x 
= (fog " (g@)) (g(y) gtx) Hy + ( Deg) (y) WO Be)" for x <= Jj, 


i pn (x), when x = y, 


(8.3.19) 
Vx,y €[a,bl;a>0,n= fa]. 
We may assume that (see [13], p. 3) 
(Ai (f)) @ x) — (Ar (A). DL = LF @ - F ©) (8.3.20) 


(Ff og) G@@) — (fF og) GO) < lg@) —gO)I.V x, y € Ia, BI; 
where ® > 0. 


We estimate and have 
(i) case of x > y: 


If a) - fO)- (41 (f) & GO) - gO) = 


nore (g(x) —g(t))*"' 7g! (1) (D af) @) dt— 


(g (x) — 
T(at+l) 


(Bit) @ (8.3.21) 


(by [1] p. 426, Theorem 11.43) 


1 - ae, a a 
=“T@ | / (g(x) — 9)! gO (Doug f) O — (Doug f) O)) at 
(by [9]) 
1 : a-l vs a Qa 
ae) | (g(x) -— 9) ' 9 @ | (Dif) O — (Doug f) || at (8.3.22) 


(we assume that 


| (Dsof) O — (Doy.gf) CO] <Arlg@ - gL, (8.3.23) 


Vt,x,y€[a,b]:x>t>y; 1 >0) 
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So / (gx) — 9)" g (OG) —gM)at = 
(a) Jy 
ae G@-9@r'¢ wae GOI" samy 
r @) T (a) (a+ 1) 
We have proved that 
If) — f(y) — (Ar (FP) @, vy) 9 @) -— 9 ODI S 
_ atl 
Ar (g@)-—g9)) (8.3.25) 
LP (a) (a+ 1) 
Vx,y €[a,b]:x>y. 
(ii) case of y > x : We have that 
If @) -— Ff) -— Ai (P) @, vy) G@@) - 9g ON = (8.3.26) 


If — f @) — (Ai () Oy) GO) -—g @)I = 
nore (gy) —9 "7G (1) (Doig f) Ot 


(g(y)-g@))*] 
Tia+l | 


(es gi) (y) 


= 


| gy) -—9))* "9 © (Dig f) O - (D&i.gf) 0) a 
‘ (8.3.27) 


[oo-so JO || (Def) O — (Deg f) Oat (8.3.28) 
I (a) x righ x+ig a: 


I (a) 


(we assume here that 


| (Pgf) O - (Deygf) | <2lgO-IOI. (8.3.29) 


Vt,y,x € [a,b]: y>t>x;r.>0) 


<= 2 i gy) -9@)°' OGO)-—g@)dt= (8.3.30) 
(a) Jy 


2 (g(y) —g(x))°*! 
Tl (a) (a+ 1) 


2 f Gun _ 
| (gQ)-—g@))"°g dt = (8.3.31) 


We have proved that 
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If O-FO-Ae@yn@@)-gODI Ss (8.3.32) 


2 (g(y) —g(x))°*" | 
['@ teh? 7 eles 


Conclusion 8.6 Set \ := max (A), A2). Then 


If O-FfO)-AG)@ y@G@- gO Ss 


X Ig@- 9G)" 
I (a) (a+ 1) 


,Vx,y € [a,b]. (8.3.33) 


Notice that (8.3.33) is trivially true when x = y. 
One may assume that 


aS el, (8.3.34) 


Now based on (8.3.20) and (8.3.33), we can apply our numerical methods presented 
in this chapter to solve f (x) = 0. 

(I) In the next background again we use [5]. 

We need 


Definition 8.7 ((5]) Let a > 0, [a] = n, [-] the ceiling of the number. Let 
f € C" ({a,b], X), where [a,b] C R, and (X,||-||) is a Banach space. Let 
g € C! ({a, b]), strictly increasing, such that g~! € C” ([g (a), g (b)]). 

We define the right generalized g-fractional derivative X-valued of f of order a 
as follows: 


_] n b , 
(De, f) @) i= oo [ go-seorr'9 @ (Fog) Wena, 


(8.3.35) 
V x € [a, b]. The last integral is of Bochner type. 


If a ¢ N, by [5], we have that (D5..,f) € C ({a, b], X). 
We see that 


B(CKY" (Fog!) og) @) = (Dpyf)@.asxsb. 83.36) 
We set 
Dig f @) = (-1" ((Fog!)" og) @EC(ab], XEN, (8.3.37) 


Dj_.gf (x) = f (&), Vx € [a, B]. 


When g = id, then 
Dig f @) = Deaf (%) = Def, (8.3.38) 
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the usual right X-valued Caputo fractional derivative, see [3]. 
We also need the Taylor’s formula. 


Theorem 8.8 ({5]) Leta > 0,n = [a], and f € C" ([a, b], X), where [a,b] CR 
and (X, ||-||) is a Banach space. Let g € C! ({a, b]), strictly increasing, such that 
ec" ((g@.g())),a <x <b. Then 


n—1 


b 
f@= FO) +P GOI (f, og!) gb) + 


i=1 


i, 7? 7 ; 
Tia) —= | (g(t) — g(x))*"' 79! (1) (Dp_.g f) Oat = 


n—1 


by)! 
fOr y FO 8) (5. og!) (gb) + (8.3.39) 


i=1 


1 g(b) chi | ; 
nos g(x)" ((Dp_yf) 0g!) © dz. 


The remainder of (8.3.39) is a continuous function in x € [a, b]. 


Here we are going to operate more generally. We consider f € C” ([a, b], X). 
We define the following X-valued right g -fractional derivative of f of order a as 
follows: 


a by" 


(D°.,f) @) = a ff (90) —g@))" 14 @ (Ff og)” Godt, 


(8.3.40) 
Vx € [a, y]; where y € [a,b]; 


(DI f)@O=CD" (fog) og) @).¥x,yelabl, 8341) 
(Dy_.,f) @) = f @).¥x € [a,b]. (8.3.42) 

For a > 0, a ¢ N, by convention we set that 
(Do. f) (x) = 0, for x > y, Vx, y € [a, B]. (8.3.43) 


Similarly, we define 


a eal * n—a-l I -1\" 
(Deaf) (y) = nea (g(t) -—g()) 'd (fog 1) (g (t)) dt, 
(8.3.44) 
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V y € [a, x], where x € [a, b]; 
(Dif) 0) =D" ((Fog7!)” 09) O).V x.y € [a,b], 
(Do_.gf) O) = £0), Vy € (a,b. 
For a > 0, a ¢ N, by convention we set that 
(De. f) &) = 0, for y > x, Vx, y € [a,b]. 


We get that 


| (fo ey)" og _— a 
Tran-—atl) (g (b) — g (x)) < 


| (f ° gy og loan b nay 
Tan—at+) (g (b) — g(a))" *, Vx € [a,b]. 


|| (Dp_., f) @|| < 


That is 
(DP, f) @)=0, 


and 
(Pyof} (= (Oo ,7) (x) =0, Vx, y € [a,b]. 


Thus when a > 0, a ¢ N, both Dy gd D? at € C ([a, b], X), see [5]. 


met 


Hence by Theorem 8.8 we obtain 


n—-1 eS —1)\k) 
safe y E d wy 


k=1 


(g(x) — g(y))k + 


1 y 
mol (g(t)-g(x))*"'0'@) (DE) (t)dt, alla<x<y<b. 


Also, we have 


came ee a 
fore g ) (g(x) 


k=1 


(g(y) — g (x))K + 


1 x 
l(a) 


Pia) / G@-gW" gO (DL, f) Odt, alla<y <x <b. 
a 
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(8.3.45) 


(8.3.46) 


(8.3.47) 


(8.3.48) 


(8.3.49) 


(8.3.50) 


(8.3.51) 


(8.3.52) 
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We define also the following X-valued linear operator 


(Az (f)) @, Y) = 


n—-l 1) 


° y _ a-1 
YEE GG) — ao" = (Dh gh) 0) EET, fore < y, 


) 1 
b (fog” zh (g(x) gi) —-g (x))* 1 — (pe 7) (y) eer forx > y, 
fo (x), when x = y, 


(8.3.53) 
Vx,y €[a,bl;a>0,n= fa]. 


We may assume that ([13], p. 3) 


(Az (f)) x) — (A2 (fF), I = LF @ - FO) (8.3.54) 
< ®* |g(x)-g9Q)|, Vx, y € [a,b]; 


where ®* > 0. 
We estimate and have 
(i) case of x < y: 


If @) = FO) = 2 (A) & 9) G@) = ODI = 
a-l a 
noe (GO) — 9 CN GO (Dig f) O4t— 


(9 (y) — g (x))* 


(Deaf) (x) T(at+l 


(8.3.55) 


(by [1] p. 426, Theorem 11.43) 


moll (g (t) — g(x)! 7’ @) ((DE_., f) © — (DS. f) @)) at 


(by [9)) 


1 J a-l vs 
<ao/ (g(t) — 9 (x))*"'79' 0) | (DS. f) © — (DS_.gf) || dt (8.3.56) 


(we assume that 


| (Def) © -— (Dif) @] < lg -—9 I, (8.3.57) 


Vt,x,y € [a,b]: y=>t=x; p, > 0) 
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g(x)" 70 () (gt) —g@))dt= 


f° wea — pr Gy) —g@))*! 
7 a f (g@) —g(x))"g (dt = Pia) +) (8.3.58) 
We have proved that 
If @) — f(y) — (A2 (fF) @ y) (9 @) — 9 OI S 
a+l 
pi (g(y) — 9 @)) (8.3.59) 
T (a) (a+ 1) 
Vx,y ela, b]:x<y. 
(ii) case of x > y: 
If @®) -— Ff) -— (A2(f)) @, ») G@) - gO) = 
If — f @) — (277) @ WG) -—9g@)Il= (8.3.60) 


If (y) — fF @&) + (A2 (f)) OY) G@) -— gO) = 


lx = [ (9) — 90)! g' (Df) Dat- 
(g(x)-g ("| 


(De nf) rere ae J 


< 


raf eo-son' 70 (or gf) () — (De oink. 
(8.3.61) 


1; GH) -gO) I O(DLgf) O- (Deigf) OI] at (8.3.62) 


(we assume that 


| (Du, f) © - (DX, f) Ol < 21g @ —9 OI. (8.3.63) 


Vt,y,x €[a,b]:x >t> y; po. >0) 


< oF / g@)-—90))*'!¢ OGWM-gQ))dt= 
(a) Jy 


Tr a al (g(t) —g9())* 9 (at = (8.3.64) 
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pr (g(x)—g9 oy 


Pia) a0 (8.3.65) 
We have proved that 
If O- FO) - A277) @ y@@)-gODI S 
pr (g&)-gQ))" 
F(a) (a+1) ,Vx,yela,blix>y. (8.3.66) 
Conclusion 8.9 Set p := max (pj, p2). Then 
If @) — fy) — 42) @ GO) - 9g ODI S 
p_\g@)—g(yl*" 
F@) Gat , Vx,y € [a,b]. (8.3.67) 
Notice that (8.3.67) is trivially true when x = y. 
One may assume that 
tei, (8.3.68) 
T (a) 


Now based on (8.3.54) and (8.3.67), we can apply our numerical methods presented 
in this chapter to solve f (x) = 0. 

In both fractional applications a + 1 > 2, iffa > 1. 

Also some examples for g follow: 


g(x) =e*,x € [a,b] CR, 

g(x) = sinx, 

g(x) = tanx, 

where x € [-$ +e, 7 _ | , where ¢ > O small. 


(8.3.69) 
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Chapter 9 
Generating Sequences for Solving in Abstract 
g-Fractional Calculus 


The aim of this chapter is to utilize proper iterative methods for solving equations on 
Banach spaces. The differentiability of the operator involved is not assumed neither 
the convexity of its domain. Applications of the semi-local convergence are suggested 
including Banach space valued functions of fractional calculus, where all integrals 
are of Bochner-type. It follows [6]. 


9.1 Introduction 


Sections 9.1—9.2 are prerequisites for Sect. 9.3. 

Let B,, By stand for Banach spaces and let Q stand for an open subset of By. 
Let also U (z, p) := {u € B, : ||u — z|| < p} and let U (z, p) stand for the closure of 
U (z, p). 

Many problems in Computational Sciences, Engineering, Mathematical Chem- 
istry, Mathematical Physics, Mathematical Economics and other disciplines can writ- 
ten as 

F(x) =0 (9.1.1) 


using Mathematical Modeling [1-18], where F : 82 — Bp is a continuous operator. 
The solution x* of Eq. (9.1.1) is sought in closed form, but this is attainable only 
in special cases. That explains why most solution methods for such equations are 
usually iterative. There is a plethora of iterative methods for solving Eq. (9.1.1), more 
the [2, 7, 8, 10-14, 16, 17]. 

Newton’s method [7, 8, 12, 16, 17]: 


Kat = Xn — F’ (tn) * F in). (9.1.2) 

Secant method: : 
Xn4+1 = Xn — [%n—15 Xn3 Fl F (Xn) ; (9.1.3) 
© Springer International Publishing AG 2018 139 
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where [-, -; F'] denotes a divided difference of order one on Q x  [8, 16, 17]. 
Newton-like method: 
Xnt1 = Xn — Ey F (x), (9.1.4) 


where E, = E(F)(x,) and E : Q — L(B,, Bo) the space of bounded linear 
operators from B, into Bz. Other methods can be found in [8, 12, 16, 17] and the 
references therein. 
In the present study we consider the new method defined for each n = 0, 1, 2, ... 
by 
Xn+1 = G (Xn) 


G (Xn41) = G (tn) — Ay F Gn), (9.1.5) 


where xo € Q is an initial point, G : Bz; — Q (B3 a Banach space), A, = 
A(F) (n41,%n) = A (Xn41,%) and A: Q x Q > L(B), By). Method (9.1.5) 
generates a sequence which we shall show converges to x* under some Lipschitz- 
type conditions (to be precised in Sect. 9.2). Although method (9.1.5) (and Sect. 9.2) 
is of independent interest, it is nevertheless designed especially to be used in g- 
Abstract Fractional Calculus (to be precised in Sect.9.3). As far as we know such 
iterative methods have not yet appeared in connection to solve equations in Abstract 
Fractional Calculus. 

In this chapter we present the semi-local convergence of method (9.1.5) in 
Sect.9.2. Some applications to Abstract g-Fractional Calculus are suggested in 
Sect.9.3 on a certain Banach space valued functions, where all the integrals are 
of Bochner-type [9, 15]. 


9.2 Semi-local Convergence Analysis 


We present the semi-local convergence analysis of method (9.1.5) using conditions 
(M): 
(m,) F : Q Cc B, — Bp» is continuous, G : B; — Q is continuous and 
A(x, y) € £(B), Bo) for each (x, y) € Q x Q. 
(m) There exist 3 > 0 and Qo C B; such that A (x, y)~! € L (Bo, B;) for each 
(x, y) € Qo x Qo and 
JAG yy" <8". 


Set Qy =QN Qo. 


(m3) There exists a continuous and nondecreasing function w : [0, +00)? > 
[0, +00) such that for each x, y € Q) 


IF @)-FO)-AG, y)(G@)-GO))I < 


Bw (lx — yl. lle — xoll ly — xoll) IG @) -— GOI. 
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(m4) There exists a continuous and nondecreasing function wo : [0, +00) > 
[0, +00) such that for each x € Qy 


|G (x) — G (ro)Il_ S Yo (lx — xoll) lx — xoll - 
(ms) For x9 € Qo and x; = G (xq) € Qo there exists 7 > 0 such that 
||A G1, x0)! F (xo) || < n- 
(m6) There exists s > 0 such that 
v(m, s,s) <1, 
Yo (8) < 1 


and 


(Go) =x <s2—"— 
— 40 


where go = w (7), 5, 8). 

(m7) U (x0, 5) CQ. 

Next, we present the semi-local convergence analysis for method (9.1.5) using 
the conditions (/) and the preceding notation. 


Theorem 9.1 Assume that the conditions (M) hold. Then, sequence {x,} generated 
by method (9.1.5) starting at xy € (2 is well defined in U (xq, 5), remains in U (xo, S) 
for each n = 0,1,2,... and converges to a solution x* € U (xo, 5) of equation 
F (x) = 0. The limit point x* is the unique solution of equation F (x) = 0 in 
U (xo, S) e 

Proof By the definition of s and (ms), we have x; € U (xo, s). The proof is based on 
mathematical induction on k. Suppose that ||, — x,—1|| < ge 'n and ||x, — xol| <s. 

We get by (9.1.5), (m2) — (ms) in turn that 


IG Gays) — G Oa) = Ap’ F Gd = 
| Ar! CF Gee) — F e-1) — Ani (Ge) — G Gr-1))) | 
< Ag! || IF Gd) — F Gai) — Ani (G Ge) — G Ge—))I 
8-8 (late — Xe—1ll  llte—1 — oll s Ilye — xoll) IG Ge) — G @e-DII S 


(nS, 8) IG xe) — G (x1) Il = Go IG Ge) — E (Oe_1) II < 9K Ilan — x0ll < ahn 
(9.2.1) 
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and by (m6) 
Ilxx+1 — Xoll = IG (xe) — Xoll S IG Gx) — G Co) Il + IG (0) — oll 
S Yo (lle — xoll) Ilex — xoll + IG Go) — xoll 
< Yo (s) 5 + [IG (xo) — xoll < ». 


The induction is completed. Moreover, we have by (9.2.1) that for m = 0, 1, 2, ... 


Ixeem — aell $< gh, 

— 90 

It follows from the preceding inequation that sequence {G (x;,)} is complete in a 

Banach space B, and as such it converges to some x* € U (xo, 5) (since U (xo, 5) 

is a closed ball). By letting k — +00 in (9.2.1) we get F (x*) = 0. We also get 

by (9.1.5) that G (x*) = x*. To show the uniqueness part, let x** € U (xo, 5) bea 

solution of equation F (x) = 0 and G (x**) = x**. By using (9.1.5), we obtain in 
turn that 


\|x** — G Gea) || = |x — G @) + Ap! F Gx) — Ap! F (x) | < 
4c" |F (e) — F Gx) — At (6) - Gen)|| s 


B" Babo (x* — xe]. lea — oll. Ilex — xoll) |G (x) — G @)I| < 


’ 


qo |G (x) — G Gx) S40" x — x0 


so lim x, = x**. We have shown that lim x, = x*, sox* = x**. | 
k—>+00 k—>+00 
Remark 9.2 (1) Condition (mz) can become part of condition (m3) by considering 


(m3) There exists a continuous and nondecreasing function ¢ : [0, +00)? > 
[0, +00) such that for each x, y € Q 


|A (x, y) | LF @&) — FQ) — A(x, y) (E(x) — Gy) < 
yg (IIx — yl], lx — xoll, lly — xoll) IG @) — GG)I. 


Notice that 
yp (U1, U2, U3) < W (U4, U2, U3) 


for each uy > 0, ua > O and uz > 0. Similarly, a function y; can replace 7 
for the uniqueness of the solution part. These replacements are of Mysovskii-type 
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[7, 12, 16] and influence the weaking of the convergence criterion in (776), error 
bounds and the precision of s. 
(2) Suppose that there exist @ > 0, 3; > O and L € L(By, Bo) with L~! € 
L (B2, B,) such that 
Je“ | <8" 


|A@,y)-Llil sf 


and 


(p= 8B =. 


Then, it follows from the Banach lemma on invertible operators [12], and 


|Z IAG@ »- Lis 8A = <1 
that A (x, y)~! € £ (Bo, B,). Let B = —_ Then, under these replacements, condi- 
tion (m2) is implied, therefore it can be dropped from the conditions (M). 


Remark 9.3 Sect.9.2 has an interest independent of Sect. 9.3. Itis worth noticing that 
the results especially of Theorem 9.1 can apply in Abstract g-Fractional Calculus as 
illustrated in Sect.9.3. By specializing function 7, we can apply the results of say 
Theorem 9.1 in the examples suggested in Sect.9.3. In particular for (9.3.21), we 
choose for u; > 0, u2 > 0, u3 > 0 


AT 


W (uy, U2, U3) = AP Ww+ 1’ 


if |g (x) — g (y)| < 4 for each x, y € [a, Db]; 


Au3 


W (uy, U2, U3) = BPO Wt’ 


if |g (x) — 9g (y)| < &2 |x — y|| for each x, y € [a, b] and pz = & |b — al; 


AUS 


WY (uy, U2, U3) = Ar wywW+1)’ 


if |g (x)| < & for each x, y € [a, b] and 3 = 263, where X, v and F are defined in 
Sect. 9.3. Other choices of function w are also possible. 

Notice that with these choices of function 7 and f = F and g = G, crucial 
condition (m3) is satisfied, which justifies our definition of method (9.1.5). We can 
provide similar choices for the other examples of Sect. 9.3. 
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9.3. Applications to X-valued g-Fractional Calculus 
of Canavati Type 


Here we deal with Banach space (X, ||-||) valued functions f of real domain [a, D]. 
All integrals here are of Bochner-type, see [15]. The derivatives of f are defined 
similarly to numerical ones, see [18], pp. 83-86 and p. 93. 

Here both needed backgrounds come from [5]. 

Letv > 1, v € N, with integral part [v] = n € N. Let g: [a,b] ~ Rbea 
strictly increasing function, such that g € C! ({a, d)), g°| €C" ([g(@, 9 (b))), and 
let f € C” ({a, b] , X). It clear then we obtain that (f fo) g") E 
C" ([g (a), g(b)], X). Leta:=v—-[v)=v—-—n(O<a< 1). 

(1) See [5]. Leth € C ([g (a), g (b)], X), we define the X-valued left Riemann- 
Liouville fractional integral as 


(JPh) @) = a (g—t)""h(t)dt, (9.3.1) 


rv) 


for g (a) < zo < z < g (b), where I is the gamma function. 
We define the subspace Coa) (Lg (a), g (b)], X) of C” (Lg (a), g (b)], X), where 
x € [a,b]: 
C%y (Ig @),. 9), X) = 


{heC™ Ug @).9@)].X): IQA ECLV@).g@I.X}. 03.2) 


So leth € C 90) (Lg (a), g (b)], X); we define the X-valued left g -generalized 
fractional derivative of h of order v, of Canavati type, over [g (x) , g (b)] as 


D'oh := = (En) (9.3.3) 


g(x) 


Clearly, for h € C” 


vex) (Lg (a), g (b)], X), there exists 


Vv — 1 a : _ 4)7@ p(n) 
(Diayh) @ = Tdow a y. (2—t) *h™ (t) dt, (9.3.4) 


forall g(x) <z<g(b). 
In particular, when f o g7! € C 108 (Lg (a) , g (b)], X) we have that 


Zz 
& 


: ¢ =a —1)(n) 
ce (fog™) (dt, (9.3.5) 


(Diy (Fog) @= 
for allz: g(x) <z<g(b). 


We have that D?.,. (fog™') =(fo gy and Do... (fog!) = fog”. 
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From [5] we have for (f co) g') E Coa) (Lg (a), g (b)], X), where x € [a, b], 
(X-valued left fractional Taylor’s formula) that 


n—1 S = (k) 
fo-fe@=rl : J 2 6Gh26ys (9.3.6) 


k=1 


1 a v-1 V —l 
To o (gy) — 1)" (Diy (Ff o9 )) @ dt, forall y € [a,b]: y>x. 


Alternatively, for (f og!) € Cry (Lg (a), g (b)], X), where y € [a, b], we can 
write (again X-valued left fractional Taylor’s formula) that: 


n-1 og7! (k) 
fo-7os> <= J 8) ge) —g Oy + (9.3.7) 


k=1 


1 g(x) 
PY) Joc) 


(g (x) —t)”"! (DY, (fog™')) ® dt, forall x € [a,b]: x > y. 


Here we consider f € C” ({a,b], X), such that (fog') € Crea (Lg (a) , 
g (b)], X), for every x € [a,b]; which is the same as (f ) g') E Coy Ig (a), 
g(b)], X), for every y € [a, b] (i.e. exchange roles of x and y); we write that as 


(fog ') € Cr, (g(a), 9 (b)], X). 
We have that 


(Diy (F297) @= Fara a 


1 d [* = ~1)(@) 
/ Gn * (fog) (dt, 9.3.8) 
gy) 
for allz: g(y) <z<g(b). 
So here we work with f € C” ([a, b], X), such that (f ° g') €E Cie (Lg(a) , 
g (b)], X). 
We define the X-valued left linear fractional operator 


~1\ (4) 
4 (fe (x _ 
sitet Foo) 9 (yy — g (yy! + 
yo vl 
(Dy. (f 0971) (GO) LR y > x, 


(A (f)) @&, y) = ae (fea) "@o) CC ae (9.3.9) 


= yyy-! 
(Do (fo c')) (9 (x) maser ae x>y, 


if" Qi.2ey 
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We may assume that (see [13], p. 3) 
Ar (A) G2) = (A (PO YUL = FC) -— F OD = 
IPP or) Ga)-(FP og) GO| <Flg@—-gO)I, 3-10) 
where ® > 0; for any x, y € [a, b]. 
We make the following estimations: 


(i) case of y > x : We have that 


lf &) - f &) — (Ai (f)) &@, vy) gy) -— g (x) = 
1 gy) a ? 2 
foe (9) — 9° (Dien (F097) Oat 


(g(y) — 9 (x))” 


(Di (fog ')) GO) Tort 


(by [1], p. 426, Theorem 11.43) 


[ (g(y)—1)”"! ((Px5 (fo ')) (t) — (Dy (fo c')) iC ))) a 


. (9.3.11) 
(by [9]) 
=F a 90) = 9°" | (Diy (F297) © = (Phen (Fe97')) @ | at 


(we assume here that 
| (Dhoo (F 097) © = (Dhey (F297) GO] <ALe= GOI, 3.12) 


for every t, g (y), 9 (x) € [g (a), g (b)] such that g (y) > t = g(x); A; > 9) 


dM" gy) i 
= (g(y)— 1)" (gy) —tdt = (9.3.13) 
L(Y) Jax) 
A pe Ar (g(y) = g @))* 
—— —t)’dt= : 9.3.14 
ie TO CUS ean 
We have proved that 


A (g(y) — 9 (@))"t! 
Ifo -f@- (Ai (f/f) @ yg) -—9g@)Il< rw Oat ; 
(9.3.15) 
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forall x, y € [a,b]: y>x. 
(ii) Case of x > y : We observe that 


Ifo) — f &) -— (Arf) & 9) GO) — 9g @)I = 
If a) - fO)-(1(f) &. GO) - gO) = 


(g (x) — 1)" (Dj (F 0. 97')) Mat— 


| 1 g(x) 
PW) Ig) 


(g (x) — 9 (y))” 


V =] = 
(Diy (F29"')) GO) = a i= (9.3.16) 
1 [~ a@ij—1y-* (7 (fog!) (t) — (Dy (fog!) (g(x))) dt 
PY) to) iY) 9) 
1 gi) ae 7 a4 . ] 
=F Ih) 279 [Pion (FD) © ~ (Pion (F271) oon] a 


(9.3.17) 
(we assume that 


[Poo F297) O— Dgoy F097) GON] Ss Arle —G@I, 0.3.18) 


for all t, g (x), 9 (y) € [9 (a), g (B)] such that g (x) > t = g(y); A2 > 0) 


No g(x) 
= (g(x) —1)""! (g(@) 1) dt = (9.3.19) 
PY) Jay) 
dy 9 ya, 2 YG@-gQ))"! 
TW i WGN a= a ace ty 
We have proved that 


2 (g(x) — 9)" 
If) — f @) -— (Ar (f)) @ y) GO) — 9 &))Il < rw) w+) 


(9.3.20) 
for any x,y € [a,b]: x>y. 


Conclusion 9.4 Set \ := max (1, A2). Then 


X |gQ)-g@|""" 
If) — f &) -— (Ai (f)) @, y) gy) — 9g &)I S rw wi , 
(9.3.21) 


V x,y € [a, b] (the case of x = y is trivially true). 
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We may choose that _ <1. 
Also we notice here that yy + 1 > 2. 
(II) See [5] again. Leth € C ([g (a), g(b)], X), we define the X-valued right 


Riemann-Liouville fractional integral as 
(UP h) (z) := : [ (t —z)’ | h(t) dt (9.3.22) 
zo—P) (2) = 5 Qi Zz ; 3) 


for g(a) <z< 2% <g(b). 
We define the subspace CGA ([g (a), g (b)], X) of C” ([g (a), g (b)], X), where 
x € [a,b]: 
Cry (Lg a). 9 (b)], X) = 


fh eC™ (ga). g@1.X): Ij _h EC (Ig@).g@)1,X)}. 0.3.23) 


So let h € Ca (Lg (a) , g (b)], X); we define the X-valued right g -generalized 


fractional derivative of h of order v, of Canavati type, over [g (a) , g (x)] as 


Dif <= (-1"! (ga): (9.3.24) 


Clearly, for h € Cee (Lg (a) , g (b)], X), there exists 


(-1)""! d sy spas 
(Dio h ies Pla) de (t—z) “h” @)dt, (9.3.25) 


forall g(a) <z< g(x) < g(). 
In particular, when f og7! € Cia (Lg (a), g (b)], X) we have that 


=| n—-1 d g(x) ‘ 
Oho FeO VO= Tamaya fl CHI Fer)” eode 


ld -—a) dz 
(9.3.26) 
for all g(a) <z< g(x) <g(b). 
We get that 
(Dio (fog ) @=EN" (fog), (9.3.27) 
and 
(Dyw- (fog) @ = (fog), (9.3.28) 


for all z € [g (a), g (x)], see [5]. 
From [5] we have, for (fog ) € Cag (Lg (a), g (b)], X), where x € [a, b], 


vy > | (X-valued right Pan orl Taylor’s formula) that: 
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ye: g') Go) 


ai (g(y) — g(x) + 


iQl=f R= 


k=1 


g(x) 
os (¢-—g9())”" (Diy (f09")) Oat, dlasy<x. (9.3.29) 
gy 


Alternatively, for (f 0 g~') € Choy (Lg (@) , 9 (b)], X), where y € [a, b],v > 1 


(again X-valued right eal Taylor’s formula) that: 


3 (Fog ')® GW) 


7 (g(x) — gy) + 


fa@-fO= 


k=1 


gy) 
rol... (t —g(x))”* (Di, (F og')) Oat, allasx<y. (9.3.30) 
g(x 


Here we consider f € C” ([a,b], X), such that (f fe) g') € Coay_(g (@), 
g (b)], X), for every x € [a,b]; which is the same as ta fe) ge) ec’ 26)— (Lg (a), 
g (b)], X), for every y € [a,b]; (i.e. exchange roles of x and y) we write that as 


(fog') eC? (lg@),g(b)], X). 
We have that 


—1)" 1 d gy) : 
Cin KO NOTE, CI (For)” wade 
; (9.3.31) 


for all g(a) <z <g(y) <g(b). 
So here we work with f € C” ({a, b], X), such that (f ° g*) € Co_ (g(a), 
g (b)], X). 


We define the X-valued right linear fractional operator 


og! = 
— ciao Ge) G6ysae= 
—atyyyrv—1 
(Dia (fog =) (9 (y)) oe x>y, 


(A2 (f)) 9) = | act (fea) oo» Gay-gO> = (9.3.32) 


(g)-9))" 
(Dy, (fog) Ge) LB, y > x, 


i? @) 22% 


We may assume that ([13], p. 3) 


(Az (f)) @&, x) — (A2 (f)) OY, I = | F @ -— fF O || < O*lg@ -—g OIL, 
(9.3.33) 
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where ®* > 0; for any x, y € [a,b]. 
We make the following estimations: 
(i) case of x > y : We have that 


If @)—- f &)- (42°77) & yg) -gO)I = 
If ) — f &) — (A2 (f)) & GO) -— 9g @)I = 


(9.3.34) 
If (y) — f @&) + (A2 (f)) & 9) (G @) -— 9g ODI = 
1 g(x) 4 7 “i 
| TD Jos (t-— 90)" (Diy (F 0.97’) Oat—- 
v zs (g@)-—g9())”’ 
(Di (f og ')) (g(y)) “Toth (9.3.35) 


(by [1], p. 426, Theorem 11.43) 


1 
~ Tv) 


g(x) 
[ @— 9)" ((Dgey- (F 9.97')) O — (Dgw- (F 2 97'))) 
gy. 


(g (y)) dtl| 
(by [9)) 


= 


[2 e209" (Bian (F209) ©~ (Dhay- or) on a 


(9.3.36) 
(we assume here that 


| (Dx (F097) O = (Diy (Fe97')) GO) < prlt-—gO)], 0.3.37) 
for every t, g (vy), g (x) € [g (a) , g (b)] such that g (x) > t = g(y); pi > 9) 


i / . (t—g(y))” "(tt -gQ)) dt = 
rv) gy) on “— - 


= 


pi g(x) 


pr (g(x) -— gy)" 
= “dt= 9.3.38 
FO) Joo (t — g(y))" at ro w+) ( ) 
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We have proved that 


pr (g(x) -—g (yt! 


If) — Ff) — A2) @ IG @)-— I ODI S 


~ ry) (v+1) 
(9.3.39) 
Vx,yela,b]:x>y. 
(ii) Case of x < y : We have that 
If &) — f O) — (A2 (f)) @, y) (9 @) — 9 OY) = 
If &) — f(y) + (A2 ()) @, (9) — 9 @))I = (9.3.40) 


1 gy) : - 
lFo i (9 @))"" (Dyyy- (fF 097!) Oat— 


(gy) — 9 @))” 
Tw+l) 


(Di (fog-')) G@) 


1 
r(v) 


i G-7o)y"" ((Dyo)- (fo o')) (t) — (Diy) (fo 9')) (g (x))) dt 


— — pr (t-—g(x))”! | (Diy) (fo o')) (t) — (Diy (fo o')) é («| a 
(9.3.41) 
(we assume that 


[Poo Fog") © = (Diw-(F 297") GOD) |< p2lt— 9 @)I, 0.3.42) 


for any t,g(x),9() €lg@, 9g): 90) =t = 9); po > 9) 


< ef gy ¢ gat = 
rae g(x g(x = 
p2 g(y) : _ 
rw) [. (t—g(x))" dt = (9.3.43) 


p2 (g(y)—g(x))”*! 


rv) (y+ 1) 


(9.3.44) 


We have proved that 


p2_ (gly) — 9 @)) 
lf) -—fO)- (2(f) & y G@)-gO))I < ro w+) 


(9.3.45) 
Vx,y ela, b]:ix<y. 
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Conclusion 9.5 Set p := max (p1, p2). Then 


y+l1 


p |g(x)-g()I 
lf @) — f &) — (42 (P)) @, y») 9G @) -— 9 ODI S TW) w+ 


’ 


(9.3.46) 
V x,y € [a, b] (9.3.46) is trivially true when x = y). 


Pp 
One may choose To < 1, 


Here again y+ 1 > 2. 


Conclusion 9.6 Based on (9.3.10) and (9.3.21) of (I), and based on (9.3.33) and 
(9.3.46) of (II), using our numerical results presented earlier, we can solve numeri- 


cally f (x) =0. 


Some examples for g follow: 


g(x) =e*,x € [a,b] CR, 
g(x) =sinx, 
g (x) = tanx, 


where x € [-5 Pays e| , with e > 0 small. 
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Chapter 10 
Numerical Optimization and Fractional 
Invexity 


We present some proximal methods with invexity results involving fractional calcu- 
lus. It follows [3, 4]. 


10.1 Introduction 


We are concerned with the solution of the optimization problem defined by 


min F(x*) (10.1.1) 


s.t, x*ED 


where F : D C R” —-» Ris a convex mapping and D is an open and convex set. 
We shall study the convergence of the proximal point method for solving problem 
(10.1.1) defined by 


Xpi1 = argmin [Fay rs 4d (in. x*)| (10.1.2) 


x*eD 


where xo € D is an initial point, y > 0 and d is the distance on D. 

The rest of the chapter is organized as follows. In Sect. 10.2 we present the con- 
vergence of method (10.1.2) and in Sect. 10.3 we present the related application of 
the method using fractional derivatives. 


10.2 Convergence of Method (10.1.2) 


We need an auxiliary result about convex functions. 
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Lemma 10.1 Let Do © D be an open convex set, F : D —> Rand x* € D. 
Suppose that F + 4d?(.,.x*) : D — R is convex on Do. Then, mapping F is 
locally Lipschitz on Do. 


Proof By hypothesis F + 44°(., x*) is convex, so there exist L;,7; > O such that 
for each u, v € U(x", 1r)) 


Fw) 4 7d(u,3*) ~ (FW + Adv, x)| < Lidtu, 0), (10.2.1) 


x") 


It is well known that the mapping ab is strongly convex [9]. That is there exist 
Ly, r2 > 0 such that for each u, v € U(x*, 12) 


* d2( *) a _) 
= u,x)— = ee 4 
2 2 


< Lod u, v). (10.2.2) 


Let 
r=min{r},™} and Lo =L,+yL. (10.2.3) 


Then, using (10.2.1)—(10.2.3), we get in turn that 
Y 42 x Y 72 x 
Fu) — FQ)| = [Fw + 2d, x) = (FO) + 30, 2) 
4, Faw, x") = 2a(v, x") 
< Lyd(u, ¥) + Loyd (u,v) = Lod(u, v): (10.2.4) 


a 
Next, we present the main convergence result for method (10.1.2). 


Theorem 10.2 Under the hypotheses of Lemma 10.1, further suppose: 


— oo < inf F(x*), (10.2.5) 
x*eD 
Sy = {x* € D: F(x*) < F(y)} CD, inf F(x*) < FQ), (10.2.6) 
x*eD 


the minimizer set of F is non-empty, i.e. 


Paix 2 FO) = inf Fw} #9, (10.2.7) 
x*eD 
| F(x*) — x*|| < Ls, (10.2.8) 


f= U4 Sate 21, (10.2.9) 
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Then, the sequence {Xn} generated for x9 € S* := Sy 1 U(x*,r*) is well defined, 
remains in S* and converges to a point x** € T, where 


L3 


= : 10.2.10 
r i-L ( ) 


Proof Define the operator 
: ai * 
G(x) := F(x) + 5 |x —x*|. (10.2.11) 
We shall show that operator G is a contraction on U (x*, r*). Clearly sequence {x,} 
is well defined and since xo € S, we get that {x,} C S, foreachn = 0,1,2,.... 


In view of Lemma 10.1 and the definitions (10.2.8)-(10.2.11) we have in turn for 
u,v € U(x*,r*) 


1 1 
IG(u) — G0v)| < |F@) — FQv)| +7 50 ux") - 52 wx) 


< (Lo + yL2)dtu, v) = Ld(u, v) (10.2.12) 
and 


|G(u) — x*| < |G) — Ga*)| + |G@*) — x*| 
< Ld(u, x*) + |F(x*) — x*| 


< Ld(u,x*)+L3 <r". (10.2.13) 


The result now follows from (10.2.9), (10.2.12), (10.2.13) and the contraction map- 
ping principle [1, 5-8]. a 


10.3. Multivariate Fractional Derivatives and Invexity 


Let X = II [a;, bj]. 
i=1 


1. Let 0 < a < 1, we consider the left Caputo fractional partial derivatives of f 
of order a : 


O° f (x) _ 1 fo po OF Gets ay os Bs Satie Si) 
~Td-a)Ja 


dt;, (10.3.1 
Ox? Ox; mnet) 


U 
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where x = (x1,...,X) € X,i = 1,2,. em and 2A ssstn)) E Les la;: bi). 
i=1,2,...n.Here I stands for gamma daniction: Note that 


FIO) 2 [mw “dt, 
Ox? =z —a) 
| Of (1, x2, seey Xi-1, +) Xi41,--- Xn) 4 
Oxi 00, (abi) 
EA | OF Gig Rica tig MENG oy MORI eee A) 
T2-a) Ox; 00, (a; ,bi) 
<O, 
for alli = 1,2,...,n. Therefore, YL) exist for alli = 1,2, . 


Ox a 

Now we jaeae ihe left fractional Gradient of F of order a, 0. <a<l1: 
Of (x*) Of (x*) 

Ox? = Oxe . 


Vt f(x) = ( (10.3.2) 


We replace in Definition 10.3.1 of Verma [9], Vf (x*) by Vi f(x*). 
2. Let 0 < a < 1, we consider the right Caputo fractional partial derivatives of f 


of order a: 
6° —1 bi fo) Xs 4 cag Lig aweg By 
fe) _ | (x) 2 OF ou me to) yan eas 
Ox? ) OX; 
where x = (x1,...,%9) € X,i = 1,2,...n and 2G) € 7. (a;, bj), 
i= 1,2,...n. Note that 
OP f(R)| 
i —t;) “dt; 
axe | — <I [ i" 
| OF Kit Nyse MAT ag Mids 3 Kn) 4 
Oxi 00, (a; ,b;) 
_ or Gr Oe ats st Re as vee Ma)) 
r(2-a) Ox; 00, (a; ,b;) 
< ©, (10.3.4) 
for alli = 1,2,...,n. Therefore, — © exist for all i = 1,2;. 


Now we sonsiden the right fractional Gradient of F of order Q, ie <a<l: 
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Ve f(x") = (Ee. _ af). (10.3.5) 


a 
Ox} Oxo 


We replace in Definition 10.3.1 of Verma [9], V f (x*) by Vif (x*). 
3. Define fork EN: Vi f = Vi... Vif, k— times composition of left fractional 


gradient, 1.e., 
ko * kau * 
vir=(° fo") 2 =), (10.3.6) 


a 0 
Ox} Oxe 


kee oH 
where 2 ao = 2 


ae ~ f, k—times composition of left partial fractional 
‘ : ka ¢ . . 
ot exist for alli = 1,2,...n 


derivative,i = 1,2,... 


4. Define fork e N: Vint = Vv, .Wy f, k— times composition of right frac- 
tional gradient, i.e., 


_ Oke * aka * 
” ahi fo), (10.3.7) 
x} x 
ake f(x) ae _ e “ae : : : 
where —2- = Oat * at, k—times aa of right partial fractional 


a 


exist for alli = 1,2,. 


derivative, i = 1,2,. 


5. Let a > 1, we eanitier: the left Caputo fF rebral partial denivatives - f of 
order a (fa] =m EN, [.] ceiling of the number [2]): 


OW f(x) _ [ Oe; — 4)" god Midge: -%n) a4 
Ox? a —a) ' as 
(10.3.8) 
i=1,2,...n. We set an) equal to the ordinary partial derivative g of CO) We 
assume that an 
aan a ar) oney Hn) € Lo (aj, Di) 
Ox! 
i.e., 
eo” 
| ea < © 
Ox 00, (a; bi) 
for alli = 1, 2,...,n. Note that 
O° f (x) (Xj = a;)"— Qa O” f 
| Ox? = Ox” 7 (XI te aygtts Mp) < Ow, 
x T'(m —a+1) x; eta 
(10.3.9) 
for alli = 1,2,...n. Therefore, ae exist for alli = 1,2,...n. Now we 


consider the left facsonal gradient of f of ordera,a>1: 
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O° f(x*) ~) 


oe 
Ox} Oxe 


Vv fin )= ( (10.3.10) 
We can replace in Definition 10.3.1 of Verma [9] (or similar invexity formula), 
V f (x*) by Vo" f(x"). 

6. Let a > 1, we consider the right Caputo fractional partial derivatives of f of 
order a ([a@] = m): 


ae air OF eRe ids incase 
f(x) = (-1) i: (x; — 7 i f(xy X2 x D at, 
Ox? Pim — a) Jy, Ox}" 

7 (10.3.11) 
i=1,2,...n.Weset ot = (-1)” es (where ve is the ordinary partial). We 
assume that gn 

cae wen yey oes Xp) © Leg (Gi, Dj) 
Ox} 
for alli = 1, 2,...,n. Note that 
O° f(x), Gi —x)"* OF 
| aI < op (Xda een ey ees Xn) < oO, 
Ox; Tim — a+ 1) || Ox’ seta Gs 
for alli = 1,2,...n. Therefore, ote) exist for alli = 1,2,...n. Now we 


consider the right fractional gradient of f of ordera,a>1: 


Vz f@%= (>. ee ae). (10.3.12) 


a a 
Ox} Oxe 


We can replace in Definition 10.3.1 of Verma [9] (or similar invexity formula), 
V f (x*) by Vo 7 f (x*). 


Next we follow [4]. 
The discrete minmax fractional programming problem is 


(P) Minimize max fi @) 


: 10.3.13 
I<i<p gj (x) ( ) 


subject to G; (x) <0, jeg, A(x) =O0,k er, x eX, 
where X is an open convex subset of IR” (n-dimensional Euclidean space), fj, g;, 
ie p=({I,2,..., p}, Gj, 7 € gq, and AH, k € r, are real-valued functions defined 
on X, and for eachi € p, g; (x) > 0 for all x satisfying the constraints of (P). 
Consider a function f : X — R with fractional order derivatives, see [2]. In this 
section we give the definition based on the work [10] on several classes of generalizes 
convex functions. 
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Definition 10.3 ((4]) The function f is said to be (¢, 7, p, 8, m)-invex at x* for the 
left Caputo fractional partial derivative of order a, a > 1 if there exist functions 
@o:R—7>R,1: XxX > R",p: Xx X > R,andé: X x X — R", anda 
positive integer m such that for each x € X (x # x*), 


(FO) —F (x*)[) = (VET (a). 7") + 0 (2) A ( 2") |", 
(10.3.14) 
where 


eI are). (10.3.15) 


ae oe tn 
Vo f() = ( axe dies axe 


Definition 10.4 ((4]) The function f is said to be (¢, 7, p, 6, m)-invex at x* for the 
left Caputo fractional partial derivative of order a, 0 < a < | if there exist functions 
o:R—> R41: XxX > R",p: Xx X > R,andé: X x X — R’, anda 
positive integer m such that for each x € X (x 4 x*), 


(IF 0) =F (6°))) = (WEF (*) (0.2) 0 (x27) [8 (s.)", 10.3.16) 


Vif (x*) = (Geers ms aro). (10.3.17) 


a 
Oxy Oxe 


where 


Similar concepts hold for V7 f (x*) and V7 ~ f (x*). 
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